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Chapter 1 

Introduction 



In section 11.11 of this introduction, we place the subject of this thesis in 
the general context of statistical mechanics. In section II. 2| we discuss the 
motivations for our contribution to statistical mechanics in the emerging field 
of nanosciences. In section ri.3| we indicate the link between our approach 
of quantum statistical mechanics and the fundamental problem of quantum 
measurement. In section [m we expose the main goal of this thesis. Finally 
in section 11.51 we set out the plan of this thesis. 



1.1 Statistical mechanics 



An intuitive understanding of the natural laws occurs at the spatial scales, 
temporal scales and energy scales which are directly accessible to the sen- 
sory organs of the human being. These typical physical scales, which do 
not require a sophisticated experimental setup to make the measured in- 
formation accessible to us, were naturally the first historically investigated 
and broadly define what physicists call the macroscopic world. Thermody- 
namics is of course the most widely and best defined branch of macroscopic 
physics which also includes geometrical optics, fluid dynamics, mechanics 
and many others. Progress in the systematic study of nature came with the 
improvement of experimental techniques and has made possible the investi- 
gation of physical scales which are sometimes far less intuitive to the human 
being. For example, the investigation of very small spatial scales has led to 
the discovery of atoms and molecules while the investigation of very large 
spatial scales has opened the way to the discovery of solar systems, galaxies 
and galaxy clusters. Each of these new fields of investigation has given rise 
to physical laws and principles formulated in terms of the relevant quan- 
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titles i.e. the most easily measurable or understandable quantities within 
each level of description. In astrophysics, the dynamics of planets is, in first 
approximation, the dynamics of massive points interacting through Newto- 
nian gravitational forces. In particle physics, the interactions between the 
smallest matter constituents are often studied through the cross sections of 
colliding particles. 

The rational character of the scientific approach motivates a unified un- 
derstanding of these different physical laws occurring at various scales and 
gives a particularly important status to the physics of small scales, also 
called microscopic physics. In fact, if matter is made of fundamental ele- 
ments, say atoms and molecules, it is natural that the macroscopic physics 
describing this matter is the consequence of the microscopic physics de- 
scribing these fundamental elements. One should be able therefore to re- 
construct the macroscopic physics from the microscopic physics instead of 
building it on phenomenological arguments. This procedure allows a much 
deeper understanding of the macroscopic physics and the explicit derivation 
of its phenomenological elements. We know of course today that atoms and 
molecules are not the fundamental constituents of matter because they have 
an internal structure made of electrons and nucleons, the latter being again 
made of quarks, etc. However, if one can be sure that the internal structure 
of an atom will not contribute to some macroscopic phenomenon that one 
wishes to describe, the dynamics of the atom can be considered fundamental 
regarding this macroscopic phenomenon. This means that the microscopic 
theory that one uses depends on the nature of the macroscopic theory that 
one wishes to build. It is not clear at present if a definitive fundamental mi- 
croscopic theory (which cannot be reduced anymore to anything else) will 
ever be found. 

We can notice that an substential amount of information required for a 
microscopic description (for example all the positions and velocities of the 
atoms or molecules of a gas) is not relevant to understand a macroscopic 
phenomenon already captured by much fewer variables (such as the volume, 
pressure and temperature for an ideal gas). A huge reduction of the unnec- 
essary information has to be carried out to understand how the macroscopic 
world can emerge from the interaction between the microscopic constituents 
of matter. This is precisely the main goal of the branch of physics called 
statistical mechanics. It should be noticed that this body of knowledge is 
of great importance not only in physics but also in all the other fields in- 
vestigated by the human intelligence such as biology, economy and social 
behavior, where one wishes to understand the behavior of macroscopic vari- 
ables that consist in a large collection of microscopic variables. 
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In this thesis, we will focus on the specific part of statistical mechanics 
which deals with the physical systems described by non-relativistic quan- 
tum mechanics. As we know, the microscopic laws of nature are described 
by quantum mechanics. However, for many studies in statistical mechanics, 
one can suppose that the equations of motion of the microscopic constituents 
of a physical system are classical. This approach has a well-defined range 
of applicability (the de Broglie wavelength has to be shorter than the mean 
free path of the microscopic particle) but has the advantage of being more 
easily tractable. An important step to construct statistical mechanics is to 
take the so-called thermodynamic limit of a system. This means that one 
supposes that the number of microscopic variables in the system becomes 
very large so that the law of large numbers of probability theory applies. 
Classical statistical mechanics is concerned with systems which can be de- 
scribed in the classical limit before taking the thermodynamic limit. When 
the quantum effects manifest themselves, statistical mechanics has to use a 
quantum mechanical description. Quantum statistical mechanics therefore 
formally consists in taking the thermodynamic limit without any classical 
limit. 

The field of statistical mechanics can be roughly subdivided into equi- 
librium and nonequilibrium statistical mechanics. 

In equilibrium statistical mechanics all the thermodynamic fluxes in the sys- 
tem are vanishing and the macroscopic quantities describing the system are 
time independent. It is well known that the equilibrium probability distri- 
bution of the microscopic degrees of freedom of the system depends on the 
type of system one considers (microcanonical distribution for an isolated 
system, canonical distribution for a system exchanging energy with its en- 
vironment and grand-canonical distribution for a system exchanging energy 
and particles with its environment [2Q\ \^0\ I61j ) . 

In nonequilibrium statistical mechanics, there are three main ways of putting 
a system out of equilibrium. On can impose nonzero thermodynamic fluxes 
through the system using nonequilibrium boundary conditions (i.e. escape- 
rate formalism j29j^. This method is essentially used to study the nonequi- 
librium steady states of the system. One can also apply time-dependent ex- 
ternal perturbations to a system which induce local nonzero thermodynamic 
fluxes in this system (i.e. linear response theory [H^). This method is often 
used to study the response of a system to perturbations of its equilibrium 
state. Finally, one can choose an initial condition putting the system out 
of equilibrium (i.e. Boltzmann's kinetic theory jSI). This induces transient 
nonzero thermodynamic fluxes in the system which relax at long times to 
equilibrium. The nonequilibrium probability distribution of the microscopic 
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degrees of freedom of a system is more difficult to construct than its equilib- 
rium probability distribution and there is unfortunately no general method 
for this construction. Furthermore, the nonequilibrium distributions obeys 
kinetic equations which manifest irreversible properties. The origin of this 
irreversibility is an important problem at the core of nonequilibrium statis- 
tical mechanics. 

The work presented in this thesis deals with nonequilibrium quantum 
statistical mechanics. The kind of systems on which we focus our study of 
nonequilibrium quantum dynamics consist in a quantum subsystem which 
is in contact with a quantum environment. In most cases, the initial condi- 
tions of the total system composed of the subsystem plus the environment 
are initially out of equilibrium and one studies the relaxation of the total 
system toward equilibrium through the subsystem dynamics. Because of the 
interaction with the environment, the subsystem dynamics is then expected 
to become irreversible, relaxing toward some equilibrium distribution. The 
kinetic equations we will establish for the description of the irreversible dy- 
namics rely on the weak-coupling assumption between the subsystem and 
its environment. 



1.2 Emerging irreversible processes in nanosystems 

The fundamental assumptions of statistical mechanics arc based on proba- 
bility theory. They rely on the concept of statistical ensemble. A statistical 
ensemble is a collection of identical systems submitted to the same con- 
straints. However, the initial condition of each member of the ensemble is 
slightly different and a given probability is associated with each of these 
different initial conditions. The fundamental assumption necessary to make 
the statistical description of a system useful is that an individual member of 
the ensemble has, with a very high probability, the same behavior as the one 
averaged over all members of the statistical ensemble. The use of such argu- 
ments requires either a large number of degrees of freedom in the system in 
order to exploit the laws of large numbers or a chaotic dynamics which intro- 
duces the dynamical randomness necessary for the probabilistic description 
even if the system only possesses a few degrees of freedom. With the recent 
progress in experimental techniques, one can study in great detail single 
systems ranging from macroscopic scales down to the scale of the nanome- 
ter. This ability to study individual systems of size varying from one or a 
few atoms to thousand or more atoms allows us to investigate the crossover 
between a dynamical behavior very far from the one expected from kinetic 
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theory and a statistical behavior well described in terms of the concepts of 
nonequilibrium statistical mechanics. It is a great challenge to refine the 
statistical mechanical description in order to be able to study systems in 
these intermediate regimes. It is precisely this domain that nanotechnolo- 
gies are nowadays reaching. Notice that such systems are also very useful to 
investigate the importance of chaos (classical or quantum) for the emergence 
of statistical behavior in few-degree-of-freedom systems. 

We will now give some examples of typical systems in which such ques- 
tions can be studied. It is known that semiconducting devices can range from 
macroscopic scales (standard electronics) down to the scale of micrometer 
(microelectronics) or even to the nanoscopic scales (artificial atoms, quan- 
tum dots). These devices are characterized by their electronic conduction. 
Thanks to nanotechnology, it is nowadays possible to give different shapes 
to such electronic devices at the nanoscale |^ |^ EH EH] . Because they can 
have different shapes, giving rise to either classically integrable dynamics 
or classically chaotic dynamics, these quantum dots are ideal candidates to 
study the emerging statistical behaviors due to the dynamical complexity. 
In such systems one can also focus on the spin of the electron instead of its 
charge. This new emerging field studying the spin dynamics inside quantum 
dots is called spintronics jBHl EZl ,.2^ • The study of spin decoherence, 
which is a statistical phenomenon, in such devices is of particular impor- 
tance for potential future applications (i.e. quantum computers). Another 
example of quantum systems which can be studied today are single atoms 
in microwave cavities |76j . These spectroscopic experiments provide a de- 
tailed study of the atom-cavity entanglement and of the decoherence effect 
induced by the interaction of the atom-cavity with its surrounding. The fem- 
tochemistry study of IVR (intramolecular vibrational relaxation) in small 
molecules is also a very rich topic for statistical mechanics. Depending on 
the vibrational dynamical complexity (chaoticity) , the energy initially local- 
ized by a laser on a given chemical bond can flow or not in a statistical way 
through the other bonds and reach equipartition Besides, the different 
microscopy techniques (FIM, STM, AFM, . . . ) [331 142j play a fundamental 
role to allow the nanoscale investigations. Let us give some examples of 
nanometric systems in which emerging statistical behaviors can be studied 
by such techniques. FIM can allow the real time study with atomic res- 
olution of chemical reactions over catalytic surfaces which can give rise to 
nonlinear collective statistical phenomena such as surface wave propagation, 
concentration pattern formation and explosions ^Sj. STM allows the study 
of electronic transport properties of single molecules (nanotubes, conducting 
polymers and other molecules) [7511^ ITUn] . Besides, mechanical properties 
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of a single macromolecule can be studied by AFM in experiments which 
require statistical arguments for their interpretation |48j . 
These few examples from the nanoscience world show us the great number 
of possible applications of statistical-mechanical theories in this expanding 
field. There are numerous theoretical as well as technological challenges 
opened in the field of nanosciences that need to be investigated. 



1.3 Quantum measurement as a nonequilibrium 
process 

The measurement postulate of quantum mechanics, which tells us that the 
state of a measured system is projected with a given probability in one of 
the eigenstates of its measured observable, if applied to single quantum tra- 
jectories, is in apparent contradiction with the unitary evolution of quantum 
mechanics. It is natural, especially in the context of statistical mechanics, 
to consider that the measurement apparatus as well as the measured system 
can be described by the laws of quantum mechanics without introducing any 
dichotomy in the physical description between macroscopic and microscopic 
objects. Therefore, the unitary quantum evolution has to apply to the total 
system (measurement apparatus and measured system). The measurement 
process in this scheme automatically introduces a statistical description of 
the system dynamics. In fact, due to the non-separability of quantum me- 
chanics, even if the system and the measurement apparatus are described 
by pure states before the measurement, once the measurement process has 
started the system cannot be described anymore by a pure state but rather 
by a statistical mixture. This means that, according to the measurement 
postulate in its actual formulation, entropy is produced when a measure is 
performed on a quantum system. We see here an important analogy between 
the case of a subsystem interacting with its environment in nonequilibrium 
quantum statistical mechanics and the case of a measured system interact- 
ing with a measurement apparatus in quantum mechanics. This analogy 
already appears at the level of the total Hamiltonian, which is the sum of 
the subsystem Hamiltonian, the environment or the measurement apparatus 
Hamiltonian, and an interaction term describing the coupling between these 
two systems. Adopting this view of the problem, the delicate question to 
solve becomes to understand how the total unitary evolution can statisti- 
cally induce a measurement-like effect on the dynamics of the subsystem. 
We want here to point out the importance that nonequilibrium quantum 
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statistical mechanics has for the fundamental understanding of quantum 
mechanics and its measurement postulate. 

1.4 Goal of the thesis 

The general goal of this thesis is to contribute to the theoretical description 
of irreversible processes occurring in quantum systems. We would like more 
specifically to investigate three major questions. 

The first one concerns the application of the standard description of ir- 
reversible quantum processes to nanosystems. In fact, in the standard ap- 
proaches, one studies the dynamics of a small quantum subsystem in contact 
with an infinitely large environment ^. However, in nanosystems, the en- 
vironment cannot always be considered infinite. We therefore would like 
to elaborate a new theory permitting the description of the irreversible dy- 
namics of a subsystem interacting with a finite environment by taking into 
account the mutual effects of the environment on the subsystem and vice- 
versa. The statistical ensemble relevant to study such problems is no longer 
the canonical ensemble as for infinite environments, but the microcanonical 
ensemble. 

The second important question we want to answer in this thesis is related 
to the necessity of having a system with a continuous spectrum to manifest 
true irreversible quantum dynamics. We would like to understand how and 
under which conditions an almost-kinetic process can take place in a system 
with a discrete but dense spectrum. Understanding the transition between 
purely reversible quantum dynamics and kinetic quantum dynamics is also 
an important issue in nanosciences. 

Finally, the third major question to be solved in this thesis concerns quantum 
diffusion in spatially extended systems. We would like to understand how, 
and under which conditions, a diffusive transport process can take place, due 
to the interaction with the environment, in spatially extended quantum sys- 
tems. This problem of quantum diffusion is once again of great importance 
in nanosciences where transport processes can be studied in many different 
spatially extended nanosystems. 

1.5 Plan of the thesis 

In this thesis, after having introduced the basic concepts of quantum statis- 
tical mechanics (in chapter I^J and after having derived the most important 

^System sizes refer here to the energy scales rather than to spacial scales. 
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weak-coupling quantum kinetic equations known at present (in chapter 
we will mainly focus on three important aspects of the emergence of statis- 
tical behaviors in nanosystems. 

The first aspect (in chapter E}, is the derivation of a kinetic equation for 
the relaxation dynamics of a quantum subsystem weakly interacting with a 
quantum environment of finite size whose energy is modified by the subsystem- 
environment energy transfer. In this way, we generalize the Redfield equa- 
tion, which under its actual form only applies to quantum subsystems in- 
teracting with an infinite environment (heat bath) which is not affected by 
the subsystem. This is an important issue in nanosystems, where the non- 
relevant part of the degrees of freedom of the total system, playing the role 
of an environment, cannot always be considered infinite. The application of 
this kinetic equation to a two-level subsystem interacting in a non-diagonal 
way with a general environment, allows a better understanding of the phys- 
ical processes behind the general formalism and the explicit derivation of 
the subsystem relaxation rates as well as of the subsystem equilibrium dis- 
tribution. The second aspect of the emergence of statistical behaviors in 
nanosystems is studied (in chapter [SJ by specifying the environment op- 
erators of the precedent model in terms of random matrices. Doing this, 
one completely specifies the total system and allows a detailed comparison 
between the predictions of the kinetic equation and the exact dynamics of 
the total system. Such an analysis permits the important discussion of the 
applicability of the kinetic equation to finite quantum systems. This specific 
model, we call spin-GORM model, allows us to derive an analytical solution 
for the subsystem relaxation in the strong-coupling regime. 
Finally, the third aspect of the emergence of statistical behavior in nanosys- 
tems investigated in this thesis (in chapter El) concerns the study of quantum 
transport on a quantum unidimensional chain due to the interaction of the 
chain with its environment. We show in this model that the coherent bal- 
listic transport on the isolated chain disappears, due to the interaction with 
the environment, and is replaced by a diffusive transport. This study is 
an important step towards the understanding of the emergence of transport 
phenomena in quantum nanosystems. 



Chapter 2 



Basic concepts in 
nonequilibrium quantum 
statistical mechanics 



We establish in this chapter the main ideas necessary for the derivation of 
quantum kinetic equations. We start in section ITTI bv defining the concept 
of statistical ensemble in quantum mechanics in terms of density matrices 
and by deriving the von Neumann equation which rules the dynamics of 
these ensembles. On this ground, we discuss the main properties of the time 
evolution of a quantum statistical ensemble. In section 12.21 we formally 
derive a dynamical equation for the relevant part of a statistical ensemble 
using projection operators. We show that the main concepts that we need 
for the derivation of quantum kinetic equations in the following chapters, 
are already present at this formal level. In section 12. 3| we discuss the prob- 
lem of the positivity preservation of quantum kinetic equations and present 
the most general form of a quantum master equation preserving positivity 
(Lindblad equation). 



2.1 Quantum dynamics of isolated systems 

We introduce in this section the basic concepts necessary for the statistical 
description of an ensemble of isolated quantum systems ruled by a Hamil- 
tonian H. 
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2. Basic concepts in nonequilibrium quantum statistical mechanics 



2.1.1 Single quantum system dynamics 

In quantum mechanics, the state of a single system is described by the 
time-dependent state \tp{t)) evolving in the state space S according to the 
Schrodinger equation 

where the Hamiltonian H is a Hermitian operator in 6. The formal solution 
of the Schrodinger equation is 

-im 



|V'(t)) = e-i^*|V(0)), (2.2) 



where e R^* is the evolution operator. The eigenvalue equation of the 
Hamiltonian reads 

H\n)=En\n), (2.3) 

where {En} are the real eigenvalues of the Hamiltonian and {|n)} its eigen- 
vectors. The eigenvectors constitute an orthogonal basis in the state space 
S: {n\n') = 6nn'- The closure equation reads J2n — where I is the 

identity operator. The spectral decomposition of the evolution operator in 
£ is given by 

|V,(t))=5^e-*^"*c„(0)|rz). (2.4) 



2.1.2 Quantum statistical ensemble dynamics 

The fundamental concept of statistical ensemble has been introduced by 
Gibbs in classical statistical mechanics. A statistical ensemble is a collection 
of noninteracting identical copies of a system, the copies being in different 
states at a given time. A density matrix describes the ensemble by assigning 
a probability P^^ to the different states found in the ensemble (the probability 
corresponds to the fraction of copies in the ensemble which have the same 
state and therefore = 1). The density matrix is the central quantity 

in quantum statistical mechanics and reads 

m = T.nr{t)){r{t)\. (2.5) 

By construction, the density matrix is a Hermitian operator: = p. The 
normalization implies Trp = 1. Furthermore, a density matrix describing a 
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pure state {P^ = 6i^) has the property: p'^ = p, which is not the case for a 
density matrix describing a statistical ensemble: ^ P- In a given basis, 
the diagonal elements of the density matrix, which are always positive real 
numbers, are called the populations and are interpreted as the probability 
of being in a given state. The non-diagonal elements of the density matrix 
are called the coherences and have no classical interpretation. A detailed 
presentation of quantum statistical ensembles can be found in Ref. j25j . 
The density matrix and therefore the quantum statistical ensemble evolves 
according to the von Neumann equation 

^ = -^[H,m] = hit), (2.6) 

which is a direct consequence of the definition of the density matrix and of 
the Schrodinger equation. The von Neumann equation is the fundamental 
equation of nonequilibrium statistical mechanics. All the derivations of ki- 
netic equations start from the von Neumann equation which contains the 
full microscopic information of the statistical ensemble. 
A superoperator is a mathematical object which maps an operator into an- 
other operator. The superoperator C which generates the evolution in the 
von Neumann equation is called the Liouvillian. The formal solution of the 
von Neumann equation is given by 

p(t)=e-i^*p(0)ei^* = e^*p(0), (2.7) 

where e^* is the evolution superoperator. The Liouvillian eigenvalue equa- 
tion reads 

tp" = s,p^, (2.8) 

where {s^} are the Liouvillian eigenvalues and {p^} the Liouvillian eigen- 
vectors. The spectral decomposition of the evolution superoperator is given 
by 

/3(t) = j;a,(0)e^^y. (2.9) 
Using the density matrix 1)2. 5|) with the spectral decomposition 1)2. 4() . we get 

P^) = EE^^^-(0)<'(0)«"*^''""''"'^*l^)("'l- (2.10) 

n,n' i^ 
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By identifying (|2.9() with (|2.1U|) we conclude that the coefficients of the 
spectral decomposition (|2.9|) are given by 

a.(0) = a„„,(0) = ^P^<(0)<,(0), (2.11) 



that the Liouvillian eigenvectors are given by 

p^ = p"'^' = |n)(n'|, (2.12) 
and that the Liouvillian eigenvalues are pure complex numbers given by 

Su = -iuJu, (2-13) 
where uJi, are the Bohr frequencies of the system given by 

En -En' rniA\ 
= ^nn' = 7 • (2.14) 

n 

Notice that if the Hamiltonian operator of the system, and therefore also 
the evolution operator, is represented by a x matrix (two-index object) 
in £, the Liouvillian superoperator and the evolution superoperator are rep- 
resented hy aNxNxNxN tetradic (four-index object) in £. 

Let consider A an observable (Hermitian operator) in £. Its evolution is 
given by 

(i(t)) =TYy5(t)i = ^e-i(^"-^"')*(n|/5(0)|n')(n'|i|n). (2.15) 

n,n' 

Such an observable is an almost-periodic function of time. If the energy 
eigenvalues En are rationally dependent, it is even a periodic function of 
time. The almost-periodicity implies a recurrence theorem (which is the 
quantum analogue of the Poincare recurrence theorem) . The theorem states 
that there always exists a recurrence time around which the observable 
reaches a value arbitrarily close to its initial condition This recurrence 
theorem, which reflects the reversibility of the full dynamics, limits the ap- 
plicability of kinetic equations to quantum finite systems. In fact, kinetic 
behaviors will only be observable on time scales shorter than the Heisenberg 
time, which is defined as tn = huavie), where nav{e) is the mean density 
of states at total energy e. However, the Heisenberg time grows very fast 
with the size of the quantum system. The greater the quantum system, the 
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denser its density of states and the longer its Heisenberg time. This discus- 
sion shows that the kinetic behavior in a finite quantum system is intimately 
related to its density of states. Mathematically, of course, a true irreversible 
quantum dynamics only exist in the thermodynamical limit where the size 
of the system goes to infinity and therefore its spectrum becomes contin- 
uous. The quasi-continuous limit is a mathematical tool which sends the 
Heisenberg time to infinity. This limit can be performed by assuming that 
the eigenvalues of the spectrum are so densely distributed that the energy 
can be considered continuous En = e. Furthermore, one has to assume 
that the Hamiltonian eigenbasis representation of the relevant operators is 
a continuous function of the energy (n|A|n') = A{€,e') with e = En and 
e' = En'- Under such circumstances, the sum over the quantum states in 
the relevant quantities can be replaced by integrals over the continuous en- 
ergy J de n(e) where n(e) is the density of states of the system at 
energy e. The quasi-continuous limit of (|2.15|) would give 



{A{t)) =Ttp{t)A= de de'n(e)n(e')e-^(^-^')Vo(e,e')^(e',e). (2.16) 



Now, the calculation can exploit the power of infinitesimal calculus, the 
Heisenberg time is infinite, and a pure irreversible dynamics can occur. 



rium value of the observable, is an important quantity in nonequilibrium 
quantum statistical mechanics. It is defined by 




The time average of an evolving observable A{t), which defines the equilib- 





where we used the property 




(2.18) 



By rewriting H2.17() in the following way 



TrpooA, 



(2.19) 



we define the equilibrium density matrix 




(2.20) 



We notice that the time average of the observable is equivalent to its en- 
semble average over the invariant equilibrium density matrix p^o- Therefore, 
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according to O. Penrose j7Uj . one could assert that all finite quantum systems 
with a non-degenerate spectrum are ergodic. Unfortunately, this concept of 
ergodicity is so trivially satisfied, that it is far less important for quantum 
statistical mechanics than for classical statistical mechanics. Notice also 
that the equilibrium density matrix depends on the initial conditions. 
With the definition of the equilibrium density matrix and with the recur- 
rence theorem we have the keys to understand the dynamics of statistical 
ensembles of finite quantum systems. In a small quantum system containing 
only a few levels, the dynamics of an observable will strongly oscillate and 
the equilibrium value of the observable has no relevance to characterize the 
system. For such systems a statistical description in terms of kinetic equa- 
tions has no relevance. But when the spectrum of the considered system 
gets denser and the number of levels implied in the dynamics increases, the 
observable dynamics can display a relaxation behavior which reaches the 
equilibrium distribution on time scales shorter than the Heisenberg time. 
On this time scale, a kinetic description in terms of reduced statistical en- 
sembles is relevant. On larger time scales however, recurrences appear so 
that the kinetic description fails. The observable may sometimes leave its 
equilibrium value, even if the observable spends most of its time close to its 
equilibrium value. 

2.1.3 Liouville space 

We will end this section by a brief description of quantum statistical me- 
chanics in the Liouville space. Because it is rarely used, the notation that 
we now introduce will not be used further in this thesis calculation. How- 
ever, we want to introduce this formalism because we think it gives a nice 
intuition of the way in which superoperators act. 

The isomorphism between the matrix algebra for operators in Hilbert space 
and the tetradic algebra for superoperators in Hilbert space allows us to 
represent the operators as vectors and the superoperators as operators in a 
new space called the Liouville space. A quantum operator is represented by 
a matrix in the Hilbert space 





n,n' 



where A^n' = (n|A|n'). We can represent this operator by a ket \A ^> in 
the Liouville space. The adjoint of this operator is represented by the bra 
<C ^|. The scalar product between two operators in the Liouville space is 
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defined by 

< B\A >= TrB^A. (2.22) 



The norm in the Liouville space is given by||4|| = Y<c4|^^. A pro- 
jector and its adjoint, which arc superoperators in the Hilbert space, are 
respectively represented by the ket 

|n,n' »= |n)(n'| (2.23) 

and the bra 

<^n,n'\ = \n'){n\ (2.24) 

in the Liouvihe space. The ensemble of kets {|n, n' 3>} constitute an or- 
thonormal basis in the Liouville space 

< n,n'|m,m' >= 5„„i5„/to'- (2.25) 
The identity in the Liouville space is given by 

I = ^ |n,n' »< n,n'|. (2.26) 

n,n' 

An operator A is therefore represented by 

i = ^ < n, n'\A > |n, n' >= ^ Ann'\n, n' >, (2.27) 

n,n' n,n' 

and a superoperator S by 

S = <C n, n'\S\m, m! |n, n' m, m'\ 

n,n' m,m' 

= ^nn',mm'\n,n' »< m,m'\, (2.28) 

n,n' m,m' 

in the Liouville space. The adjoint of a superoperator in the Liouville space 
is defined as 

< n, n'\S^\m, m' >=< m, m'\S\n, n! >* . (2.29) 
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A Hermitian superoperator is therefore defined by = S. 
The Liouvillian can be represented in the Liouville space by 

C = <C n, n'|^|m, m' 3> |n, n' 3><C m, m'l, (2.30) 

n,n' m,m' 

where 

<^n,n \jC,\m,m, »= Cnn',mm' = —^{HnmSn'm' - Hn'm'^nm)- (2-31) 

We see that the Liouville superoperator is anti-Hermitian = —C. 

If the n and m indices correspond to the eigenstates of the Hamiltonian, we 

have 

(2.32) 

The von Neumann equation in the Liouville space can be written 

!^ = i|*)», (2.33) 
where the formal solution takes the form 

|/)(t)»=e^*|/)(0)». (2.34) 
The spectral decomposition of the evolution superoperator now reads 

\pii) >= XI H n,n'|e^*|m,m' »< m,m'\p{Q) » |n,n' » . (2.35) 

n,n' m,m' 

If the n and m indices again correspond to the eigenstates of the Hamilto- 
nian, we have 

\p{t) »= ^ e"""'* < n, n'|p(0) » |n, n' » . (2.36) 

n,n' 

2.2 Projected dynamics 

We will introduce in this section the Nakajima-Zwanzig projection operator 
formalism and derive from it the generalized quantum master equation. This 
method is a relatively general method to construct the relevant statistical 
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operator to describe a kinetic process by starting from the full density matrix 
describing a system. We notice that this formalism is restricted to linear 
transformations for the construction of the relevant statistical operator. In 
fact, non-linear transformations could also be considered, but we shall not 
use such ones in this thesis. By applying the projection technique to the 
von Neumann equation, one is able to obtain a general master equation 
which is exact and is the starting point of many quantum derivations of 
kinetic equations [23 EZl • We derive this generalized master equation 
because, even if it is a formal derivation, it shows how non-Markovian effects 
arise as a natural consequence of the dynamical description. We would also 
like to mention that the Nakajima-Zwanzig projection operator formalism 
shares many common points with the method of nonequilibrium statistical 
operators by Zubarev and coworkers In a sense, the Nakajima-Zwanzig 
projection operator formalism can be seen as a specific case of the non- 
equilibrium statistical operator method l9ll. However, they are very similar 
for our purpose and the Nakajima-Zwanzig formalism has the advantage of 
being simpler and more compact. 



2.2.1 Generalized master equation 

We start by defining the projection superoperators. V is the projection 
superoperator which selects the relevant part of the dynamics when acting 

on the full density matrix and Q the projection superoperator which selects 
the non-relevant part of the dynamics when acting on the full density matrix. 
The basic properties of these projection superoperators are 

P + Q = I (2.37) 

= P (2.38) 

= h (2.39) 

pQ=Qp = 0. (2.40) 

We want now to obtain an equation for the relevant part of the dynamics. 
To do this, we apply the projection operators to the von Neumann equation 
(Uni) 



= Ptpp{t)+Pttp{t) 

= hhhm + QtPm- 



(2.41) 
(2.42) 
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The equations for the relevant part and for the non-relevant part of the 
density matrix are coupled together. Integrating the evolution equation for 
the non-relevant part H2.42() . we get 



Making the change of variable r ^ t — r' in (12.431) and inserting the result 
in the evolution equation of the relevant part 1)2. 4ip . we get the generalized 
master equation 



This equation is exact and is simply a way of rewriting the von Neumann 
equation. The second term of the right-hand side of p. 44)1 is called the 
memory term. Because of this term, the evolution of the relevant part of the 
density matrix at a given time depends on its past evolution. The dynamics 
of a quantity having an evolution which depends on its history is called a non- 
Markovian dynamics. We can therefore make the very important observation 
that the evolution of the projected part of a density matrix following a von 
Neumann equation is non-Markovian. This plays a fundamental role in 
nonequilibrium statistical mechanics. The third term in the right-hand side 
of 1)2. 44() is inhomogeneous and depends on the initial condition of the non- 
relevant part. In many cases, the initial condition of the non-relevant part 
is assumed to be zero so that this third term vanishes. 
The generalized master equation (|2.44|) has been obtained independently 
by Nakajima |j66j in 1958 and by Zwanzig [HHl EH in I960. It should be 
mentioned that the generalized master equations, obtained respectively by 
Resibois j79| and by Montroll [HSl, both in 1961, where shown in Ref. [HH] 
to be equivalent to the Nakajima-Zwanzig master equation. 

2.2.2 Markovian approximation 

The generalized master equation can take a much simpler form if one per- 
forms the so-called Markovian approximation. This approximation is justi- 
fied if the memory kernel e'^^'~QjC decays to zero on a short time scale Tc 

compared to the rate of change in time of the relevant distribution Vp{t). 
The Markovian approximation consists in the following assumption 




(2.43) 




(2.44) 




(2.45) 
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If one furthermore neglects the third term of (|2.44|) representing the mem- 
ory of the non-relevant initial distribution Qp{0), the generalized master 
equation ()2.44|) takes the simple Markovian form 

^Pm = Mm, (2.46) 

where the Markovian dynamical generator is given by 

M = VC+ I dre^^^QC. (2.47) 
Jo 

A system for which a relevant distribution can be found such that it obeys a 
Markovian master equation like 1)2. 47() , is a system which can successfully be 
described by statistical mechanics and whose relevant distribution undergoes 
an irreversible dynamics to equilibrium. 

2.2.3 Slippage of initial conditions 

An important remark concerns the Markovian approximation. As afore 
mentioned, this approximation is good if the memory kernel decays fast 
enough (on time scales of the order of the correlation time Tc) compared 
to the evolution of the relevant part. However, on time scales shorter than 
Tc, the non-Markovian evolution is very different from the Markovian one. 
Therefore the Markovian equation, even if correct on the long-time scale, 
contains a bias coming from the wrong initial dynamics. This problem can 
be solved by the so-called slippage of initial conditions. This slippage is 
the rapid evolution during the time scale Tc. We can take this slippage into 
account by modifying the initial conditions of the relevant density matrix 
used for the Markovian evolution. This modification can be expressed by 
a slippage superoperator applied to these initial conditions. This slippage 
superoperator takes into account the non-Markovian effects occurring on 
short time scales smaller than Tc and maps them on the appropriate initial 
conditions of the Markovian equation. This procedure eliminates the bias 
coming from the initial non-Markovian dynamics on the long-time scale of 
the Markovian evolution. If we write the solutions of the generalized master 
equations for the non-Markovian (NM) and for the Markovian (M) dynamics 
as follows 

pp^^{t) = Tl'^'[Pp''^'{0)] (2.48) 
Pp^'it) = T,^'[Pp''{0)], 
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the slippage superoperator is defined by imposing that, on a long time scale 

(much longer than Tc), we have that Vp^^^{t) = Vp^\t). This implies 
that the initial condition of the Markovian equation is "slipped" by the 

superoperator S satisfying 

rf^^[Pp^^'\^)]=J'f^[hpp'''^{Q)]. (2.49) 

If a reasonably simple form of the slippage superoperator can be found, it is 
useful to apply it to the real initial condition of the system before implement- 
ing the Markovian scheme. In this way, one keeps a Markovian evolution 
equation which is much easier to solve than a non-Markovian equation and 
one eliminates the bias coming from the Markovian short time scale evolu- 
tion on longer time scales. 



2.3 Positivity-preserving evolution 

We present in this section the restriction imposed on the form of a quantum 
kinetic equation if the evolving quantity (the relevant projected density ma- 
trix) has to be interpreted as a density matrix. 

The density matrix operator, defined by equation (|2.5() . is a Hermitian 
operator which satisfies positivity (the diagonal elements in any represen- 
tation of the density matrix operator are real positive numbers) and which 
is normalized to unity (the trace of the density matrix operator is equal to 
one). Quantum kinetic equations are dynamical equations for the relevant 
part of a density matrix operator. If the relevant part of a density ma- 
trix has to be interpreted as a density matrix, the kinetic equation has to 
preserve, during the dynamics, the basic properties of the density matrix: 
normalization, Hermiticity and positivity. This is for example the case if 
one considers the reduced density matrix of a subsystem interacting with 
an environment because the reduced density matrix is the projection of the 
total density matrix on the subsystem state space. It is therefore natural 
to wonder if the master equation of the reduced dynamics still preserves 
positivity, Hermiticity and the norm. This question has been discussed by 
Lindblad in 1976 |56j who found that the more general form of a quantum 
master equation preserving these three properties as well as Markovianity is 
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given by the so-called Lindblad equation 
dVp{t) 



-i[H,Vm] (2.50) 



where H is an Hermitian operator. 

The mathematical formulation of this problem relies on the concept of quan- 
tum dynamical semi-groups which goes beyond the scope of this thesis. For 
further information, a mathematically rigorous discussion on this subject 
can be found in [T71 EHl Hi EH] • 
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Chapter 3 



Kinetic equations for weakly 
perturbed systems 

After having established the main concepts necessary to construct quantum 
statistical mechanics, we will derive in this chapter some of the most im- 
portant equations of nonequilibrium quantum statistical mechanics. These 
equations are presented in this chapter before the presentation of our new 
kinetic equation in chapter 0] where we shall need the material given here. A 
common point to all the theories discussed in the present chapter is the fact 
that they start from the same perturbative expansion of the von Neumann 
equation. Therefore, in section l3T] we start by performing the perturbation 
expansion of the von Neumann equation. In section 13.21 we present the 
linear response theory describing the response of a quantum system to a 
time-dependent external perturbation. In section the Fermi golden rule 
is derived because of its important role in the interpretation of the kinetic 
equations. In section 13.41 we derive the Pauli equation describing the dy- 
namics of transitions between the unperturbed states of a quantum system 
due to a constant perturbation term in its Hamiltonian. Finally, in section 
13.51 we derive the Redfield equation ruling the dynamics of a small quantum 
system weakly interacting with an infinitely large environment. 

3.1 Perturbative expansion of the von Neumann 
equation 

In this section, we develop the perturbative expansion of the von Neumann 
equation. This expansion will be the staring point of all the kinetic equa- 
tions derived in this thesis. 
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The weak-coupling assumption, even if restrictive, is one of the most pow- 
erful tools to study nonequilibrium phenomena. However, some qualitative 
arguments indicate that such an approach is not as restrictive as one could 
first think. In fact, if one wishes to put a system out of equilibrium, this has 
to be done by imposing constraints on the system which very often act at 
the surface of this system. Therefore, for reasonably large systems, because 
the surf ace/ volume ratio is low, the interaction can often be considered as 
being weak. Another argument is that even for reasonably strong interac- 
tions it is sometimes possible to incorporate the "simple" component of the 
interaction in the non-perturbed part of the Hamiltonian and to keep a pos- 
sibly weaker "complicated" component of the interaction in the interaction 
term. Finally, if the interaction is very strong, one can always renormalize 
the full Hamiltonian in such a way that the old non-perturbed part of the 
Hamiltonian becomes the new weak interaction term and that the old per- 
turbation term becomes the the new non-perturbed part of the Hamiltonian. 
Of course, such a qualitative justification to extend the application range of 
the weak-coupling approximation cannot claim to be always valid and there 
are still many situations for which non-perturbative schemes have to be used. 
However, in this thesis almost all the results are restricted to weak-coupling. 

Let us consider a Hamiltonian which is given by 

H = Ho + \V{t). (3.1) 

Hq is generally considered to be an easily diagonalizable Hamiltonian, V{t) 
is an interaction term which does not commute with Hq and A is the inter- 
action parameter which measures the intensity of the interaction term. The 
dynamics of the total system is described by the following von Neumann 
equation 

^ = tp{t) = -^[H,m] (3.2) 
= {to + \h)m = -'j:[Ho,p{t)]-X'-[V{t),m]. (3.3) 
Let us define the interaction representation of the operators as 

pj{t) = e-^°*p(t) = ei^o*/5(0e"^^°*, 
Vj{t) = e-^°^V{t) = ei^°*V{t)e-i^°K (3.4) 
Let us also define the interaction representation of the Liouvillian as 

£j{t) = e-^»*£7e^o*. (3.5) 
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In the interaction representation, the von Neumann equation H3.2() takes the 
simple form 

^ = Xtdt)piit) = -X^mt),Pi{t)]- (3.6) 
By integrating (|3.6jl . we get 

pi{t) = /5(0) + a/ dtiti{ti)p{0) (3.7) 





[ dti r dt2tl{ti)il{t2)p{t2) 

Jo JO 



p(0)-A- / dt,[Vi{t,),m] 







The perturbative expression of the von Neumann equation in the interaction 
representation is now obtained by closing (|3.7|) to order A^. 

pi{t) = p(0) + a/ dtihi{ti)p{Q) (3.8) 
Jo 

ft i-ti 



A2 f dti r dt2tl{ti)tl{t2)p{0) + 0{X^) 

Jo Jo 



m-\ J dti[Vi{ti),p{0)] 

1 2 ft fti 



I'^^'l' c^i2[V7(ti), [Vi{t2),pm] + 0{X' 
With the change of variable T = ti and t = ti — t2, we get 



pi{t) = p(0) + A/ dTCi{T)m (3.9) 
Jo 

ft rT 



+X^ f dT [ dThi{T)hi{T -t)p{Q) + 0{X^ 
Jo Jo 



10 JO 

p(o)-x'- f dT[Vi{T)rpm 



h 

^ 1^ dT £ driMn [Vi{T - r), p(0)]] + 0{X') 
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Going back to the original representation, we get the perturbative expansion 
at order two for the von Neumann equation of a system with the Hamiltonian 



A further formal step can be achieved by applying the relevant superoperator 
projector defined in section l7!^ to equation ()3.1U|) . Using the property ()2.37|1 . 
we get 



This equation will be the starting point of all our kinetic equation deriva- 
tions. The strategy is to obtain a closed equation for the relevant projected 
density matrix from equation (|3.11|) . This requires, depending on the type 
of system one studies, a careful choice of the projector. Notice that in all 
our applications, the non-relevant projected density matrix is chosen in such 

a way that Qp{0) = 0. This is a standard assumption simplifying (|3.11j) . 

3.2 Linear response theory 

Weakly perturbing a system at equilibrium in order to study its response 
to the perturbation and its return to equilibrium, is one of the main ap- 
proaches of nonequilibrium statistical mechanics. This theory, often called 
the linear response theory, has essentially been developed by Kubo in 1957 
|49j . General reviews of this theory can be found in and [01] ■ This 
approach to study nonequilibrium processes is often used to derive general 
fluctuation-dissipation relations, to compute transport coefficients and to 
prove the Onsager reciprocity relation of the phenomenological coefficients 



(ED: 



m 




(3.10) 
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relating the thermodynamical forces to the fluxes in irreversible thermo- 
dynamics. We present in this section the general ideas of this theory and 
derive the fluctuation-dissipation relation that relates the response of an 
externally perturbed system to the correlation functions characterizing its 
equilibrium fluctuations. Frequent reference will be made to appendix 
where we discuss the property of the correlation functions and their Fourier 
transforms. These results are gathered in appendix Elbecause they will play 
an important role in our later presentation of the quantum kinetic equations 
and should be easily accessible. 

Let us consider a non-perturbed system described by a Hamiltonian Hq 
at equilibrium in the state p^'^. Suppose that at time i = 0, a time-dependent 
perturbation is applied to the system. This time-dependent perturbation 
takes the system out of equilibrium and is assumed to have the form of a 
time-independent operator A multiplied by a time-dependent function f{t). 
Such a perturbation is called a mechanical perturbation. A coupling param- 
eter A measures the strength of the perturbation. The total Hamiltonian of 
the system therefore reads 

Hit) = Ho- XAfit), (3.12) 

where f{t) = for t < 0. This Hamiltonian has the form (|3.1|) with V{t) = 
—Af(t). Using the fact that the initial condition is an equilibrium density 
matrix p{0) = p^'^, the expansion (|3.11|) reads 

Pit) = p-'i + xyyTfiT)[AiT-t),p''^] + ^ J\t £ dr (3.13) 
/(T)/(r - r)e-i-^°*[i(r), [i(r - r), p"'?]]ei^o* + 0(A=^) 

where 

Ait) = e^^o*ie-:^^o*. (3.14) 
Keeping in 1)3. 13(1 only the first-order term in A, we get 

Apit) = Pit) - r = j'^ dTfiT)[AiT - t),p''] + OiX'). (3.15) 

Suppose now that we want to study the effect of the perturbation on a 
system operator 13. Its deviation from equilibrium is given by 



(ABit)) = (Bit)) - {B)eg = TiApit)B. 



(3.16) 
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Using (|3.15|) . equation 1)3.16(1 becomes 

/t fOO 

dTf{T)cPBA{t -T) = X dTf{T)xBA{t - T), (3.17) 
-oo J — oo 

where we have defined the response function (t)BA{t) 

ct>BA(t) = '-TT[A{-t),p-'^\B = '-Trp^'^im^Al (3.18) 
and the quantity 

XBAit) = e{t)(j)BAit), (3.19) 
where @{t) is the Heaviside function. 

Notice that, for a periodic forcing f(t) = /qc"*'^*, equation (|3.17|1 becomes 

{AB{t)) = 27rXfoe-''^'xBA{u;), (3.20) 
where we have defined the susceptibihty 



Xba{^) = / — e^"*XBA(t) = / —e'^^'Mt). (3.21) 

J — OO J 

We define now the correlation function of operators 13 and A which char- 
acterizes the equilibrium fluctuations of one operator with respect to the 
other: 

aBA{t) = TTp^'iB{t)A = CBA{t) + iDBAit). (3.22) 

The important properties of these correlation functions are detailed in the 
appendix It is remarkable that the susceptibility ()3.19|) and the response 
function (|3.18() can now be related to the imaginary part of the correlation 
function [see equation ()A.5|) ] 

XBAit) = eit)<PBAit) = -^eit)DBAit). (3.23) 

The Fourier transform of the susceptibility is related to the Fourier transform 
of the response function by 

XBA[uj) = — / du; V— -. (3.24) 

2 ZIT ./_^ U!' - UJ 
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Because of the relation ()3.23|) . (pBAi^) has the same properties as Dba{^) 
in equation ()A.8|) . Using these properties, one can show that 

4>BA{i^) = XBAi^^) - XAsi^)- (3-25) 

By defining the symmetric and antisymmetric parts of the response function 
and of the susceptibihty as 

4>%a{^) = {4>ba{^) + 4>ab{uj))/2 (3.26) 

4>BAi^) = i4>BA{uj) - 4>ab{^))/2 

Xba{^) = {xba{uj) +Xab{uj))/'2 
Xba(^) = (Xi?A(w) - XAi?(w))/2, 

and by using equation (|3.25|) . one finds that 

4>%A^) = 2ilm[rBA{^)] (3.27) 
0^^(u;) = 27^e[x^^(u;)], 



and that 



where 



Xba{uj) = ne[xBA{i^)] + iTm[xBA{uj)], (3.28) 



TJefo.Ml = &^'r,^vtm!^ (3.29) 

2 27r y_oo u;' -uj 



2 27r 



UJ' — UJ 



If the equihbrium density matrix is canonical {p^'^ = e /Z, where Z = 
Tre~^^°), the KMS (Kubo-Martin-Schwinger) property of the correlation 
function [see equations (|A.29() - HA.31|) ] implies that 

Cba{P; uj) = 2i^^DBA{P; u^) = -i^^MP; (3.30) 

where 

E,H = ^coth^. (3.31) 
Using equations ()3.27|) . ()3.3U|1 and HA.9|) . we also find 

ChAiP;^) = 2^^Im[rAB{P,^)] (3.32) 

UJ 

C%AiP;uj) = -2i^^ne[rABiP,^)]- 
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Equation ()3.3U|) and l|3.32() are manifestations of the so-called fluctuation 
dissipation theorem. In fact, these equations connect the response or the 
susceptibility of the system to the fluctuation properties of the equilibrium 
state. This connection is very important and occurs in many different ways 
in nonequilibrium statistical mechanics. 



3.3 Fermi golden rule 

In this section, we derive the Fermi golden rule which will give us an impor- 
tant intuitive tool for the further studies. The derivation is made starting 
from the second-order perturbative expansion of the von Neumann equation, 
in order to insist on the similarity between the different weak-coupling the- 
ories. We also show that, for a system submitted to a periodic forcing, the 
dissipated power can be related to the imaginary part of the susceptibility 
using the Fermi golden rule. 

The systems we consider here are the same ones as in the linear response 
theory. We have a non-perturbed system with a Hamiltonian Hq at equi- 
librium in the state p'^'^, to which a time-dependent perturbation is applied 
at time t = 0. This time-dependent perturbation has again the form of a 
time- independent operator A multiplied by a time-dependent function f{t) 
and the coupling parameter A measures the strength of the perturbation. 
The total Hamiltonian therefore reads 

H{t) =Ho- \Af{t), (3.33) 

where f{t) = for t < 0. The von Neumann equation ruling the dynamics 
of the system is the same as in linear response [see equation (|3.13j) ]. How- 
ever, when the operator B, that we defined in linear response theory to 
study the response of a system to an external perturbation [see equation 
(|3.1(i)l ]. commutes with Hq, the first-order expansion of the von Neumann 
equation (|3.13|) and therefore also the response function (|3.18() of linear re- 
sponse theory are zero. It is then necessary to go to the second order of the 
perturbative expansion of the von Neumann equation 1)3.13(1 . The response 
of B will in this case have the form 

{AB{t)) = TTApit)B = Y,{n\Am\n){n\B\n) = Y,^Pnn{t)Bnn, (3.34) 

n n 

where 

Apnnit) = ^ l^dT £ drfiT)fiT-T){n\[AiT),[AiT-r),n]\n). 

(3.35) 
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We choose a stationary state of Hq as equilibrium state 

p^i = \n){n\, (3.36) 

and focus on a time-periodic forcing of the form 



f{t)=cosLvt = ^ . (3.37) 



3.3.1 The derivation of the transition rates 

Using the fact that 



f{T)f(T - r) = i{(e''^2T ^-^^-i^r (3 gg^ 



where CC stands for complex conjugate, and that 

(n|[i(T),[A(r-r),p-9]]|n> = +{Y^ e^^^^^^" ~^-)-\A„^„\^ 5^„, + CC} 

n" 

-{e--|{£;„'-^;")-|A„„,|' + CC}, (3.39) 

equation I3.35i becomes 

)2 ft rT 

Ap„„{t) = Tie— dT{e"^^'^ + l) dr{ 
Jo Jo 

+ Y: (e-i ('*"+^"" -^")- + e-»("+^"-^." )-)| A„„„ p'Jnn' 

n" 

_(g-i(^+B„,-i5„)r ^ g-.l(/i^+B,.-E„,)^)|A„„,|2 }. (3.40) 



Now we make the long-time assumption which consists in taking the upper bound T of the integral 
over T in l3.4Ui to infinity, so that one can use the property 



Jo 



(3.41) 



7r(<5(a; - aj„„0 + 5{uj + uj„„,)) - iP( 1 ) 
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where tkj^^i = En — Therefore, we can evaluate the time-dependent quantities entering in 

I3.4UI as follows 

pt n oo 

Ue dr(e'"2T + l) / rf^(e-i + -B„)r_^g--i(fic^ + E„-B„,)r) (3 42) 

Jq Jo 

= +'!t(5(uJ - LUnn') + S{UJ + LUnn')) [ dT {1 + COS 2ujT) 

Jo 

1 1 /■* 

-iVi \ ) / dTsm2LuT 

= +Tr{5{uj - UJn„,) + 5{u) + UJnri>)){t + ) 

1 1 1 — cos 2ujt 

-iV{ + ) 

+7r((5(w <^„„/) + 5(<^ + w„„/))(t + ^<5((^)) 

+iV{ ^ + — ^ )V(^) 

+tt{5{w - uj„„,) + 5{u; + uj„n,))t. (3.43) 

The long-time behavior of equation (|3.4U|) is given by 

^Pnn{t) « +^{6{uj - UJnn') + S{UJ +UJnn'))\Ann'\'^t (3.44) 

-l^nn'-^ '^i^i^ - ^nn") + ^i^^ + ^^nri"))l^nn" P^- 
n" 

This expression has to be used carefully. In fact, it is valid on sufficiently 
long time scales to perform the preceding calculation, but it is not valid 
anymore on too long time scales because then the probability goes to infinity. 
It provides however a simple intuitive picture. Apnnit) is of course the 
probability to find the system in the state \n) at time t, having started to 
perturb the system at time zero in state \n') with a cosine periodic forcing. It 
is therefore quite natural to define the transition probability rate Wn'^ni'^) 
as follows 

Apnnit) = Wn'^n{i^)t. (3.45) 

If n 7^ n', the transition rate is given by the Fermi golden rule 

Wn'^ni^) = +^{^{^ -^nn') +^{^ + i^nn'))\Ann'f- (3.46) 

This rate gives the probability per unit time to jump from state \n') to |n) 
under the cosine time-dependent perturbation at frequency uj. We identify 
two possible kinds of transitions, the first coming from the absorption of 
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an energy huj by the system jumping from the initial state at energy En' 
to a state at energy En> + fiuj, and the second coming from the stimulated 
emission of an energy —huj by the system jumping now from the initial state 
at energy En' to a state at energy En' — fiuj- 

We can also define the total transition rate, which is the total probability 
per unit time to leave the initial state |n'), whatever the arrival state is, 
under the cosine time-dependent perturbation. By summing the transition 
rate from state \n') to \n) over all arrival states |n), the total transition rate 
is therefore given by 

Wn'{uj)=Y,Wn'^n{i^). (3.47) 

n 

If the spectrum is dense enough to consider that the energies vary in an 
almost continuous way in the spectrum En = e and that |A„„/p = \A(e, e')p 
is a smooth function of energy (the initial energy is denoted by e'), we can 
write that the total transition rate given by the Fermi golden rule becomes 
Wn'{u}) — > W{e',uj), where 

W{e',uj) = ^ J den{e)[S{nuj-e + e')+6{huj + e-e')]\A{e,e')\'^, 

(3.48) 

n(e) being the density of states at energy e. This expression can be rewritten 
as 



W{e',u;) = —{d^A{e',Lo)+aAAie',-u;)}, (3.49) 



where 

aAA{e,Lo) = hn{e + nu;)\A{e + huj,e)f (3.50) 

is the Fourier transform of the micro canonical correlation function defined 
in equation ()A.26|) . This is a very interesting connection between the equi- 
librium fluctuations and the response of the system to an external perturba- 
tion. Equation (|3.5nj) again shows the absorption and the induced emission 
contributions to the total transition rate. 

We now have a very nice interpretation of the Fourier transform of a micro- 
canonical correlation function: 
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• Q!Ayi(e,w) is proportional by a factor ti?/{'K}?) to the transition rate 
for a system at energy e to absorb an energy hw under a periodic 
perturbation of frequency to. 

• a^yi(e, —u}) is proportional by a factor fi? /{tt}?) to the transition rate 
for a system at energy e to emit an energy hjj under a periodic per- 
turbation of frequency u. 

If the system is in a thermal state described by a canonical distribution 
at temperature (3 before being submitted to the periodic forcing, the total 
transition rate is the thermal average of the total rate ()3.47() and reads 

W{l3,u:) = ^^Wn'{u:) (3.51) 

n' 

for a discrete spectrum and 

W{(3,uj) = J de ^ W{e,u;) (3.52) 

for a dense spectrum. 

Using the fact that the canonical correlation function is the thermal average 
of the microcanonical correlation function [see equation HA.28|) in appendix 
the total transition rate can be written as 



W{P,Lo) = —^{aAA{f3,io)+aAA{(3,-uj)}. (3.53) 



Here we have the interpretation of the Fourier transform of the canonical 
correlation function. 

• aAAiP,^) is proportional by a factor h^/{'K\^) to the transition rate 
for a system at temperature [3 to absorb an energy huj under a periodic 
perturbation of frequency oj. 

• a.AA{(3, —1^) is proportional by a factor /i^/ (vrA^) to the transition rate 
for a system at temperature [3 to emit an energy hto under a periodic 
perturbation of frequency oj. 
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System 













Figure 3.1: The system submitted to a periodic forcing of frequency u can 
absorb or emit an energy ^ at a rate given by the Fermi golden rule p.49|l . 



3.3.2 Dissipated power 

As in linear response theory, a connection can be made between the equilib- 
rium fluctuations of a system, its response to an external perturbation and 
the dissipation to which it is submitted. This shows again the many ways 
in which the fluctuation-dissipation theorem occurs. 

The power dissipated by the system submitted to a periodic external forcing 
is defined as 

P.taMlM (3.54) 

t—>oo t 

where 

{AHo{t)) = TTApit)Ho = ^Apnnit)En (3.55) 

n 

= ^ Apnn{t)En + A/J„/„/ (*)£"„' . 

This means that using (|3.44|) . we can write 

'l^i^'^^ - ^nn') + + ^nn'))\Ann'?{En - En'). (3.56) 



44 



3. Kinetic equations for weakly perturbed systems 



Assuming as afore mentioned that the spectrum is dense enough, we finally 
have that 

de^{5{nw -e + e') + 5{nw + e - e'))|^(e, e')\\e - e') 

= -7rA^a;Tm[xAA(e',w)], (3.57) 

where we used equations p.'27p . H3.23|) . (|A.8|) and the property (|A.9|) to ob- 
tain the last line. 

We see that the power dissipated by a system submitted to a periodic ex- 
ternal forcing is directly connected with the imaginary part of the suscep- 
tibility. As we have seen in the linear response theory, the imaginary part 
of the susceptibility of the system is related to the response function [see 
equation (|3.27|) ]. which in turn is related to the real part of the correlation 
function [see equation 1)3. 23(1 ]. All these connections are manifestations of 
the fluctuation-dissipation theorem. 



3.4 Pauli equation 

We will derive in this section the first quantum kinetic equation derived 
by Pauli in 1928 [OHl. In 1955, this equation has been studied in great de- 
tail by Van Hove j91j . General results on this topic can be found in Refs. 
|5H I94| 19 7j . In this chapter we will derive the Pauli equation starting from 
the perturbative expansion H3.9() of the von Neumann equation in order to 
insist on the similarity between the Pauli equation and the other weak- 
coupling kinetic equations derived in this chapter. 

The Pauli equation applies to systems whose Hamiltonian is given by 
equation (|3.ip . where the perturbation term is time- independent: 

H = Ho + XV. (3.58) 

Let us first define the notation by writing the eigenvalue equations of the 
non-perturbed Hamiltonian Hq and the full Hamiltonian H, respectively, as 

Ho\n) = Si°V) (3.59) 
H\a) = Ea\a). (3.60) 
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3.4.1 The projector 

We define now the form of the projection operators needed 

{n\Pp{t)\n') = {n\p{t)\ri)Snn' (3.61) 
{n\hm\n) = {n\p{t)\n){l-6nn')- (3.62) 

This means that the relevant density operator and the non-relevant density 
operator can be written as 

Pp{t) = J2{n\p{t)\n)\n){n\ (3.63) 

n 

km = E {n\mn')\n){n'\. (3.64) 

The relevant projector has the effect of keeping only the populations of 
the density matrix while loosing the information regarding the coherences. 
The Pauli equation is obtained using this projection operators to close the 
projected perturbative expansion of the von Neumann equation (|3.1H) where 

the terms containing the projector Q are neglected. Notice that this last 

assumption is exact if the initial coherences are zero: Q/5(0) = 0. 

3.4.2 The derivation of the PauH equation 

Writing the projected perturbative expansion of the von Neumann equation 13.111 where S/3{0) = 
in the eigenbasis of Hq, we get 

{n\Vp{t)\n') = (n|Pe^o*T'p(0)|n'> + A / dT(n\Ve^o*CiiT)Vp{0)\n') (3.65) 

Jo 

+\^ [ dT [ dT{n\Ve'^°'^Ci[T)Ci{T - T)Vp{ld)\n') + 0{\^). 
Jo Jo 

Using GMJ. 

we have 

{n\p{t)\n) = {n\p{0)\n) +X [' dT{n\ti{T)Pp{0)\n) (3.66) 

Jo 

j dT j dT{n\ti{T)ti{T -T)Pp{0)\n) +0{X^). 
Jo Jo 

We now need to evaluate {n\Ci{T)Vp{0)\n) and {n\Ci{T)ti(T ~ T)Pp(0)\n). 

{n\Ci{T)'Pm\n) = -Un\lV(T),Ppm\n) (3.67) 
n 

= -^^{n\V{T)\n'){n'\Pm\n) 

n' 

+ lY.'^n\Pm\n'){n'\V(T)\n) 

n' 

= 
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{n\ti{T)ti{T~T)Pm\n) = -^{n\[V{T),[V{T - r),Ppm]\n) (3.68) 



n' 

+ ^Y.'^n\V(T)\n'){n'\m\n'){n'\V{T - r)\n) 

n' 

+ ^Y.'^n\V(T - r)\n'){n'\m\n'){n'\V(T)\n) 

n' 

-^Y.'^n\m\n){n\V{T - r)\n'){n'\V(T)\n) 

n' 

■Y,\{n\V\n')\\n\m\n) cos " 



n' 

+ ^ E \{n\V\n')f{n'\m\n') cos 



Equation I3.66i then becomes 



{n\p{t)\n) = {n\m\n) (3.69) 



-2^ J2 \{n\V\n')\mn\m\n) - (n'|p(0)|n'>} 



rdT-''-(^"-^"')^/"+0(A3). 
Differentiating this expression with respect to time, we get 

'^^"'^^f = -2^^|(n|V'|„'>|2{(n|p(0)|n>-(n'|p(0)|n'>} (3.70) 

n' 

sin(E„ - E„>)t/h ^ ^^^3^ 

En — E„i 

In order to close the equation, we will have to assume that {n\p(t)\n) = (n|/3(0)|n) + 0{X^). 
(n\p{0)\n) can be replaced by {n\p{t)\n) in the right-hand side of I3.70i without affecting the 
equation on lower order than 0{X^). This assumption is the main assumption of this derivation 
and is common to many derivations of weak-coupling kinetic equations. Simplifying the notation 
and defining fv^J^m' = — E„i , equation I3.7UI becomes equation I3.71i . 

The non-Markovian version of the PauU equation for a discrete spectrum 
reads 

= -^^E\^nn'f^^^{Pnn{t)-p^,^,{t)}. (3.71) 
at ^ UJnn 

n' 

On long times scales, because limf^oo = 7r5(ci;„„/), equation (|3.71|) 

can be simplified to 

'^^'^'^^ ^ = X] \Vnn'?'5{uJnn'){Pnn{t) - Pn'n'{t)]- (3.72) 

n' 
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Equation H3.72() is the Markovian version of the Pauh equation for a discrete 
spectrum. The Pauh equation is however typicahy written for a system with 
a continuous spectrum. If a denotes the eigenvalues of ah the operators 
that commute with Hq (all quantum numbers of the unperturbed system 
except energy), the continuous version of H3.71() and (|3.72|) is obtained in 
the following way 

Pnnit) ^ P{e,a,t), (3.73) 
\Vnn'\^ ^ \{e,a\V\e',a')\^ 

J de' J da'n{e ,a'). 

n' 

The correlation function avv{^, ^) defined in equation (|A.26|) becomes avvi^, ct, a', w). 
Similarly, the total transition rate given by the Fermi golden rule W{e',uj) 
defined in equation (|3.49j) becomes W{e' ,a,a' ,u;). We now have 

avv{^: a', w) = hn{e + huj,a')\{e + huj, a'\V\e,a)\'^ (3.74) 

and 

W(e' ,a,a' ,oj) = -—7^{avvU',a,o^',^) + ctvv{'^',Ci,a',—uj)}.(3.75) 

Taking the continuous limit in equation ()3.71|1 . we get the non-Markovian 
version of the Pauli equation for a continuous spectrum 

dP{e,a,t) _ /rfe' /■rf«'^(e>')|(e,a|y|e>')|2 (3.76) 



dt h 



fT^, X I I XT sin(e — e')t/^ 

{P(e,a,t) -P(e',a',t)}- ^ 



e-e' 

On long time scales, the sum in ()3.71|) is limited to states for which |e — 
e| = 0{h/t). The number of states in this domain is roughly n{e)\e — e'\. 
This number has to be larger than unit to convert the sum in an integral 
t/{n{e)h) <^ 1. We therefore conclude that the continuous limit taken to go 
from H3.71() to 1)3. 76() is valid for time scales shorter than the Heisenberg time 
tff = n{e)h. This result applies to all the other kinetic equations derived in 
this thesis. 

Equation 1)3. 76() can be rewritten in terms of the Fourier transform of the 
correlation function using the change of variable huj = e — e' as follows 

dP{e,a,t) f , f , ,~ , I ^ /o^^vx 

— "^T^ \ duj \ da ayyxe^a^a ^—Lo) ("J-'') 



dt ff 



{P(e, a, t) - P{e - huj, a',t)}- 
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Using the long-time property lim^^oo ^i^j^ = TT6{io), we get the continu- 
ous version of equation (|H.72|) . which is the Markovian version of the Pauh 
equation for a continuous spectrum. This equation is the original equation 
derived by Pauh in 1928 [HE]: 



dP{e, a, t) 
Jt 



J daW{e, a, a, 0){P(e, a, t) - P{e, a, t)}, (3.78) 



where 

W(e, a, a ,0) = iTT-^avvi^, a ,0) = 271— nie, Oi)\ie,a\V\e,a)\^ (3.79) 
is the total transition probability given by the Fermi golden rule at frequency 



zero [see equation (|3.49j) ]. 

Notice that the long-time limit can be justified when t is large enough to 
assume that W {e^a^a\0{h/t)) T^(e, a, a', 0). Let us call t^ the minimal 
value of t which satisfies this relation. 

The Van Hove limit jHJ is often cited in the literature because it is the limit 
for which the Pauli equation is exact. This limit consists in successively 
taking the following limits: 

• N ^ oo 

• A ^ and t — > oo with A^t fixed. 

The Van Hove limit is physically not very useful. When this limit is not 
taken, the Pauli equation is valid for times t such that < t < tn- 

The Pauli equation is an important equation. It is the simplest quantum 
kinetic equation and a simple interpretation in terms of stochastic dynamics 
is possible. The dynamics can be seen as resulting from random transitions 
(given by the Fermi golden rule) between eigenstates of the non-perturbed 
Hamiltonian, caused by the presence of a perturbation term in the total 
Hamiltonian. The stochasticity of the Pauli equation is a consequence of the 
reduction of the entire density matrix describing the deterministic dynamics 
to the population dynamics. The Pauli equation is an irreversible equation 
which allows to describe quantum non-equilibrium processes. This equation 
has however a limited range of applicability. In fact, in order to be able to 
describe coherent processes, the coherences of the system have to be included 
in the description. This is what the Redfield equation achieves as we will 
see in the next subsection. 



3.5 RedHeld equation 
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a 



Figure 3.2: The Pauli equation describes transitions occurring between levels 
of different quantum numbers a belonging to the same narrow energy shell. 

3.5 Redfield equation 

As already mentioned in the introduction, one of the common ways to study 
the relaxation to equilibrium of a system is to let it interact with its environ- 
ment. We will derive in this section the Redfield equation which is the best 
known quantum kinetic equation describing the dynamics of a small sub- 
system interacting with an infinitely large environment. In this scheme, the 
irreversible subsystem dynamics arises as a consequence of the fact that the 
description only focuses on the reduced dynamics of the subsystem which is 
a projection of the full dynamics on the degrees of freedom of the subsystem. 
The Redfield equation has been first derived semi-empirically by Redfield in 
1957 in the context of NMR |77j . Many systematic derivations have been 
performed since then [111120113111131^113 although we find equations that 
not always bear the name of Redfield in these different references. This 
equation has often been used in various fields ranging from NMR jHH |H3 to 
quantum chemistry [HHl , electron transfer theories , and condensed-phase 
reaction dynamics |52j . 

The systems we consider in this section are composed of two parts, a 
small quantum subsystem with a Hamiltonian Hs and an environment with 
a Hamiltonian Hb- These two parts are interacting with each other through 
an interaction term. The magnitude of this term is measured by the coupling 
parameter A. The interaction term can be written as a sum of products of 
Hermitian operators of the subsystem with Hermitian operators of the 
environment iJ^. The total Hamiltonian of such systems reads 




(3.80) 
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This Hamiltonian is equivalent to the Hamiltonian H3.1|) , where Hq = Hs + 
Hb and V = Sk-Bk- We use the notation 

Hs\s) = Es\s) (3.81) 
HB\b) = Eb\b). (3.82) 

The dynamics of the total system is given by 

^ = tm = {h + tB + xtj)p{t) (3.83) 

In the interaction representation the coupling operators become 

B^{t) = e-^B^B^ = ei^^^B^e-i^B\ (3.84) 
and the interaction representation of the Liouvillian reads 

3.5.1 The projector 

The subsystem dynamics is described by the reduced density matrix of the 
subsystem 



PS 



{t) = TrBpit) = E \')(''\ E^^^I^WI^'^)- (3.^ 



Let us define the following projector over the relevant part of the density 
m at r i^c 

s,s' b 

and the projector on the non-relevant part of the density matrix as 

^P(t) = E E {sb\mWb')\sb){s'b'\. (3.88) 

s,s' b,b'^b 
Equivalently, one can write that 

{sb\Pm\s'b') = {sb\p{t)\s'b)6w, (3.89) 
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and 



{sb\Qpit)Wb') = {sb\mWb'){l - Sbb'). 



(3.90) 



The reduced density matrix of the subsystem is obtained from the relevant- 
part projector as follows 



This reduced density matrix contains all the information about the system, 
whereas all the information about the environment has been traced out and 
is therefore lost. The reduced density matrix of the subsystem is the central 
quantity in the Redfield theory and the standard derivation of the Redfield 
equation consists in getting a closed equation for the reduced density matrix 
of the subsystem from the von Neumann equation. All our derivations start 
from the projected weak-coupling expansion of the von Neumann equation 

p.llf) . As in the Pauli equation derivation, in order to eliminate the Qp{0) 
terms in the equation (|3.11|) . we always choose an initial condition for which 
this quantity is zero. In the Redfield theory, this corresponds to take an 
initial condition of the total system which is the tensor product of a general 
initial condition for the subsystem with an invariant density matrix for the 
environment (i.e. a density matrix which is diagonal in the eigenbasis of Hb 

and therefore with the property C-bP^b ~ *-*)• This initial condition reads 



TYBVp{t)=TTBp{t)= Ps{t)- 



(3.91) 



p(0) = P5(0)p 



(3.92) 



Therefore, Vp{<d) = Ps{'^)p'b and Qp{<d) = 0. 



3.5.2 The derivation of the Redfield equation 

We now explicitly start the derivation. 

The projected von Neumann equation lij.lli becomes 




(3.93) 



We first need to evaluate ti(T)Vp{0) and ti(T)Ci{T - t)Vp{0). We have that 



Ci{T)vm 



J2 S.{T)ps{0)B^{T)p^^ +'-Y. PsmS.mpl'B^T) 



(3.94) 



52 3. Kinetic equations for weakly perturbed systems 

and 

tiiT)ti(T-T)Pm = -LY^[s,{T)B^{T),[S^,{T-T)B^,{T-T),Pm]] 

PS{0)S^,{T - r)S4T)pl^B^,{T - r)B4T), (3.95) 

SO that equation 13.931 becomes 

psit) = e^^*/3s(0)Trsp^« (3.96) 
A . 



-'I E J dTS4T)ps{0)TrBB4T)p%^ 



.A . 



+'jiJ2^^''J dTpsiO)S.{T)TrBp'j^B^{T) 



~ E /' '^^ J'^ dTS^{T)S^,{T - T)ps(0)TrsB4r)B«'(T - t)p 



Jo 



~ E fy^^ drS^{T)psmS^,{T - r)T^rBB^{T)p^^B,,(T - r) 

A2 



e'^st / dT / drS^.iT - t)ps{0)S4T)TvbB^, (T - t)pI^B4T) 



JO 

E [ '^^ [ drpsmS^, (T - r)54T)TVsp«'B,, (T - t)B,{T) + ©(A^). 

DifTcrcntiating equation 13.961 with respect to time, we get 
dps it) 

= I-. 

. A T — \ A y r ^+ . .fn^ .A 



CsPsit) (3.97) 



E 54e^^*Ps(0)}TrsP^''B« + i- X^{e^s*ps(0)}5«TtsP^''B« 
J2 f dTSj,,{-T){e^stp^^o)}"CrBP%'B^(t)B^,(t-T) 

E /J'i^S«{e^^'ps(0)}4'(-^)Trs/3^'B,-(t-r)B4t) 

-^E /Jdr5,,(-r){e^s*ps(0)}5.Trs/3^^B4t)S«'(t-r) 

■^E f dT{e^stp^^Q^ys^,{^T)S^TrBP%'B^,it-T)B4t) + 0{\^). 
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In order to close the equation, we use the fact that ps(t) = e'"S*pg(o) + C'(A^). This does not 
affect equation i3.97i to orders lower than 0{X'^). 

The final result is the Redfield equation which reads 



i:5P5(t)-4E['^-'5^(*)K^-) + ^E r^^{ (3.98) 

+oi*KK'{T)SKPs{t)SK'{-T) + ai^^i{T)Si^'{-T)ps{t)S^ }, 



where the expectation of the environment couphng operator is given by 

{B^) = TiBp'^B^, (3.99) 

and where the environment correlation function by 

a^^^ir) = TrBp'^B^'{T)B^ = aL,(-r). (3.100) 

Equation (|3.98|) is the non-Markovian version of the Redfield equation. This 
equation preserves trace and Hermiticity. The Markovian version of the 
Redfield equation (as originally derived by Redfield) is simply obtained by 
replacing dr by dr in equation (|3.98jl . 

The correlation function appearing in H3.98() has the same form as the one 
we found in linear response theory [section 13.21 equation 1)3. 22() ]. The gen- 
eral properties of this correlation function are detailed in appendix This 
correlation function characterizes the equilibrium fluctuations of the envi- 
ronment. 

Notice that if {Bf^(t)) 7^ 0, with the redefinition 

Hs ^ Hs + \Y,S^{B^) (3.101) 

K 

Y^S.B^ ^ ^5«S,-A^5«(4), (3.102) 
the 0(A)-term can be eliminated from equation H3.98() . 
3.5.3 Slippage of initial conditions 

An important remark is that the Markovian version of the Redfield equation 
does not in general have the Lindblad form (see section \'2.'6^ and is known 



dpsjt) 
dt 
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to violate positivity if the initial condition is near the border of the set of 
physically admissible density matrices [21 69 , 90 . We notice however 
that the trace and Hermiticity are preserved by the the Markovian Redfield 
equation. The non-Markovian version of the Redfield equation is believed 
to preserve positivity. There is no general proof for this statement, but 
we do not know any work having infirmed this statement. Therefore, the 
non-preservation of positivity in the Markovian Redfield equation is meant 
to arise as a consequence of not having taken into account the initial non- 
Markovian dynamics which corresponds to the stage of the dynamics where 
correlations are established between the subsystem and its environment from 
the uncorrelated initial condition (|H.92|) . This highly non-Markovian initial 
stage of the dynamics occurs on time scales of the order of the environment 
correlation time, which is the typical time necessary for the correlation func- 
tion of the environment to decay to zero. A slippage (displacement) of initial 
conditions has been proposed by various authors I3f)[ 190) to correct this ini- 
tial highly non-Markovian dynamics which is not taken into account by the 
Markovian equation. The general method of the slippage of initial condi- 
tions has been mentioned in section 12.2.31 In this way, it can be shown 
that when the Markovian evolution starts from an initial condition modi- 
fied by the slippage superoperator, it preserves the positivity and gives the 
same results as the non-Markovian equation on time scales longer than the 
correlation time of the environment j3()| I9()j . 

3.5.4 Delta correlated environments 

An important observation can be made when the environment has a very 
fast dynamics compared to the subsystem dynamics. In this case, the envi- 
ronment correlation function decays on a very short time scale compared to 
the subsystem dynamics and can therefore be considered proportional to a 
delta distribution in time 

a«,,(r) = 2D,,,5(r). (3.103) 

In this case, the Markovian Redfield equation takes the Lindblad form and 
therefore always preserves positivity. 

This can be understood by the following argument. D^^' is Hermitian and 
can be diagonalized as D = IJ^dU where d is a diagonal matrix with eigen- 
values dfi. By introducing the Lindblad operator defined as follows 

L = Vd'nY.U--'^-'' (3-104) 

k' 



3.5 RedHeld equation 



55 



the Markovian Redfield equation becomes of Lindblad type H2.50() 
dpsit) i 



[Hs,m] (3.105) 

2 
^2 



dt h 



3.5.5 Markovian approximation 

To simplify the notation, we shall eliminate the summation over k and con- 
sider a single coupling term. This does not affect the discussion and the 
summation is easy to restore. The Redfield equation 1)3. 98(1 becomes 

kit) = hsPs{t)-i^[S,ps{t)m + ^ l^dT{ (3.106) 

-air)SS{-T)psit) - a*iT)ps{t)Si-T)S 
+a*{T)Sps{t)S{-T) + a{T)S{-T)ps{t)S }, 

To study the population and the coherence dynamics, we will project the 
Redfield equation in the eigenbasis of the subsystem. 

Using the notation {s\ps{t)\s') = Pssi{t), equation 13.1061 becomes 

ss' 

(r)e- * E,, )r 

+Ss-sP-s-s> (t)S-s,,, a* (r)e- i^'^^' -'^^'^^ 
+5««p,,,(t)5,-,,m(r)e-i(-E=-E=^)- }. (3.107) 

Using a{t) = J duja{Lu)e~^'^* , the property a{uj) = a*{iLj), and defining tiui^ = {Es — E^i), we 
get 







-SssS,,,p,,,,{t)aiu;)e''-^^'-'^^^ 
-pss{t)S,s'S,,,,a{u;)e-'(^!is'-^^^ 
+SssPss'{t)S,>s'a{'^)e-'^'^^'^'-'^^^ 
+Ss,p-ss'{t)S,,,,a{Lo)e''-^^--^^^ }. (3.108) 
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Taking the Markovian limit and using 







^^g±i"T ^ T,s{uj), (3.109) 



we get equation l3.11Ui . 

The Markovian version of the Redfield equation in the eigenbasis of the 
subsystem reads 

Pss'{t) = -iuJss'Pss'it) - i^'^{SssPss'{t) - Pss{t)Sss'){B) (3.110) 

s 



A r 

+^ ^{ SssSss'Ps's'{t){ -7ra{ujs's) - i J dwV 



LUgig — UJ 



+Pssit)Sss'Ss's'i -T^Oi{^ss') +i j dojV 

dujV ) 

Ldgigi — LO 

/(y (cj) 
dujV ) }. 
UJss - ^ 



No further approximation have been done with respect to the Markovian 
Redfield equation. 



Population dynamics 

To study the population dynamics we will now neglect the contribution of 
the coherences to the population dynamics. We get 

Pss{t) = -2TT^^^\Sss?Ol{-UJss)Pss{t) (3.111) 

s 

+ 27r-^ \Sss?Ol{LOss)Pss{t). 
s 

This master equation describes the population dynamics of the subsystem. 
This equation is commonly used is quantum optics 16 . Its simple structure 
allows us to easily understand the physical processes it describes. Accord- 
ing to our results of section 13.31 the quantity a{—0Jss) is proportional to 
the Fermi golden rule transition rate for the emission of an energy fujOgs by 
the environment. This energy is therefore absorbed by the subsystem which 
undergoes a transition from the state s to the state s. Since this transi- 
tion depopulates the state s, this term appears as a negative contribution 
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to the Pss{t) dynamics. a{ujss) is proportional to the transition rate for the 
absorption of an energy fkVss by the environment. This energy is emitted 
by the subsystem in the state s which jumps to the state s. This process 
contributes positively to the pss{t) dynamics. 

An important remark, which will become fundamental in chapter |1J is that 
the state of the environment is not affected by the emission or by the ab- 
sorption of a quantum in the Redfield scheme. This becomes a problem if 
the environment is finite and of small size. We finally notice that equation 

Environment 




Figure 3.3: The two processes of emission (on the right) and of absorption 
(on the left) of energy from the environment to the subsystem are schemat- 
ically represented, b schematically represent the state of the environment 
which is either of given energy in the microcanonical case or of given tem- 
perature in the canonical case. 

H3. 111(1 can satisfy detailed balance. If pti denotes the stationary density 
matrix, the detailed balance condition reads 

a{-uJss)p2 = c^{lOss)pTs. (3.112) 



The environment correlation functions generally used in the Redfield equa- 
tion are canonical. In this case, if the temperature is 1//3 and using the 
KMS property ()A.3fl|) . the detailed balance condition allows us to find a 
Boltzmann distribution for the subsystem equilibrium distribution 
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If one applies the detailed balance condition to the Redfield equation with 
microcanonical environment correlation functions (at energy e), we get the 
following equilibrium distribution for the subsystem 

% = 3.114 

When the environment is large, the equivalence between the microcanonical 
and the canonical ensemble should allow us to understand how the Boltz- 
mann subsystem canonical distribution ()3.113|) can emerge from the micro- 
canonical one (|3.114j) . However, the Redfield theory provides no answer to 
this question, indicating that the Redfield theory, when applied to environ- 
ments in microcanonical states, may contain some non-physical features. In 
chapter m we will see how our new kinetic equation solves this problem. 

Coherence dynamics 

To study the coherence dynamics, we isolate one coherence and neglect its 
coupling to the other coherences and populations. We get an equation of 
the form 



where 



where 



Pss'ii) ~ {-^'^ss' -lss')Pss'{i)-. (3.115) 



Eg — Egi 

Uss' = r , (3.116) 



Es = E, + ^{B{t))Sss + ^ E l-^-l' f du;V^^ (3.117) 
n J uJss — ^ 

is the Bohr frequency modified by the so-called Lamb shifts and 

7.,.' = -Tij^iSss- Ss,s.fa{fi) (3.118) 
+71"^ ^^{\Sss\'^a{uiss) + \Ssis?a{u:s's)) 

s^s 

the damping rate. 

We thus have two kind of effects on the dynamics of the coherence of the 
subsystem induced by the environment interaction. The first is a Lamb 
shift which modifies the frequency of the free oscillation of the subsystem 
coherence and the second effect is the irreversible damping of the coherence, 
this last effect being at the origin of decoherence. 



Chapter 4 



New kinetic equation for 
finite environment 

In the present chapter, we derive one of the main results of this thesis. We 
have seen in the introduction that the advances in nanotechnologies raise 
the question of the vahdity of the standard methods of statistical mechanics 
at small sizes. We here present a new kinetic theory to describe the dynam- 
ics of a subsystem interacting with an environment. The standard theory 
to describe such systems is the Redfield theory. However, in this theory, 
the environment is implicitly assumed to be infinite and therefore not af- 
fected by the subsystem dynamics. This assumption is probably not valid 
for many nanosystems for which the environment of the subsystem is itself 
of finite size. Examples of such nanosystems are atomic or molecular clus- 
ters in which a two-level impurity is enclosed. Our new formalism permits 
to derive a more general kinetic equation than the Redfield equation, which 
takes into account the possible finiteness of the environment. This equation 
has been derived recently by Esposito and Gaspard in Ref. j22j . 

Our new equation is derived in section 14.11 The population and the 
coherence dynamics are studied in section 14.21 In section I4.ISL we show 
that our new equation transforms into the Redfield equation in the infinite 
environment limit. 

4.1 Our new kinetic equation 

As in the Redfield theory, we consider a system composed of two parts, a 
quantum subsystem whose Hamiltonian is Hs and an environment whose 
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Hamiltonian is Hb- These two parts are interacting with each other via an 
interaction term. The magnitude of the interaction is characterized by the 
couphng parameter A. The interaction term is made of a sum of products of 
Hermitian operators of the subsystem with Hermitian operators of the 
environment -B^ . The Hamiltonian of the total system reads 

H = Hs + Hb + XY1 (4-1) 

K 

This Hamiltonian has the form where Hq = Hs + Hs and V = 

S^B^^. We use the notation 

Hs\s) = Es\s) (4.2) 
HB\b) = Eb\b). (4.3) 

The dynamics of the total system is ruled by the von Neumann equation 
^ = tm = {^s + tB + xtj)p{t) (4.4) 

hi 

In the interaction representation, the operators become 

pi{f) = e-(^«+^«)V(i) = e^(^^+^^)*/5(t)e-i(^^+^^)*, 
Ut) = e-^sts^ = ^fists^ifist^ 

B,{t) = e-~^^'B^ = ei^^'B^e-i"^\ (4.5) 
and the von Neumann equation (|4.4() reads 

^ = \h{t)pi{t) = -x'-Y^[S,(t)Mt),pi{t)l (4.6) 

K 

where we have defined the interaction representation of the Liouvillian 
4.1.1 The new projector 

The central quantity P{HB,t) of our new theory is defined by the following 
projection operator 

s,s' b ^(^S) 
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where n{HB) is the density of states of the environment. The projector on 
the non-relevant part of the density matrix therefore reads 

s,s' h,b'^b 

Equivalently, we can write that 

{sb\Pp{t)\sh') = {sh\p{t)\s'h)6w ^ (4.10) 

and 

{sh\hp{t)\sh') = {sh\p{t)\s'b'){l - 5w). (4.11) 

Our new kinetic equation is an evolution equation for the quantity P{HB,t). 
The only information lost by describing the total system dynamics in term 
of the quantity P{HB,t) is the information regarding the environment co- 
herences. However, the information on the subsystem populations and co- 
herences as well as the information on the environment energy is kept. The 
reduced density matrix of the subsystem is obtained from P{HB,t) by the 
tracing over the environment states 

ps{t) = TTBVm = TvB^^^^. (4.12) 

Similarly as for the Redfield equation, the initial conditions which eliminate 

the terms containing Qp{0) in equation ()3.11|) are initial conditions which 
are diagonal in the eigenbasis of Hb- These initial conditions read 

P(0,.™. ,«3) 

n{HB) 

We therefore have Vp{0) = P{HB,0)/n{HB) and Qp{0) = 0. 



4.1.2 The derivation of our kinetic equation 

We start now the explicit derivation of the new kinetic equation. 

The projected weak-coupling expansion of the von Neumann equation 13. Hi becomes 

^ p^its+ts)t mB,0) ^^p^Cts+ts)t r dTCj{T f^^^'°^ (4.14) 
n{HB) n(HB) Jo n(HB) 

Jo Jo n(HB) 
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Projecting this equation in the eigenbasis of the environment we get 

n{Et) n{Et) Jo niHg) 

Jo Jo n(HB) 

We now need to evaluate (6|i;j(T)(P(i5's, 0)/n(_H's))|b> and {b\ti(T)ti{T-T){P{flB,0)/n{HB))\b). 
We have that 

(b\UT)^^^^\b) = -'-J2{b\lS4T)B.{T),^^^§2^]\b) (4.16) 
n(Hs) n(HB) 

V "-(^b) n{Eb) 



h 



and 



{b\tiiT)ti{T-r f^^^'°\ b) (4.17) 
n{HB) 

= -i E (fe|[S.(T)B.(T), [S^,{T - t)B^.{T - r), 

= E E 54T)5«'(T- r):^i§l^(fe|B«|fe'>(b'|B«'|fe>e^(^^-^''''- 

+iEE^«(^)^^^§^^«'(^--)<^i^«i^')<^'i^«'|6>ei<^^-^''')^ 
+iEE^'''(^--)^^^S^^''(^)(f'i^'='i^')<^'i^''|f')-"^^''^"''''''^ 

Equation I4.15i becomes 



n(S6) n{Eb) h ^ Jo n{Eb) 

n ^ Jo ri{Ei,) 



lXe^stJ2 dTS4T)^^^{b\B.\b) (4.18) 
h ^ Jo n{Eb) 



K.n' b' 



Jo 



-S^{T)S^,{T-t) ^^^'"^^ (b\BJb')(b'\B^,\b)e-k'^'^''-'^'''^^ 
niEb) 

+54T)^^^§^5.,(T-r)(6|B«|6'>(6'|B„,|;'>e*(^''-^'.')- 
+S.'(T - T)^^^f^S.iT)(b\B^,\b')(b'\B.\b)e-i^^^-''^'^^ 



n{Eb) 



-S,,(T-T)S4T)(fe|B,,|6'>(fe'|S.|fe>e"«^ } + Ci(A^)- 
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Differentiating equation I4.18i with respect to time, we get 



n(Ei,) h ^ 

K,,K,' b' 

n(Eb) 

n(_Bi,') 

n{Ey) 

-^^^^^S^,{t-r)S^{t){b\B^,\b'){h'\B^\b)e-i'^^^-^y^-} + 0{W'). 



fn order to close the equation, we use the fact that P{HB,t) = e'^st p(^Hb ,0) + O(X^). This does 
not affect equation I4.19i to lower order than 0{X'^). 

Finally, our new kinetic equation reads 



,P(E^ = ,4.20) 

at n[Eb) n[Eb) n ^ n{Eb) 

\2 ft 







n{Eb) 

n{Eb') 
n{Eb) 
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If the spectrum of the environment is dense enough to be considered quasi- 
continuous, performing the quasi-continuous hmit as follows 



n{Eb) ^ n(e) (4.21) 
PiEb,t) ^ P{e,t) 

{h\B,\b') ^ i?«(e,e'), 
our kinetic equation for quasi-continuous environments becomes 

>(e,t) = tsP{e.t)-'-\Y,[S^,P{e,tmB^\e) (4.22) 

K 



-S^S^,{-T)P{e, t)n{e')B,{e, e')B^,{e\ e)e 
+Sj{e',t)S^,{-T)n{e)B,{e,e')B,,{e',e)ef^'-^> 
+SA-r)P{e ,t)SM^)B^.{e,e')B^{e' ,e)e-^n^'-'"^^ 
-P{e,t)S^,{-T)SM^')B,,{e,e')B^{e',e)e-i^'-''^^}. 



Using the correlation function [defined by equation HA.25|1 ] 



{e,t) = Jde'n{e)B^{e,e')B^,{e',e)eTA^-^'^' (4.23) 



and its Fourier transform [defined by equation HA.26|) ] 

a^K'i^i^) = hn{e + hLo)Bf^{e,e + tiuj)Bi^'{e + huj,€), (4.24) 
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we can rewrite equation H4.22() as 

Pie,t) = tsPie,t)-^XY,[Sn,Pie,t)]{e\B^\e 

K 

k,k' 



'—l)ei(^-^'^^Sj^,{-T)Pie,t) 



(4.25) 



+aL,(6',i^)ei(-^>5.P(6',i)^,,(-r) 
+a..'{e','—^)e-i^^-^'^-S,,{-T)P{e',t)S^ 



' '^e-i(^-^>P(e,t)S«,(-T)5,}. 



With a change of variable in the energy integral, our kinetic equation for a 
quasi-continuous environment becomes 



P{e,t) 



hsP{e,t) - ^XY,[Sn,P{e,t)]{e\B^\e) + ^Yl f '^^ I i 



d,^,{e,oj)e-''^^Sj^,{-T)P{e,t) 
,,{e,oj)e'''^P{e,t)S,,{-T)S^ 



huj, oj)e"^^S^P{e - huj, t)S^> (-r) 



+a««,(e - hw,uj)e-''^^S,:{-T)P{e - hw,t)S^ }. (4.26) 



This equation preserves trace and Hermiticity. The Markovian version of 
this equation is obtained by replacing the time integral dr by dr in 
equation ()4.26|) . 

In the case of a quasi-continuous environment, using equation ()4.12p . the re- 
duced density matrix of the subsystem can be constructed from the solution 
of equation (|4.26l) as follows 



Ps{t) 



de 



( Pii(.e,t) Pu{e,t) 



PiNsie,t) \ 
-P2Afs(e,i) 



(4.27) 



\ PNsl{^,t) PNs2{e,t) ... PNsNs{^,t) / 
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where Ns is the number of levels of the subsystem. 



Notice that our new kinetic equation is a closed evolution equation for the 
quantity P{HB,t). However, no closed evolution equation can be obtained 
for the reduced density matrix of the subsystem. Stronger assumptions are 
needed to derive a closed equation for the reduced density matrix of the 
subsystem, as in the Redfield theory. 

We see that the dynamic of our new kinetic equation implies different en- 
vironment energies. To better understand the energy transfers between the 
subsystem and the environment, we will project our kinetic equation in the 
eigenbasis of the subsystem in the next section. 



4.2 Markovian approximation 

To simplify the notation, we will eliminate the summation over k in our 
kinetic equation. This does not affect the discussion and the summation is 
easy to restore. Equation 1)4. 26() becomes 

Pie,t) = tsP{e,t)-'-X[S,P{e,t)]{e\B\e) + ^ l^dr I duj{{A.28) 

-a{e,u;)e-''^^SS{-T)P{e,t) 

-a{e,uj)e''^^P{e,t)S{-T)S 

+d(e - huj,uj)e''^^SP{e - huj,t)S{-T) 

+a{e - hw,uj)e-''^^S{-T)P{e - huj,t)S }. 

Now, we will project this equation in the eigenbasis of the subsystem and 
perform the Markovian approximation. 

Projecting equation i4.28i in the eigenbasis of the subsystem, we get 



(4.29) 



3 



-a(e, c^)e-»("=-=-' -")"P.,(e, t)S,,,S,,,, 

+a{e- nw,uj)e''-'^^''-'^^''SssPss'{t- ti^,t)Ss'8' }, 



where ui^^i = [Eg — Egi)/h. Remembering the property 



JO 



^^giic^r ^ _,_ 7r5(aj), (4.30) 
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s's' 



we can now perform the Markovian approximation in 14.291 . 

The Markovian version of our kinetic equation in the eigenbasis of the sub- 
system reads 

i ' 

Pss'{e,t) = -iuJss'Pss'{e,t) - -X{e\B\e) '^{SssPss'{e,t) - Pss{e,t)Sss') + ^ { 

s ss' 

+SssSss'Ps's'{(^,t) ( -7Ta{e,uJs's) - i [ duV ] 

+Pss{e,t) ( -Tra{e,iL>ss') +i [ duV "^'^'^-^ ] Sgs'Ss's' 
V J ll>ss'-ujJ 

+Sss vri^jg/(e - hujs's',t)a[e - hiOs's'-,^s's') - « / dujP -^s' 

V J ^s's' J 

,cfvf k \, ■ Pss'{e-hw,t)a{e-nuj,uj)\ 
+bss [t^Pss'{(^- nwss,t)a{e- fujss,u}ss) + i / dujV h^s's' }• 

V J UJgs - J 

(4.31) 

4.2.1 Population dynamics 

To study the population dynamics, we now neglect the contribution of the 
coherences to the population dynamics. We get 

Pss{e,t) = -27r-^^|5,5pa(e,-a;5,)P,,(e,t) (4.32) 

s 

A^ 

+27r^ ^ |5s5pa(e - huJss,^^ss)Pss{e - hjjJss,t) 

s 

This master equation describes the population dynamics of the subsystem 
and can be physically interpreted in a similar way as the Redfield population 
equation (|3.11H) . 

The quantity d(e, —u)ss) is proportional to the Fermi golden rule transition 
rate for the emission of an energy hjj ss by the environment at an energy e. 
This later energy is therefore absorbed by the subsystem in the state s which 
jumps to the state s. Because this transition depopulates the state s, this 
term appears as a negative contribution to the dynamics of Pss(e,t). This 
term is similar to the negative contribution term of the Redfield population 
equation (|3.111|) if one applies it to microcanonical environment correlation 
functions. The quantity a.{e — hwss^^ss) is proportional to the transition 
rate for the absorption of an energy fuxiss by the environment at an energy 
e — hujss- This energy is thus emitted by the subsystem in the state s 
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Environment 



Subsystem 



b 





b 



b 



Figure 4.1: The two processes of emission (on the right) and of absorption 
(on the left) of energy from the environment to the subsystem are schemat- 
ically represented, b and b are the states of the environment which can be 
either energy in the micro canonical case or temperature in the canonical 
case. 

and allows it to jump to the state s. This process contributes therefore 
positively to the dynamics of Pssi^, t). Notice that this term is different from 
the positive contribution term of the Redfield population equation 
applied to microcanonical environment correlation functions. In fact, in our 
equation, this terms takes into account the fact that the environment has 
an energy e — fiiOss before absorbing the quantum htOss- This is consistent 
with the principle of conservation of the total system energy which states 
that the energy gained (respectively lost) by the environment has to be lost 
(respectively gained) by the subsystem. The Redfield population equation, 
if applied to microcanonical environments, assumes that the environment 
is always at the same energy e and is therefore in contradiction with the 
total system energy conservation. However, this violation can sometimes be 
neglected as we will see in the next section. The absorption and emission 
processes are schematically represented in figure UTTl 

We finally notice that equation 1)4. 32() can satisfy detailed balance. If Pssi^) 
denotes the stationary distribution, the detail balance condition reads 



a(e, -uJss)P!!{e) = a{e - /iwgs, Wss)-P||(e - hujss)- 



(4.33) 



4.2 Markovian approximation 
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Therefore, if the detailed balance condition is satisfied, we get [using the 
property (|A.33j)] 

"^^^ (4.34) 



-hwss) n{e-tuvss)' 

We will see in the next section that this relation leads to the subsystem 
Boltzmann distribution for large heat capacity environments. 

4.2.2 Coherence dynamics 

To study the coherence dynamics, we isolate one coherence and neglect its 
coupling to the other coherences and populations in equation 1)4. 31() . We get 
an equation of the form 

Pss'{e,t) ^ {-iooss' -lss')Pss'{e,t) (4.35) 

with the frequency 

= ^A_M, (4.36) 

where 

Es = E, + ^{e\B\e)S,s + l'^-"!' / dioV^^^ (4.37) 
and with the damping rate 

7s,s' = -^-^(5ss-5.'s')'a(e,0) (4.38) 

s^s 

The environment interaction has two kind of effects on the coherence dy- 
namics. A Lamb shift effect which modifies the free coherences oscillations 
of the subsystem and a decoherence effect which irreversibly damps the co- 
herences. We can notice that the approximate coherences equation ()4.35p 
obtained from our new equation and the one (|3.115|) obtained from the 
Redfield equation are the same. However, the approximations performed 
to derive these two equations are crude and the real coherences dynamics 
predicted by both equation are not the same. 



70 



4. New kinetic equation for finite environment 



4.3 The infinite-environment limit 

In this section, we will investigate the limit of large environments of our 
new theory and we will show that our new equation simplifies to the Red- 
field equation in this limit. We will also discuss the greater generality of our 
new theory with respect to the Redfield theory. 

Let us start by integrating our new equation ()4.26|) over energy in order 
to obtain the reduced density matrix of the subsystem by using the fact that 
psit) = J deP{e, t). Making the change of variable of the energy e ^ e' — huj 
for the two last order terms, the energy integrated equation becomes 



No assumption has been made till now. If the environment is large compared 
to the subsystem, on can assume that the environment quantities {e\B^\e) 
and aKK'(e,w) do not vary significantly on energy scales of the order of the 
subsystem energy: 



where A represents the typical energy scale implied in the subsystem dy- 
namics. Therefore, if the environment is initially in a microcanonical state 
at energy e, the assumption H4.40() means that it will stay in this state with- 
out being affected by the subsystem emitted or absorbed quanta. Equation 




(e± A|5«|e± A) 
otKK'i^ ± A, a;) 




(4.40) 
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H4.39() can therefore be written as 

Ps{t) = tsMt)-lx{e\B^\e)Y,[S.,Ps{t)] + ^J2f^^J^''i 

-al,,{e,u)e'--ps{t)S^'{-T)S^ 

+d^,,ie,u;)e-''^^S^,i-T)psit)S, }. (4.41) 

We see that this equation is the Redfield equation for a microcanonical en- 
vironment. 

The greater generahty of our new equation compared to the Redfield 
equation can now be physically understood. In our new equation, the state 
of the environment is correlated to the state of the subsystem. The quanta 
which are emitted from the subsystem to the environment an vice versa 
have an effect on the subsystem as well as on the environment. In the 
Redfield equation, the Fourier transform of the environment microcanonical 
correlation functions which appears in the equation is always evaluated at 
the same energy. This means that the emitted or absorbed quanta have 
no effect on the environment in the Redfield equation. The only way to 
justify this assumption is to suppose that the environment is large so that 
its energy-dependent quantities vary on energy scales which are much larger 
than the typical energies emitted or absorbed by the subsystem. This is pre- 
cisely the meaning of the approximation (|4.4fl|) that we made to reduce our 
new equation to the Redfield equation. This point is particularly explicit 
if one compares the population dynamics obtained from our new equation 
1)4. 32(1 to the one obtained from the Redfield equation H3.111|) . If one applies 
both equations to a subsystem interacting with a small environment, our 
new equation takes into account the energy conservation of the total sys- 
tem (subsystem plus environment) although the Redfield equation does not 
(see figure . This shows that our new equation is more general than the 
Redfield equation and more pertinent to describe small total systems such 
as nanosystems. 

Another way to see the problem of energy conservation if the Redfield 
equation is applied to a finite environment, is to consider the expectation 
value of the total energy. When the total system evolves, the energy of the 
subsystem, the energy of the environment and the interaction energy vary 
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4. New kinetic equation for finite environment 




Figure 4.2: Representation of the energy exchanges between the environ- 
ment and the subsystem described respectively by the Redfield population 
equation (|3.111j) and by our new energy distributed (e.d.) population equa- 
tion (|4.32|) in the case of a two-level system. The subsystem energy Egys 
increases along the y-axis and the environment energy -Benw along the x-axis. 
The total energy of the system is represented by Etot- The initial condition 
corresponding to the subsystem at energy Eg and the environment at en- 
ergy e is denoted by the full black dot in the (x, y)-plane. The transitions 
that preserve the energy of the total system have to occur along the total 
energy line E = Eg + Ef, (the thick full line) in the (x, y)-plane. One can see 
that only the transition described by our energy distributed (e.d.) equation 
(dashed line) satisfies this condition. The Redfield equation describes tran- 
sitions (dotted line) that occur along a constant environment energy line 
and is therefore erroneous when the subsystem energies are of the order or 
larger than the typical energy scale on which the Fourier transform of the 
correlation function of the environment varies. 
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in time while keeping the total energy constant 

{H) = {Hs)t + {HB)t + X{SB)t (4.42) 

If the coupling parameter is small, the contribution of the interaction energy 
to the total energy can be neglected 

{H) {Hs)t + {HB)t. (4.43) 

In the Redfield theory, the environment does not evolve and the total density 
matrix is assumed to have the form 

PT{t) = ps{t)PB- (4.44) 

Therefore the subsystem energy evolves, but the environment energy does 
not. 

{H)^='Y,Es{s\psit)\s)+ I den{e)e{e\p'^\e) (4.45) 

s 

In order to avoid a violation of the total energy conservation, one has to 
suppose 

{H) « {Hb) » {Hs)u MSB)t. (4.46) 
In our new theory, the total density matrix is assumed to have the form 

prit) = (4.47) 
n{HB) 

Therefore the total system energy is kept constant but energy is transferred 
between the subsystem and the environment 

{H) Y.^- [ {e,t)+ f deeY,Pss{e,t). (4.48) 

s s 

This theory is valid if 

(H) » X{SB)t, (4.49) 

which is much less restrictive than the condition (|4.46|) . This again shows 
the larger range of validity of our theory compared to the Redfield theory. 
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One can also show that our new theory ahows a consistent understanding 
of the emergence of a canonical Boltzmann equilibrium distribution of the 
subsystem when the environment in a microcanonical state becomes large. 
The temperature of the Boltzmann distribution is then given by the micro- 
canonical temperature of the environment. In fact, in our new theory, the 
equilibrium distribution (which is obtained by using the detailed balance 
condition) is given by equation ()4.34|) and has a simple form because of the 
property found in equation ()A.33|) . Saying that the environment is large 
enough so that its energy-dependent quantities do not vary on energy scales 
of the order of the typical subsystem energies is equivalent to saying that the 
environment heat capacity has to be large. Using the results of appendix IbI 
for systems with a large heat capacity, equation ()B.7|) can be written as 

n{e - huJss) ~ n(e)e~^™-(')'^'^^ . (4.50) 

where Pmici^) = ^/{kBTmici^))- Therefore the detailed-balance condition 
given by equation (|4..S4j) becomes 



PjKe - nwss) 

Integrating equation 1)4. 51() over energy and using the property HB.8|) . we get 

{s\ps^\s) = j deef^^^^^'^^'^Pllit-hj-ss) (4.52) 

« e^— (^)^^"^(s|/5f |s) (4.53) 

Therefore, we find that the subsystem equilibrium distribution is given by 
the Boltzmann distribution at the microcanonical temperature of the envi- 
ronment 

/ I "^''1 \ 

\^\PS ^ g/3mic(e)fe^ss_ (4.54) 



Here, we have shown how, in our new kinetic theory, an equilibrium canon- 
ical Boltzmann distribution can arise on a subsystem in contact with an 
environment in a microcanonical state if the environment has a large heat 
capacity. This treatment cannot be made in the Redfield theory. In fact, 
in the Redfield theory, the microcanonical condition ()A.33|) cannot be used 
to simplify the subsystem equilibrium distribution 1)3. 114(1 . This indicates 
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that, contrary to our new equation, the Redfield equation apphed to a micro- 
canonical environment contains non-physical features which make the mi- 
crocanonical environment description inconsistent with the canonical equi- 
librium distribution of the subsystem. If the total system (subsystem plus 
environment) is an isolated finite system, our new theory becomes essen- 
tial. When the environment becomes very large, our theory reduces to the 
Redfield theory. Numerical evidence for this will be obtained in chapter [21 
devoted to the spin-GORM model. 

A summary of the different equations that we discussed is compiled in 
figure 14. HI 

4.4 Application to a two-level subsystem 

We now apply our new theory to describe the dynamics of a two-level subsys- 
tem interacting in a non-diagonal way with a general environment. Thanks 
to a simplification, the equations of our theory take in this case a sufficiently 
simple form to allow analytical solutions and an easy identification of the 
physical processes in play. This simplification consists in neglecting the con- 
tribution of the diagonal elements of the environment coupling operator in 
the basis diagonalizing the Hamiltonian of the environment. It should be 
noticed that this simplification can be the feature of many systems. For 
instance, it is exactly satisfied for any random matrix model in which the 
environment coupling operator B has a vanishing ensemble average. In this 
case, the average of the matrix elements (e|i?|e) over the ensemble of random 
matrices B also vanishes for all eigenvectors |e) of Hb- 

The Hamiltonian of a two-level subsystem interacting in a non-diagonal 
way with an environment is given by 




(4.55) 



where 

• ^x, o'y, and (Tz are the 2x2 Pauli matrices, 

• y(5"2 is the Hamiltonian of the two- level subsystem, 

• A is the energy spacing between the two levels of the subsystem, 

• Hb is the Hamiltonian of the environment, 

• (Jx is the coupling operator of the subsystem, 

• B is the coupling operator of the environment. 
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von Neumann 



NM energy-distributed eq. 



M energy-distributed eq. 




infinite environment 



infinite environment 



NM Redfield 



M Redfield 



Coherence Redfield 




Population 
energy-distributed eq. 



infinite environment 



Population Redfield 



Figure 4.3: Tree of possible approximations starting from the von Neumann 
equation to derive a quantum kinetic equation for the dynamics of a sub- 
system weakly interacting with an environment. The most general of these 
kinetic equations is the non-Markovian version of our new energy distributed 
equation (|4.26)) ("NM energy distributed eq."). By performing the Marko- 
vian approximation on this equation, one gets the Markovian version of it 
("M energy distributed eq.") which is valid on time scales longer than tc (the 
typical environment correlation function time scale). If instead, the infinite 
environment limit is performed on the non-Markovian version of our new 
equation, we get the non-Markovian Redfield equation (|3.9(Sj) ("NM Red- 
field"). If one further performs the Markovian approximation on this last 
equation, one gets the Markovian Redfield equation ("M Redfield"), which 
is equivalently obtained by performing the infinite-environment limit on the 
Markovian version of our new equation. The energy distributed population 
equation (|4.32l) ("Population energy distributed eq.") can be derived from 
the Markovian version of our new equation by neglecting the coherences 
contribution to the population dynamics. The Redfield population equation 
p.lllj) ("Population Redfield") is obtained either by the same procedure 
from the Markovian Redfield equation or by performing the infinite envi- 
ronment limit on the energy distributed population equation. Finally, an 
approximated coherence equation [ (|4.3<S|) or (|3.115j) ] ("Coherence Redfield") 
is obtained from the Markovian version of our new equation, as well as from 
the Markovian Redfield equation, by neglecting the other coherences and 
the population contribution to the selected coherence dynamics. 
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• A is the coupling parameter between the subsystem and the environment. 
The subsystem eigenvalue equation reads 

Hs\s) = ^a,\s) = s^\s), (4.56) 

where s = ±1. 



The two-level subsystem density matrix is, in our new theory, distributed 
in energy as follows 

The energy-distributed elements of the subsystem density matrix evolve ac- 
cording to our new non-Mar kovian equation 1)4. 26() . In the subsystem eigen- 
basis, equation (|4.29|) is given in the present case by 

P±±{e,t) = -'-X{e\B\e){P^±{e,t) - P±^{e,t)) (4.58) 



A2 



dr J duj { 

-P±±{e, t)a{e, -w)(e^(±^+Mr + ^-^(iA+Mr) 
+P^^{e - nuj,t)a{e - hw,uj){ei^^^+^^^ + e-^(±^+^)-) } 
for the distributed population and by 

P±^{e.t) = T'j:AP±^{e,t) - ^X{e\B\e){P^:f{e,t) - P±±{e,t)) (4.59) 

-P±^ (e, t)a{e, w)(ei(±^-^)^ + ei(±^+'^)") 

for the coherences. Performing the Mar kovian approximation on equation 
(|4.58j) and (|4.59|) or, equivalently, using the Markovian equation (|4.31|) . we 

get 

P±±{e,t) = -lA(e|S|e)(P^±(e,t)-P±^(e,t)) (4.60) 
-2Trj^a{e,±A/h)P±±{e,t) 

\2 

+27r-^d(e ± A, ^A/h)P^^ (e ± A, t) 
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for the population and 

P±^(e,t) = T^AP±^(e,i)-^A(e|^|e)(P^^(e,t)-P±±(e,t))(4.61) 



A A^ /" , ^a(e — hio,Lo) „ , ^ ^ 



A^ 

-7r-^P±T(e, t)(a(e, ±A) + a{e, ^A)) 
A^ 

+vr-^(P:^±(e T A, t)a(e ^ A, ±A/h) 

+P^±(e ± A, t)d{e ± A, tA/^i)) 

for the coherences. We see that the population and the coherence dynamics 
are coupled together by the first-order term in A in equation (|4.6flj) and 
()4.61|) . As mentioned before, to simplify the analysis, we will decouple them 
by assuming that the coupling operator of the environment is diagonal in 
the eigenbasis representation of the environment Hamiltonian : (e|i?|e) = 
for all e. Doing this, all the first-order terms in A in equation 1)4. 6U() and 
(|4.HH) disappear. Let us now define 

de'Z{e',t) (4.62) 
de'X{e',t) 



z{t) 


= Trp(t)(T^ 


x{t) 


= Trp(t)(T^ 


y{t) 


= TTp{t)ay 



de'Y{e',t), 



where 



C(e,t) ^ P++(e,t) + P_(e + A,t), (4.63) 

Z{e,t) = P++(e,t)-P__(e + A,t), 

X{e,t) = P+_(e,t) + P_+(e,t), 

Y{e,t) = iP+^{e,t)-iP^+{e,t). 

Using equation (|4.6fl)) and (|4.61|) . we find that the population dynamics can 
be studied with the variable Z(e, t) 

Z{e,t) = -2TTj^(d{e,A/h)-a{e + A,-A/h))Cie,t) (4.64) 
-27r-^(d(e, A/h) + d(e + A, -A/h})Z{e, t), 
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and Y{e,t) 



and the coherence dynamics with X{e,t) 

n^.t) = -^(^^^2^^|^vJ^^jYi,,t) (4.65) 

-Tr^{a{e,A/h) + a{e,-A/h))X{e,t) 
+7r-^a(€ - A, A/h)X{e - A, t) 
+7r-^a(e + A, - A/n)X(e + A, t) 

Y{e,t) = t(l + 2A^/^'-^^^r^)^(^'^) (4-66) 

-7r-p(a(e,A/n) + a(e,-A/n))y(e,t) 
-7r^a(e - A, A/;i)y(e - A, t) 

-7r-^a(e + A, -A/n)y(e + A, t). 

We also find that C(e,t) = is a constant of motion. If one chooses the 
initial condition of the total system to be the product of a general subsys- 
tem density matrix with a microcanonical distribution at energy e for the 
reservoir, we have that 

m = p5(0)fe^ (4.67) 
n{HB) 

P{Hb,0) 

n{HB) 

and, therefore, the reduced density matrix of the subsystem reads 

ps{0) = TrBp(O) (4.68) 
= J de'ps{0)S{e - e') 

= J de'P{e',0)5{e-e'). 
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We conclude that psi^) 
Z{e',0) 



Z{e' -A,0) 

X{e',0) 
Y{e',0) 

If the initial condition of the coherences is strictly "local" in the energy of 
the environment (i.e. a true delta distribution) and if the Markovian approx- 
imation is exact, the non-local terms in the coherence evolution equations 
(j4.65|) and (|4.66j) do not contribute to the coherence dynamics and can be 
neglected. However, the initial condition is never strictly "local" in the real 
world. Furthermore, the non-Mar kovian initial effects also contribute to 
give a width to the distribution. Therefore, even if these non-local effects 
are generally small, neglecting these terms in the coherence evolution has 
to be considered as an approximation in order to get analytical solutions 
for the subsystem dynamics. Notice that the terms we neglect are non-local 
only if A is not too small compared to the typical energy scale of variation 
of the micro canonical correlation function a(e it A,u;) ~ a{e,uj). If A is 
not too small, our assumption is not valid and our new equation reduces to 
the Redfield equation (see appendix [n]) . However, for the population, the 
non-local character of the dynamics is strong and cannot be neglected. This 
is precisely the improvement of our new theory with respect to the Red- 
field theory. Performing the local approximation on the coherence dynamics 
()4.65|) and (|4.66|) . the dynamics of the populations and of the coherences 
can be rewritten 

Zie,t) = 7pop(e)(Z(e,oo)-Z(e,t)) (4.70) 

X{e,t) ^ -(A/^ + r(e))y(e,t)-7coh(e)^(e,i) 
Y{e,t) « {A/h + T{e))X{e,t)-j,oi.{e)Y{e,t), 



= P(e, 0). This means that 



= C(e',0) = (+|ps(0)|+)(5(e-e') (4.69) 
= _c(e'-A,0) = -(-|/35(0)|-)5(6-6') 



= {{+\psm-) + {-\psm+m^ - 

= x(0)(5(e-e') 

= ii{+\Psm\-)-{-\psm+me-e') 

= y{0)6{e-e'). 
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where 



\2 

7pop(e) = 2Trj^{a{e,A/h)+a{e + A,-A/h)) (4.71) 
7coh(e) = TT-^{a{e,A/h) + a{e,-A/h)) 



r(e) = 2^ / diuV 



A2 - (;2u;)2 
7l)_ 

n(e) - n(e + A) 



- a{e + A,-A/h)+a{e,A/hf^''^> 



n(e) + n(e + A) ^ ^ 
The solutions of the equation (|4.7()j) are given by 

C(6,t) = C(6,0) (4.72) 
Z(e,t) = Z(e,oo) + [Z(e,0) - Z(e,oo)]e-^p°p(')* 

X{e,t) w {X(e,0)cos(A/;i + r(e))t-y(e,0)sin(A/;i + r(e))t}e-^-h(e)t 
y(e,t) w {X(e,0)sin(A/;i + r(e))t + y(e,0)cos(A/;i + r(e))t}e-^-h(^)*. 

Using (|4.62|) . we obtain the final result for the two-level subsystem reduced 
dynamics which reads 



zit) 


= Z{€,t) + Z{e- A,t) 


(4.73) 


x{t) 


= X{e,t) 




y{t) 


= Y{e,t). 





The two-level subsystem Markovian dynamics obeys a biexponential relax- 
ation of the populations and an oscillatory exponential relaxation for the 
coherences. 

Our results are consistent with the expected result at equilibrium for a large 
environment. In fact, the coherences disappear on long time scales and 
the populations reach an equilibrium value which is given by the canoni- 
cal Boltzmann distribution. This is shown using the property (|4.5fl|) which 
becomes in case of the two-level subsystem 



n(e- A) « n(e)e-^™-(^)'^. 



(4.74) 
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Using now the fact that C(e, 0)+C(e— A, 0) = 1, we find that the equilibrium 
value of the population becomes 

^(oo) = Z(e,oo) + Z(e - A,oo) (4.75) 

g-/3mic(€)A/2 _ g/3„ic(e)A/2 
~ e-/3mic(«)A/2 _^ g/3„ic(e)A/2 

/3mic(e)A 

= coth . 

2 

This is precisely the canonical expectation value of the population at a 
temperature given by the microcanonical temperature of the environment. 
Once again, we see that our new theory recovers the expected results from 
statistical mechanics when the environment becomes large enough. 
The Redfield equation applied to the two-level system is derived in appendix 
ini One can notice that, in general, the solutions predicted by our new 
theory and those predicted by the Redfield theory are not the same. The 
population dynamics is very different. The relaxation rate as well as the 
equilibrium value is different. But the coherence dynamics is similar. Both 
theories predict an exponential damping of the oscillating coherences at the 
same rate but the frequency of the oscillations differs only by the Lamb 
shift. One can assume for a very large environment that the microcanonical 
environment correlation function is almost constant on energies of the order 
of the subsystem typical energies around its microcanonical energy a{e it 
A, uj) « a(e, uj). As we have seen in section in this case our new equation 
reduces to the Redfield equation. In fact our population equation [z{t) in 
equation ()4.73|) ] reduces in this case to the Redfield population dynamics 
[z{t) in equation ()(L9|) ]. However, we had to do the non-local approximation 
to derive an analytical expression for the coherence dynamics. But this 
approximation fails when the environment becomes too large. Therefore, in 
this case, our approximated coherence equations [x{t) and y(t) in equation 
()4.73|) ] have to be replaced by the Redfield coherence equations [x{t) and 
y{t) in equation ()('.9|) ]. 
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Figure 4.4: Representation of the energy exchanges between the environment 
and the two-level subsystem described by the population equation (|4.73|1 . 
The subsystem energy Egyg increases along the y-axis and the environment 
energy Eenv along the x-axis. The initial condition is denoted by the two 
full black dots in the {x, y)-plane and is the product of a general two-level 
subsystem density matrix with the environment in a microcanonical state 
at energy e. The transitions due to the energy exchanges have to occur 
along the two thick black lines. These two lines correspond to the region 
of the plane preserving the quantity C(e — A,0) respectively C(e, 0) which 
are constants of motion. The dynamics of the population along these lines 
is described respectively by Z{e — A,t) and Z{e,t) of equation (|4.73j) . 
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Chapter 5 

The spin-GORM model 



In the previous chapter, we have presented a weak-couphng kinetic theory 
extending the Redfield theory to subsystems interacting with a finite en- 
vironment. The conditions of vahdity of these quantum kinetic equations 
remain however httle known. It is indeed very important to understand what 
are the conditions (in particular the minimum size) under which a quantum 
system can display an irreversible relaxation towards an equilibrium state. 
It is the purpose of this chapter to contribute to the clarification of these 
questions by studying the simple model of a two-level system or spin^ in- 
teracting with a random-matrix environment. Indeed, work done during 
the last decade has shown that typical quantum systems with a classical 
chaotic dynamics presents properties of random matrices on small energy 
scales (71 1351 I87j . Here, we consider the Hamiltonian of the environment 
as well as the operator of coupling between the spin and the environment 
to be given by random matrices taken in a statistical ensemble of Gaus- 
sian orthogonal random matrices. This defines a model that we called the 
spin-GORM model I24j . Its total Hamiltonian can be numerically diago- 
nalized and, therefore, its physical properties can be computed numerically. 
Furthermore, for this model many results can be obtained analytically. It is 
therefore a prototype model to test and to compare the different kinetic the- 
ories. Similar models using Gaussian orthogonal random matrices [321 154j 
or using banded random matrices El EH lEH have been studied. It 
has been shown that the dynamics of a subsystem interacting with an en- 
vironment trough an interaction term modeled by random matrices present 
common features with other types of interaction terms found in the litera- 
ture [57 . 

^Spin means in what follows a spin 1 /2 and is therefore equivalent to a two-level system. 
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5. The spin-GORM model 



The spin-GORM model differs from the spin-boson model by the density 
of states of the environment. Instead of monotonously increasing with energy 
as in the spin-boson model, the density of states of the environment obeys 
Wigner's semicircle law in our model and is thus limited to an interval of 
energy with a maximum density in between. We notice that such densities of 
states appear in systems where a spin is coupled to other (possibly dissimilar) 
spins as, for instance, in NMR, in which case the density of states of the 
other spins forming the environment also present a maximum instead of a 
monotonous increase with energy. Our model may therefore constitute a 
simplification of such kind of interacting spin systems. Our main purpose 
is to understand the conditions under which a kinetic description can be 
used in order to understand the relaxation of the spin under the effect of 
the coupling with the rest of the system, which we refer to as a complex 
environment. 

Our study of the spin-GORM model is separated in four parts. In section 
15.11 we define the model. In section 15.21 we study the spectral properties of 
the total system. In section [5^ we study the different regimes of the spin 
relaxation dynamics. Finally in section [5.41 we investigate the equilibrium 
distributions of the subsystem and of the total system. 



5.1 The model 

We are interested in the study of a total system composed of a simple sub- 
system (with a few discrete levels) interacting with a complex environment 
(with many levels). A prototype of such system is given by a two-level sub- 
system interacting with an environment defined in terms of random matrices. 
The von Neumann equation of this system reads 

dt n 
where the total Hamiltonian is given by 

H = ja, + flB + Aa^i (5.2) 

where 

• o'x, d'y, and dz are the 2x2 Pauli matrices, 

• Y'^'z is the Hamiltonian of the two- level subsystem, 

• A is the energy spacing between the two-levels of the subsystem, 

• Hb is the Hamiltonian of the environment, 
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• dx is the coupling operator of the subsystem, 

• B is the couphng operator of the environment, 

• A is the couphng parameter between the subsystem and the environment. 

The well-known spin-boson model |3Ul 1551 l9Uj is a particular case of this 
total system where corresponds to an infinite harmonic oscillator lattice 
and B is linear in the degrees of freedom of the environment. Here, we want 
to define a new model, the spin-GORM model, also described by the total 
Hamiltonian (|5.2j) and for which Hb and B are Gaussian orthogonal random 
matrices (GORM) (see appendix^ for some basic properties of GORM). 

Such a model can represent for example an electronic spin traveling in 
a semiconducting quantum dot in the shape of a chaotic billiard and, fur- 
thermore, submitted to a magnetic field perpendicular to the plane of the 
dot. The quantum dot contains a localized ferromagnetic impurity with a 
classical magnetic dipole. The two spin levels of the traveling electron are 
split by the magnetic field and interact through a magnetic dipole-dipole 
interaction in a non-diagonal way with the classical magnetic dipole of the 
localized ferromagnetic impurity. Since the position and momentum opera- 
tors of the electron are associated with a classically chaotic dynamics, they 
can be represented by Gaussian orthogonal random matrices according to 
the Bohigas-Giannoni-Schmit conjecture "T. This system is schematically 
depicted in figure [5Tl 




Figure 5.1: Schematic representation of an electronic spin in a chaotic bil- 
liard quantum dot with a magnetic field perpendicular to the plane of the 
dot interacting through magnetic dipole-dipole potential with a ferromag- 
netic impurity described by a classical spin. This physical system can be 
modeled by a Hamiltonian similar to the spin-GORM model. 

Let us now discuss the random matrix aspect of the spin-GORM model 
in more detail. As we said, we want to model a two-level subsystem that 
interacts with an environment that has a complex dynamics. Here, complex 
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is used in a generic way. The complexity can come, for example, from the 
fact that the corresponding classical system is chaotic j3 OSl OHl EZI like in 
a quantum billiard or for the hydrogen atom in a strong magnetic field. It 
can also come from large number of couplings in an interacting many-body 
system as in nuclear physics |35| or in interacting fermion systems such as 
quantum computers jSHlEl- In the sixties, Wigner O ESI ESI the first to 
develop random-matrix theory for the purpose of modeling spectral fluctua- 
tions of complex quantum systems containing many states interacting with 
each other. This tool has now become very common in many fields from nu- 
clear physics to quantum chaos. This is the reason why we consider random 
matrices to characterize the complexity of the environment operators. The 
environment operators of the spin-GORM model, Hb and B, are defined by 



B = a^j,X , (5.3) 



where X and X' are two different Gaussian orthogonal random matri- 

ces with mean zero. Their non-diagonal (diagonal) elements have standard 

deviation crj^^, = 1 (fx^ = \/2)- X and X are two different realizations of 
the same random-matrix ensemble and have therefore the same statistical 
properties, cr^f) and are the standard deviations of the non-diagonal 

elements of Hb and B, respectively. For these random matrices, the width 
of their averaged smoothed density of states is given by 



VHb = a^f,Vm, 

VB = a^aVm (5.4) 



(see Appendix Id]). 

It is interesting to define the model in such a way that, when N is increased, 
the averaged smoothed density of states of the environment increases with- 
out changing its width VHb- The width can be fixed to unity. This is 
equivalent to fixing the characteristic time scale of the environment. For 
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doing this, it is necessary to rescale the parameters as follows: 

' a = a^fjVm , 
t = at , 



A = A^, 



The von Neumann equation of the spin-GORM model becomes 

dm 



dt h 
with the rescaled total Hamiltonian 



[H,m], 



(5.5) 



H = Hs + HB + \a.^B 
A 1 ^ 



-.X . 



(5.6) 



As announced, we have now VHb = T^B = 1. 

In what follows, without loss of generality, a will always be taken equal to 
unity. Notice that, to model an environment with a quasi-continuous spec- 
trum, the random matrix must be very large (A^ — > oo). 
In order to get ensemble averaged results, one has to perform averages over 
the different results obtained for each realization of ()5.6|) . When we use 
finite ensemble averages, the number of members of the ensemble average 
will be denoted by x- 

We see that the Hamiltonian ()5.6|) is characterized by three different pa- 
rameters: A, A, and A^. We define three different domains in the reduced 
parameter space corresponding to a fixed A^ in order to facilitate the follow- 
ing discussion. These three regimes are represented in figure EI21 domain A 
with 1 > A, A; domain B with A > 1, A; domain C with A > 1, A. 



5.2 Spectral analysis 

In this section, we study the spectrum of the complete system for the differ- 
ent values of the parameters. This study is important in order to understand 
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A 



1 
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A,,1<A 




A,l<^ 



domain C 



k,A <1 



domain A 







X 







1 



Figure 5.2: Representation of the three different domains in the space of the 
reduced parameters A and A of the model for a fixed number N of states. 

the different dynamical behaviors that we encounter in the model. 

Let us begin by defining the notations in the simple case where there is 
no coupling between the two parts of the total system (A — ^ 0). The isolated 
subsystem has two levels separated by the energy A: 



where s = ±1. The environment has the standard spectrum of a GORM 



where b = 1,2,..., N/2. The Hamiltonian of the total system without inter- 
action between the subsystem and the environment is thus 



(5.7) 



HB\b) = E^\b), 



(5.8) 



Ho = Hs + Hb 



(5.9) 



and the spectrum is therefore given by 



(5.10) 



with n= 1,2, N and 




(5.11) 
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The eigenvectors are tensorial products of both the subsystem and environ- 
ment eigenvectors: 

|n) = \s)\b). (5.12) 

Let us now define the notations in the opposite case, when tlie coupling term 
is so large that the Hamiltonian of the subsystem and of the environment 
can both be neglected (A — oo). Using the unitary matrix U acting only on 
the subsystem degrees of freedom 



U 



1 1 ■ 

1 1 
V2 



the total Hamiltonian becomes 

^0 = Aa^S. (5.13) 

E^n and {ki]) = \k) <Si \rj) are, respectively, the eigenvalues and eigenvectors 
of the Hamiltonian 

Ho\Kri) = XE^r,\i^ri) = \KEr^\Kri), (5.14) 
where 77 = 1, -y and k = ±1. 

After having defined the notation in the two extreme cases A ^ and 
A ^ 00, we start the study of the spectrum with interaction A 7^ 0. The 
total spectrum is given by the eigenvalues {E^} that are solutions of the 
eigenvalue problem: 

H\a)=E^\a), (5.15) 

where a = 1, 2, A^. It is very difficult to obtain analytical results for this 
problem. We will therefore study the total spectrum using a method of 
numerical diagonalization of the total Hamiltonian. 



5.2.1 Smoothed density of states 

In order to have a quantitative understanding of the global aspect of the 
spectrum (on large energy scales), we will study the total perturbed aver- 
aged smoothed density of states. 

The environment-averaged smoothed density of states obeys the semi- 
circle Wigner law (see equation 
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where e is the continuous variable corresponding to the environment energy 
. Therefore, when A = 0, the total averaged smoothed density of states 
is the sum of the two environment semicircular densities of states which 
correspond to both states of the two- level system (see equation (|5.1fljl ): 



n{e) = n-(e--) + n-(e + -) 

if (i-#)<N<(i + #)>(5-17) 
elsewhere 



where e is the continuous variable corresponding to the total energy E^. The 
semicircular densities n'"(e — y) and n'^{e + ^) are schematically depicted 
in figure 1^31 for different values of A. The numerical density of states of the 
total system corresponding to A = are depicted in figure |S3t a)- (c). 
When A — > oo (meaning that the coupling term becomes dominant in the 
Hamiltonian) , the averaged smoothed density of states of the total system 
(see equation (|5.14|) ) is given by 



n(e)=n-(i) + n-(-i) 







if 
if 



|e| < I 
lel>4 



where e is the continuous variable corresponding to the total energy E^t]- 
This result can be observed in figure EUfa) for A = 10 (because A, 1 <C A). 
When A 7^ 0, the total averaged smoothed density of states is also plotted 
in figure 15.41 The main observation is that there is a broadening of the 
complete spectrum when one increases A. In figure IKTif a). we see that the 
averaged smoothed density of states changes in a smooth way from (|5.17j) 
to H5.18() . But in figure I^Tlf bl and much more in figures [Ofc) and (d), the 
two semicircular densities n"'(e — y) and n'"(e + y) seem to repel each other 
as A increases. This is due to the fact that the levels of a given semicircular 
density do not interact with each other but only interact with the levels of 
the other semicircular density. This is a consequence of the non-diagonal 
form of the coupling. Therefore, having in mind the perturbative expression 
of the energies (see equation ()E.6|) of Appendix 0, one understands that 
when A is nonzero, the eigenvalues that are repelling each other with the 
most efficiency are the ones closest to the center of the total spectrum. 



5.2 Spectral analysis 



93 



Smoothed density of state 



N/2 levels 




Energy 



A«l 





Energy 



Energy 




A»l 




Energy 



Figure 5.3: Smoothed densities of states n'^{e 



and n^{e + 4) for 



different values of A. The total smoothed averaged density of states of 
the non-perturbed spectrum is obtained by the sum of them (see equation 



5.2.2 Eigenvalue diagrams 

The global effect of the increase of A on the eigenvalues has been studied 
with the average smoothed density of states. But in order to have an idea 
of what happens on a finer energy scale inside the spectrum, it is interesting 
to individually follow each eigenvalue Ea as a function of A on an eigenvalue 
diagram. 



The first thing to note (see figure 1^3)) is that the increase of the coupling 
induces a repulsion between the eigenvalues. The result is the broadening 
of the total spectrum as we already noticed on the smoothed averaged den- 
sity of states. If one looks closer inside the fine structure of the eigenvalue 
spectrum, we see that there is no crossing between the eigenvalues. This 
is a consequence of the fact that there is no symmetry in the total system. 
Therefore, the non-crossing rule is always working. Each time two eigenval- 
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Figure 5.4: Total smoothed averaged density of states obtained numerically 
for different values of A and A. (a) corresponds to A = 0.01, (b) to A = 0.5, 
and (c) and (d) to A = 5. For all of them A'' = 500 and x = 50. Notice that 
the A = and the A = 0.1 curves are not distinguishable. 
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ues come close to each other, they repel each other and create an avoided 
crossing. One can notice that there is a large number of avoided crossings 
inside the region where the two semicircles n^{e— -y) and n"'(e + overlap 
(see figure E21 ill order to visualize the overlapping zone that extends from 
y — 1 to 1 — y). But outside this overlapping zone, there are very weak 
avoided crossings (and of course no crossing) and all the eigenvalues appear 
to move in a regular and smooth way. The regions seen in figures [S3t a), 
(b) and (d) are inside the overlapping zone and the one in figure EHI^) is 
outside. In figure 1^31 c) . the lower part of the diagram is inside the over- 
lapping zone and the upper part is outside. This phenomenon is due to the 
fact that the eigenvalues of a given semicircular density do not interact with 
the eigenvalues of their own semicircular density but only with those of the 
other semicircular density. It is the consequence of the non-diagonal nature 
of the coupling in the spin degrees of freedom. When the coupling becomes 
large enough (A > 1), the avoided crossings also disappear inside the over- 
lapping zone. Around A = 1 and inside the overlapping zone, there is a 
smooth transition from an avoided-crossing regime that gives rise to a tur- 
bulent and complex A evolution to another regime without much interaction 
between the eigenvalues that gives rise to a smooth A evolution. 

5.2.3 Spacing distribution 

The eigenvalue diagrams only give us a qualitative understanding of the fine 
energy structure of the spectrum. For a more quantitative study, it is inter- 
esting to look at the spacing distribution of the spectrum. 

When A = 0, each eigenspace s = it-j, corresponding to a different 
subsystem level, has a Wignerian level spacing distribution 

p«'(s) = |se~f^'. (5.18) 

But the total spectrum is obtained by the superposition with a shift A of 
two of such spectra. The superposition creates a Poissonian component to 
the total spacing distribution in the overlapping zone (y — Ito— y + 1). A 
Poissonian distribution is given by 

pp(s) = e-^ (5.19) 

Therefore, we choose to fit the total spacing distribution by the mixture 



pf^\s) = CiP'"{s) + (1 - Ci)PP{s). (5.20) 
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2 4 6 8 



X 

Figure 5.5: Different parts of the eigenvalue diagrams with N = 500, cor- 
responding to (a) and (b) A = 0.01, (c) and (d) A = 0.5, and (e) A = 5. 
They represent the eigenvalues of the total Hamiltonian 1)5. 15|) as a function 
of A. 
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This choice of the form of the fit is empirical but reasonable because the 
correlation coefficient of the fit is always close to one (between 0.963 and 
0.982). 

We computed the spacing distribution and made the fit (|5.2fl|) in order to 
compute the mixing coefficient Ci for different values of A in the case where 
the overlapping zone covers almost the whole spectrum (A <C 1). The results 
are plotted in figure EH We see that there is a specific region of A values 
where the total spacing distribution is close to a pure Wignerian one. This 
region corresponds to the situation where A^A^ = 0(1), i.e., when the typical 
intensity of the interaction between the non-perturbed levels [ which is O(A^) 
since the first order correction in perturbation theory is zero in our model ] 
becomes of the order of the mean level spacing 0{j^) in the total system. In 
this region, the level repulsion is maximal and effective among almost all the 
states. In the limit A — > oo, the spectrum is again the superposition of two 
semicircle distributions, one for the Ei^^^s and one for the i?_i^'s according 
to equation ()5.14|) . As a consequence, one gets a Poissonian type of spacing 
distribution. The other A cases will have a Wignerian spacing distribution 
outside the overlapping zone and a mixed one (like for the case A ^ 1) 
inside. This can be visually seen on the eigenvalue diagram that we studied 
before. 

5.2.4 The shape of the eigenstates (SOE) 

We want now to have some information about the eigenstates of the total 
system inside the overlapping zone of the two semicircular densities. 

Following ^1], we define the quantity 

e(e,e°) = ^\{a\n)fS{E^ - e)5{E^ - e°) (5.21) 

We call local density of states (LDOS) the function ^{e^e^) for fixed and 
varying e. 

We call shape of the eigenstates (SOE) the function ^{e,e^) for fixed e and 
varying . 

We here focus on the SOE. The SOE tells us how close a perturbed eigen- 
state (A / 0) at energy e is from the non-perturbed eigenstate (A = 0) at 
energy e*^. If the SOE is a very narrow function centered around e = , 
the concept of a non-perturbed eigenstate is still useful. This regime corre- 
sponds to very small coupling, for which the interaction intensity is lower 
than the mean level spacing < A^A^ < 1, and will be called to the localized 



98 



5. The spin-GORM model 




X 

Figure 5.6: (a) Spacing distribution for different values of the coupling pa- 
rameter A. (b) Fitted coefficient Ci of equation (|5.2fl|) . The closer is Ci 
from unity, the closer is the spacing distribution to the Wigner spacing dis- 
tribution. In the two figures, A = 0.01, = 500, and x = 50. 
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regime. In the limit N ^ oo, this regime disappears. If one increases the 
couphng, the interaction intensity between the non-perturbed states begins 
to be larger than the mean level spacing between the states: X'^N > 1. The 
non-perturbed levels start to be "mixed" by the interaction and the SOE 
starts then to have a Lorentzian shape with a finite width F, centered around 
e = e^. This regime is called the Lorentzian regime. If one further increases 
the coupling parameter, the SOE begins to spread over almost the whole 
spectrum. This regime is called the delocalized regime. 
In the banded random matrix model of 14.., the regimes are classified ac- 
cording to a different terminology and there is an additional regime cor- 
responding to a spreading that goes beyond the energy range where the 
coupling acts (due to the finite coupling range of the banded matrices) . The 
motivation for our change of terminology will become clear in the study of 
the dynamics. 

The Lorentzian regime can be separated into two parts. For small coupling, 
the width of the Lorentzian F is smaller or of the order of magnitude of 
the typical energy scale of variation of the averaged smoothed density of 
states of the environment 6e: n(e -|- 6e) ~ n{e). For larger coupling, the 
Lorentzian width extends on an energy scale larger than the typical energy 
scale of variation of the density of states of the environment. Therefore, one 
has F < 5e in the former case and T > 6e m the latter case. The different 
regimes are represented in figures EiZl and EiBl To represent the SOE of the 
spin-GORM model, we discretize the energy axis in small cells of the order 
of the mean level spacing and we average the SOE over x realizations of 
the random-matrix ensemble. We see in figureEHta) the typical shape of the 
SOE going from the perturbative regime (A = 0.01,0.05) to the beginning 
of the Lorentzian one (A = 0.1). In figure EHJb), we see the SOE across 
the Lorentzian regime (0.2 < A < 0.8). We also see the delocalized regime, 
when A ^ oo and the SOE gets completely flat (A = 10). Figure EU^c) 
shows the width of the Lorentzian obtained from a fit made on the SOE 
curve. The correlation coefficient of the fit helps us to determine the region 
of the Lorentzian regime where the SOE is very well fitted by a Lorentzian. 
It has been verified that the width of the Lorentzian is independent of N in 
the Lorentzian regime. Figure I5.9r d) (log-log) shows that the width of the 
Lorentzian has a power-law dependence in the coupling parameter close to 
two in the Lorentzian regime. 
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Figure 5.7: Diagram of the different regimes as a function of the coupling 
parameter A. 
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Figure 5.8: Schematic representation of the different regimes in the plane of 
the reduced parameter A and N. The Lorentzian 1 regime corresponds to 
r < (5e and the Lorentzian 2 regime to F > (5e. 
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Figure 5.9: (a) SOE in the perturbative regime, (b) SOE in the Lorentzian 
regime, (c) The width and the correlation coefficient of the fit of the SOE 
by a Lorentzian. (d) Power-law dependence of the width of the Lorentzian 
SOE in the coupling parameter. In all the figures: A = 0.01, N = 500, 
X = 50, and e = 0. 
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5.2.5 Asymptotic transition probability kernel (ATPK) 

An interesting quantity, which is close to the SOE, but which has a nice phys- 
ical interpretation, is the asymptotic transition probability kernel (ATPK). 

The transition probability kernel (TPK) gives the probability at time t 
to be in the level \n) if starting from p(0). It is defined as 

Ut{n\p{0)) = {n\e-'^^p{0)e'^*\n) (5.22) 
The ATPK is the time average of the TPK 

noo(n|p(0)) = hm i [ dmt{n\p{0)) 

= ^|(a|n)p(a|p(0)|a). (5.23) 

a 

The distribution of the ATPK in energy is given by 

noo(£°|p(0)) = ^ U^{n\p{0))6{En - e°). (5.24) 

n 

The ATPK has an intuitive physical interpretation. It represents the prob- 
ability after a very long time to end up in a non-perturbed state |n), having 
started from the initial condition p{0). The ATPK is the convolution 

noo(e°|e°') = ^ |(a|n)n(a|m)|25(^, - e'>)6{E^ - e°'), (5.25) 

n,m a 

With p(0) = Zm \m){m\6{Em - £°'). 

Since we are interested in a random matrix model and for the purpose of 
studying the dynamics, we will perform averages of two different kinds. 
The first kind of average is a microcanonical average over states belonging 
to the same given energy shell of width 5e for a given realization of the 
random-matrix ensemble used in our total Hamiltonian. The second kind 
of average is an ensemble average over the x different realizations of the 
random matrices ensemble. For the microcanonical average, a choice of the 
width of the energy shell de has to be done in such a way that it is large 
enough to contain many levels (to get a good statistics) and small enough 
to be smaller than or equal to the typical energy scale of variation of the 
averaged smoothed density of states of the environment 5e. Therefore, the 
adequate choice of the width of the energy shell corresponds to jr < Ss < 6e. 



0.8 



^ 0.6 
H 

< 0.4 



0.2 







(a) 



-0.4 



Nr=0.1 Pin=99.5% 

NX^=l Pin=94. 

N?i^=10 Pin=68. 

Pin=9.4% 
X^IO Pin=7.3% 



-0.2 







0.2 







0.4 



0.8 



0.6 



Oh 

H 

< 0.4 



0.2 



(b)' 



Nr=o.i 

NA.^=10 



-0.015 -GDI -D.005 



0.005 0.01 0.( 



!.02 



0.01 r 



(C) 



0.0075 - 



:0.005 - 



0.0025 r ~« 



-0.4 



-0.2 



0.2 



Nr=io 



0.4 



Figure 5.10: (a) ATPK from the localized regime, through the Lorentzian 
regime, to the delocalized regime. The ATPK is microcanonically averaged 
over the different initial conditions corresponding to the levels inside the 
energy shell centered at e'' =0 with width 5e^ = 0.05. "Pin" denotes the 
proportion of the ATPK that stays inside the initial energy shell after an 
infinite time, (b) and (c) ATPK for a single level as an initial condition, 
without any average. All the figures are obtained for very small system 
energy A = 0.01 with no random-matrix ensemble average x = 1 ^or 
N = 500. 
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Different ATPK are depicted in figure I^TTUl where there is no random-matrix 
ensemble average x = 1- In figure I^TUT a). the ATPK is a micro canonical 
average inside the energy shell at energy e'' and of width 5e^ . " Pin" denotes 
the total probability of staying inside the energy shell after a very long time. 
A^A^ = 1 is on the border between the localized and the Lorentzian regime 
and A = 1, 10 in the delocalized regime. We see that, until A^A^ = 1, the 
main probability stays in the initial energy shell. But when the Lorentzian 
regime starts, the probability spreads over energies larger than 6e (here 
6e^6e^'). In fi gures IKTinr b) and (c), no average has been done. The initial 
condition is a non-perturbed pure state (corresponding to an energy close 
to zero) and one can see the individual probability of being on another non- 
perturbed state for the different regimes. We can see that in the localized 
regime the probability of staying on the initial state is much more important 
than the probability of leaving it. We also see that the Lorentzian regime 
starts when the initial state loses its privileged position containing the main 
probability and, therefore, when the neighboring levels begin to have an 
important fraction of the total probability. 



5.3 Spin dynamics 

In this section, we study the different relaxation regimes of the spin in the 
spin-GORM model. This study relies on a detailed comparison between the 
numerical exact dynamics (subsection 15.3.3)) and the theoretical predictions 
in the weak-coupling limit (subsection l5.3.l]) and in the strong-coupling limit 
(subsection I5.3.2|) . The important question of the dynamics of a statistical 
ensemble versus the dynamics of a single element of the ensemble is investi- 
gated in subsection 15.3.41 

The exact spin dynamics in the spin-GORM model is obtained from 
the von Neumann equation (|5.5|) using the reduced density matrix psit) = 
TiBPit)- In this study, the total system has always a given finite and con- 
stant energy. We also always consider initial conditions for which the sub- 
system is in any state and the environment is in a microcanonical state at 
a given fixed energy e corresponding to a uniform distribution inside an en- 
ergy shell centered at e and of width 6e. 6e is chosen in such a way that it 
is large enough to contains many levels (to get a good statistics) and small 
enough to be smaller than or equal to the typical energy scale of variation 
of the averaged smoothed density of states of the environment 6e. There- 
fore, the adequate choice for the width of the energy shell corresponds to 
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< 6e < 6e. The dynamics is then averaged over the x reaUzations of the 
random-matrix ensemble. 

5.3.1 The weak-coupling regime (A ^ 1) 

In this subsection, we make the theoretical study of the dynamics of the 
spin, weakly interacting with its random-matrix environment. This means 
that we apply the weak-coupling theory that we developed in chapter |1] to 
the spin-GORM model. 

We derived in section ETTl of chapter a kinetic equation to describe the 
evolution of a subsystem (with a discrete spectrum) weakly interacting with 
its environment (with a quasi-continuous spectrum) and applied this equa- 
tion in sect Ion ian to a two-level subsystem interacting in a non-diagonal way 
with any environment. The spin-GORM model is in fact a specific case of 
the two-level subsystem model studied in section ^31 where the environment 
operators are Gaussian orthogonal random matrices as defined in section 
15.11 The central quantity characterizing the environment operators in the 
weak-coupling kinetic theories is the correlation function of the environment. 
If one is able to obtain a simple expression for the correlation function of the 
environment in the spin-GORM model, all the results of section I17H applv to 
this chapter by replacing the correlation functions by their random-matrix 
expressions. 

We therefore start by evaluating the correlation function of the environment 
and its Fourier transform in the spin-GORM model. Using the random- 
matrix environment operators that we defined in section 15.11 to write the 
spin-GORM Hamiltonian (|5.6|) . we find that the Fourier transform of the 
environment correlation function defined by equation ()A.19|) of appendix El 
is given by 

a{Eb,uj) = Y,6{uj + uJw)\Bw\^ (5.26) 

b' 

= 6{uj)\Bbb\^ + Y^5{uj + ujbb')\Bw\^ 

b'^b 

where cohh' = {E^ — Ei,i)/h. Averaging the Fourier transform of the environ- 
ment correlation function over the random-matrix ensemble, we get 

a{E,,u) = S{u;)^ + ^ h5{fuv + E, - E,,)^. (5.27) 
b'y^b 
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To simplify the notations, we will use the convention 

^=1 if x<0. ^^-28) 

Assuming that the environment has a spectrum that is dense enough to 
perform the quasi-continuous limit and using the smoothed density of states 
given by ()5.16|1 . we get 

ii(e.w) = i(„)_L + yrf/Ay/r7^j(Sj, + j_j') (5.29) 

This means that 



«(e,^)^=°°^y^-(6 + M^- (5.30) 



The environment correlation function can now be obtained by taking the 
inverse Fourier transform of (|5.3fl|) 

a{e,t) = J doje-''^^a{e,uj) (5.31) 

j de'i-yi - (e')2e-^^'*/V^*/'^. 

Performing the integral we get 



(e,t) A^z;oo -^i(;/(f)) ,..^A (5.32) 



Ji {z) is the Bessel function of the first kind. We insist on the fact that these 
two results suppose that the environment is quasi-continuous {N oo) and 
that the random-matrix ensemble average has been performed {x ~^ oo). 
Given these fundamental quantities, the results of section ^1] directly apply 



5.3 Spin dynamics 



107 



to our present study of the spin-GORM model. The non-Markovian dynam- 
ics of the subsystem is given by (|4.62|) using (|4.5<S|) and (|4.59|) . The order 
A term is neglected because in the spin-GORM model the random-matrix 
ensemble average of (e|-B|e) vanishes. The non-Markovian population dy- 
namics become 



z{t) = I de'[P++{e',t)-P_.{e' + A,t)], (5.33) 
x{t) = I de'[P+-{e',t)+P.+{e',t)] , 
yit) = i J de'[P+.ie',t)-P.+{e',t)], 



where the quantities Pss' (e, t) obey the equations 

sin(±A + huj)t 



P±±ie,t) 



A2 



duj- 



(5.34) 



P^^(e - Tuj,t)^l^ - - P^^{e,t)J^ - {e - hL,y 



and 



(5.35) 



~j da;P±^(e,t)yrr^^r^(ei(±A-M- + ei(±A+Mr) 
+ J dujP^^ie - nw,t)^ ^ - e2(e-i(±A-Mr + g-^(±A+M-) }. 

This non-Markovian spin dynamics can be computed numerically. In the 
following we will refer to this equation [ (|5.33j) with (|5.34|) and (|5.35j) ] as 
"new NM" which stands for "our new non-Markovian equation". 
Notice that the application of the non-Markovian version of the Redfield 
equation [ ()G.5|) with (|C.1|) and ()G.2j) in appendix El to the spin-GORM 
model is performed in the same way by the replacement of the Fourier tran- 
form of the environment correlation function by its expression 1)5. 29(1 and 
by neglecting the order-A term which disappears due to the random-matrix 
average. In the following we will refer to this equation as "Redfield NM" 
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which stands for " non-Markovian Redfield equation". 

We now focus on the Markovian spin dynamics. We can again simply 
apply the results of section 14.41 to the spin-GORM model. The analytical 
expression 1)4. 73() can be used if one evaluates the spin-GORM expression 
of the different environment-dependent quantities 7pop(e)5 7coh(e)) r(e) and 
Z{e,oo). This can be done easily because all these quantities depend on 
the Fourier transform of the environment correlation function which is now 
given by (|On|) . Therefore, using with (|^^ and (jUMl), the final 

solution of our new Markovian equation, to which we later will refer to as 
"new M", can be written 



z{t) = 


z{0) + 1 r 

2 r 


(e,oo) + (1 + z(e,cx)))e"^p°p(')* 


(5.36) 




z{0) - 1 


z{e - A, oo) + (1 + z{e - A, oo))e"'^p 


,p(e-A)t" 






2 






xit) « 


[x(0) cos(A/;i + r(e))t - y(0) sin(A/;i + r{e))t] e 


-7coh(<:)t 




y{t) « 


[x{0) sm{A/h + r(e))t + y(0) cos(A/;i + r(e))t] e 


~7coh(<:)t 





where the two population relaxation rates are given by 

7pop(e) = 2^-^(d(e,A/;i) +a(e + A,-A/?i)) 

h 



(5.37) 



and 



7pop(e-A) = 2Tr-^{ai€- A,A/h)+a{e,-A/h)) 

h 



while the coherence relaxation rate is 



7coh(e) = 7rj^{a{e,-A/h)+a{e,A/h)) 

2h 



i-(e-A)2 + Ji-(e + A)2 



(5.38) 



(5.39) 
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the Lamb shift affecting the coherence frequencies 



A2 - {hiof 



^^^'^t/^7:7^> (5-40) 



which after solving the integral gives 



r(.) = VA - A"''' Cretan ^ <^ + ^' + ' 

TT 



A -'^ arctan 



+ A arctan 

TT 



+ A"^^ ^arctan I ]' ~' " 1 . (5.41) 



TT 





1 

4 


(6 + A)- 


1 
2 




1 
4 


(6-A) + 


1 
2 


^(e - A)2 - 


1 

4 


(6-A)- 


1 
2 


J(6 - A)2 . 


1 
4 



The two asymptotic population values are given by 

Z(e, oo) n(e) — n(e + A) 



2;(e, oo) 



and 

z(e — A, oo) 



C(e,0) n(e) + n(e + A) 

y^+^i-(e + A)2 

Z(f - A,oo) _ n(f - A) - ;i(e) 
C(e- A,0) ~ ra(e- A) + n(e) 



\/i-(^-A)^-V^ 



e2 



i-(^-A)^ + Vi-e^ 



(5.42) 



(5.43) 



Notice that the Rcdfield Markovian equation for the spin dynamics, which 
we will refer to as "Redfield M", can be obtained in a similar way using 
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(|C.9|) . where the population as weh as half of the coherence relaxation rate 
is given by 



\2 

7 = 27T-^{a{e,-A/h)+a{e,A/h)) 



l-(e-A)2 + Jl-(e + A)^ 



(5.44) 



the asymptotic value of the populations by 

d{e,-A/h) - d{e,A/h) 



z{oo) 



d{e,-A/h) + a{e,A/h) 



(6-A)2 



(6 + Af 



(e - A)2 + 



{e + Af 



(5.45) 



and where the Lamb shift T is the same as for our new equation (|5.41|) . 
We recall that our new equation has been shown to be equivalent to the 
Redfield equation when the typical energy scale of the system (typical energy 
spacing between the system levels) can be considered small compared to the 
typical energy scale of variation of the environment correlation function 
In the spin-GORM model, this condition is reduced to n(e it A) 



n e 



5.3.2 The strong-coupling regime (A ^ 1) 

We are now interested in making a theoretical prediction for the spin dynam- 
ics when the coupling parameter A is very large and, therefore, dominant in 
front of 1 and A. 



The dynamics of the subsystem is fully described by 

0(t) = Tt p{t)a0 (5.46) 

where 9 = x,y,z. Using the following unitary transformation acting on the 
degrees of freedom of the spin 

f/ = ( J 1 (5-47) 
V V2 J 

■^The only restriction to this statement concerns the analytical expression for the Marko- 
vian coherence dynamics which has been derived by making the non-local approximation 
as explained in section Because of the fact that this assumption is not satisfied when 
A — > 0, the Markovian Redfield analytical expression has to be used in this specific case. 
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we have 

U^a^U = a^ , U^a^U = a, , U^ayll = -ay. (5.48) 
We can therefore rewrite H5.46() as 

e{t) = e^(#^-+-^s+Aa,B)tfjtp(o)f7e-^(T^-+^fl+^'^-^)*[/tCTe[/ (5.49) 
= Tr e^^(^t^-+x^s+^.B)tf/t^(o)[/e-*^(^t^-+^-^s+'^-^)*?7"f<Tef7. 

It is easy to check that 

Using the following perturbative expansion of the evolution operator to order 
zero in j 

^.aii%a.+iHs+a.B)t x-o ^-ixta^B ^ q ( \) + O ( ^) (5.51) 



and using the fact that p(0) = ps{0)p^^ , we get 
eit) Ti e'^'^^^U^ ps{0)Up''^e~'^'^^^U^agU + O + O (^j^ (5.52) 

Using the following notation 

dzBlnr]) = az\K)B\r]) = KE^\K,ri) (5.53) 

and 

HB\h)=E^\h), (5.54) 

we find that 

e{t) ~'=" Yl e^'*^'^-'''^''''(^l/5^'l^)(^|f^^/5s(0)?7|/.')(K'|t/ta,{7|K) 



K,K' ,T] 



Because is diagonal in the basis that diagonalizes Hb, we have 

{7i\pBmr])=Y,\{r]\b)\'{b\pB\b). (5.56) 
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If we perform an ensemble average over different realizations of Hb and use 
the statistics of the random-matrix eigenvectors [HI IH2| 17,'Sj , we find 



b 



b 

Equation (|5.55|) therefore becomes 

K,K',r] 

Using the relations (|5.48|) . we therefore have that 

K r} \ / \ / 



We now perform the following ensemble average over different realizations 
of 5 



Q±2iXtEr, = de n(e)e 



±2i\te 



1 



vr / 1 V 4 
= ivM^, (5.60) 

where Ji (t) is the Bessel function of the first kind and where we have used 
the fact that Ji{t) is an odd function. Performing this ensemble average 
over 1)5. 59() and using (|5.5U|) . we finally find 
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z{t) 






(5.61) 


xit) 

yit) 




x(0) 





These equations are valid up to corrections O (y) and O (^), and are av- 
eraged over the Hb and the B realizations even if we have dropped the 
explicit notation. We have therefore found a well-defined behavior of the 
system dynamics when the coupling parameter A is so large that the cou- 
pling term can be considered to contribute alone to the dynamics. We have 
again supposed that the environment is continuous (A^ oo) and that the 
random-matrix ensemble averages have been performed (x oo). 

5.3.3 Numerical results 

We will now present the numerical results on the dynamical evolution of the 
spin in the spin-GORM model. This will give us the opportunity to compare 
the theoretical predictions made in the two preceding subsections 15.3. ll and 
15.3.21 with the exact dynamics obtained numerically. 

The different parameter domains represented in figure [5?2l will play an 
important role in our discussion on the dynamics. 

Weak-coupling (domain A) 

The most important point we want to investigate is the validity of our new 
weak-coupling kinetic theory. To do this, we will now compare the non- 
Markovian version of our new equation with the numerical exact von Neu- 
mann equation, staying in the parameter domain A (see figure 15. 2|) which 
corresponds to the weak-coupling regime. Figures l5.11f a) and l5.1lT b) show 
the time evolution described by both equations at different A values using 
a A^t time scaling. The very important result is that the non-Markovian 
equation is valid not only for values of A below an upper bound, but also 
above a lower bound. The lower bound of the coupling parameter corre- 
sponds in our spectral analysis to the border between the localized and the 
Lorentzian regime of the SOE. We recall that this bound corresponds to a 
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Figure 5.11: Time evolution of the z component of the spin. In all figures, 
X = 50, e = and the von Neumann equations have been computed with 
5e = 0.05. (a), (b), and (c) Comparison between the exact von Neumann 
equation and the non-Markovian version of our new equation (" new NM" ) 
for different values of the spin-GORM model parameters, (a) and (b) cor- 
respond to A = 0.1 and (c) to A = 5 and A = 0.1. (d) Comparison for 
A = 0.01 between the exact von Neumann equation with = 500 and 
X = 50 and the Bessel strong-coupling result given by equation (|5.6H) . 
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value of the coupling parameter for which the first nonzero correction to 
the unperturbed energy levels of the total system becomes of the order of 
the mean level spacing between these coupled levels (A = yl/iV). This 
bound however disappears in the limit of a continuous environment spec- 
trum (N oo). The upper bound corresponds to the transition between 
the Lorentzian 1 regime and the Lorentzian 2 regime (see figure EHl- This 
bound is located to a value of the coupling parameter for which the width 
of the SOE in the Lorentzian regime becomes larger than the typical energy 
scale over which the density of states of the environment varies. In figure 
15. Uf a), we can see that when A is smaller than the upper bound (localized 
regime of the SOE), the non-Markovian version of our perturbative equation 
reproduces the initial dynamics but not the total relaxation process to the 
equilibrium value. One can see in figure IHTTTT b) that our perturbative equa- 
tion again looses its validity when A cannot be considered as a perturbation 
anymore. 

Non-Markovian vs. Markovian dynamics 

We have often used the Markovian approximation in the various derivations 
of our new weak-coupling master equation presented in this thesis. We 
argued that this approximation was justified on time scales longer than the 
characteristic time scale of the environment Tc (which corresponds to the 
typical decay time of the environment correlation function), i.e., if Tc is 
much smaller than the characteristic time scale of the subsystem r^. In the 
spin-GORM model, the characteristic time scale of the environment given 
by (|5.32l) is of the order r;, = 2tt and the characteristic time scale of the 
spin dynamics is given by = 27r/A. It is interesting to test this prediction 
by numerically comparing the Markovian version ("new M") with the non- 
Markovian version of our equation ("new NM"). Notice that, using the A^t 
time scaling, the Markovian equation is independent of A. The discrepancy 
between the non-Markovian and the Markovian equation is therefore a direct 
measure of the non-Markovian effects. The curves of figures I5.1ir a) and 
I5.1ir b). for which < r<j, confirm our theoretical predictions. In fact, in 
figure E^^a), the Markovian and the non-Markovian versions of our new 
equation are so close that we only plotted the second one. This is expected 
because the discrepancy is predicted to occur on a time scale 27rA^ which 
is here too short to be noticed. However when A increases on the curves of 
figure l5!T]T b). we see that the non-Markovian initial dynamics increases due 
to the increase of the predicted discrepancy region 27rA^ . 
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Pure non-Markovian dynamics (domain B) 

In the parameter domain B (see figure 1^^ . we are in tiie situation wliere 
Tfe < Tg. Here, tlie Markovian version of our equation is predicted to be 
never valid and only the non-Markovian version can be used. This situation 
is shown in figure I^TTTT c'). where we compare the Markovian and the non- 
Markovian versions of our perturbative equation with the numerical exact 
dynamics given by the von Neumann equation. As predicted, the Markovian 
version of our equation does not reproduce at all the dynamics. But the non- 
Markovian equation reproduces this behavior with a very high accuracy. 
The frequency of the oscillations corresponds to the subsystem frequency 
Ts = 0.47r and the damping of these oscillations occurs on a time scale 
corresponding to the environment characteristic time scale = 27r. Notice 
that this pure non-Markovian dynamics scales in time according to t. 

Strong-coupling (domain C) 

We now want to verify our theoretical prediction of the spin dynamics in 
the strong-coupling regime corresponding to the parameter domain C (see 
figure 15. 2j) . We see in figure I5.1ir d) that when A becomes large enough 
for the coupling term to dominate the total Hamiltonian, the dynamics 
obeys our predicted Bessel- function behavior derived in equation (|5.61|) . It 
is important to notice that this behavior scales in time according to Xt. 

Summary of the different dynamical regimes 

We summarized in figure 15.121 the different types of spin relaxation that 
one can encounter in the spin-GORM model. These different types of spin 
relaxation correspond to the different parameter domains A, B and C. Each 
of these different relaxations has a characteristic time scaling. For each of 
them we have found a valid theoretical description. 

Our new equation vs. Redfield equation 

We have motivated the extension of the Redfield theory to our new weak- 
coupling kinetic theory by the fact that when the system energy becomes 
larger than the typical energy scale of variation of the density of states 
of the environment, the energy conservation of the total system has to be 
taken into account. We now want to verify this statement by comparing the 
Redfield equation predictions (Markovian and non-Markovian) for the spin 
population relaxation in the spin-GORM model with the exact numerical 
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Equat ion^^. 


A 


B 


C 


new M 


Yes 


No 


No 


new NM 


Yes 


Yes 


No 


Bessel 


No 


No 


Yes 


time scaling 


2 

X t 


t 


xt 



Figure 5.12: Table summarizing the different types of spin relaxation that 
one can encounter in the spin-GORM model with their respective valid the- 
ory and time scaling. The parameter domains are defined in figureESl "new 
M" and "new NM" are the Markovian version and the non-Markovian ver- 
sion of our new weak-coupling kinetic equation, respectively, and "Bessel" 
refers to the strong-coupling equation (|5.6H) . 
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Figure 5.13: Comparison between the exact von Neumann equation ("von 
Neumann"), the Markovian and the non-Markovian versions of our new 
equation ("new M" and "new NM" , respectively) and the Markovian and 
the non-Markovian versions of the Redfield equation ("Redfield M" and 
"Refield NM" , respectively). Everywhere A = 0.1 and for the von Neumann 
equations = 2000, x = W and 6e = 0.05. For (a) and (b): A = 0.01, 
e = 0.25. For (c): A = 0.1, e = 0.25. For (d) and (e): A = 0.5, e = -0.25. 
For (f): A = 5, e = 0.25. 
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results and with the prediction of our new master equation (Markovian and 
non-Mar kovian) . We place ourself in the parameter domain A, where the 
weak-coupling applies, and choose an environment with a number of states 
such that the weak-coupling kinetic description is satisfied (A > \/\/N). 
Figures 15. 13r al. [5.13f cl and 15.13r d) depict the global relaxation of the spin 
population for increasing values of the energy splitting A. In accordance 
with what we argued on theoretical grounds in chapter ^ we see in these 
figures that the larger the energy splitting A of the system is, the bigger is 
the difference between the Redfield equation and our master equation. We 
also see that our equation always fits very well with the exact von Neumann 
equation, which is not the case for the Redfield equation. These results are 




t 

Figure 5.14: Population and coherence dynamics in the spin-GORM model 
for A = 0.1, e = —0.35, and A = 0.1. We compare the exact von Neumann 
equation ("von Neumann") computed with N = 2000, x = 10) s-^d e = 0.05 
with the Markovian version of our new equation ("new M") as well as with 
the Markovian version of the Redfield equation ("Redfield M"). 

a really strong support for the importance of our new theory. As expected 
and as already discussed, we see in figures I5.13f a') and I5.13( c)- I5.13f d') that 
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the Markovian equations are, on time scales longer than Tc = 27r, really 
close to the non-Markovian equations. An interesting result can be deduced 
by comparing the non-Markovian results of the Redfield equation with the 
non-Markovian results of our new equation in figures IB-LSf b) and l.^.l.Sf e)- 
IS.l^f f). In fact, on a shorter time scale of the order Tc = 2tt, there is no 
difference between the Redfield and our new equation, and both give cor- 
rect results. The discrepancy between the Redfield equation and the exact 
dynamics when A is large enough, only occurs on longer time scales. This 
can be interpreted by the fact that, in quantum mechanics, energy conser- 
vation only takes place on long time scales ^. This is particularly explicit in 
figure IKTTHT f ) where the relaxation of the spin population occurs in a pure 
non-Markovian regime (the Markovian version of the Redfield equation as 
well as our new equation predict no evolution in this case). We finish our 
comparison between the Redfield theory and our new weak-coupling kinetic 
theory by showing in figure [5TH the dynamics of the population as well as 
of the coherences of the spin. As predicted in section 15.3. 11 the Redfield co- 
herences dynamics is really close to the one predicted by our new equation. 
They only slightly differ by their Lamb shift. However, our new equation 
reproduces the exact dynamics with a higher accuracy than the Redfield 
equation. 

Energy-distributed populations and coherences 

We conclude the comparison of the exact spin dynamics with our new equa- 
tion by representing the evolution of the energy-distributed elements of the 
spin density matrix. These objects are the central elements of our new the- 
ory and are obtained, as we have seen in chapter |1J by a projection of the 
total density matrix. The very high similarity between the numerically exact 
dynamics of the energy-distributed spin density matrix element (obtained 
by the projection of the total density matrix evolving according to the exact 
von Neumann equation) and this same dynamics predicted from our new 
theory (see figure 15.151 and I5.16|) is once again a really strong support for 
the validity of our new theory. The non-local aspect (in energy) of the pop- 
ulation dynamics due to the energy conservation of the total system can be 
very well seen in figure IHTl^l The almost local aspect of the coherence dy- 
namics can be seen in figure E^Sl One notices however the slight non-local 
effect at an energy of itA from the initial energy. 



■^This fact can be easily seen on the Pauli equation. 
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Figure 5.15: Representation of the energy-distributed populations P+-|-(e, t) 

and P (e,i) which are central quantities in our new theory. The plots on 

the right have been obtained from the exact von Neumann equation and 
the ones on the left from the non-Markovian version of our new equation. 
The remarkable similarity between the right and left graphs is a strong 
support for the validity of our new theory. We here have A = 0.1025, 
e = 0.25, and A = 0.1. The exact results from the von Neumann equation 
have been obtained with N = 2000 and x = 10- The initial condition 
of the subsystem is given by (1|/55'(0)|1) = 2/3, {—l\ps{0)\ — I) = 1/3 and 
{I\ps{0)\ — 1) = {—l\ps{0)\l) = \/2/3. For a better graphical representation, 
the initial distribution of the reservoir is not exactly defined as a delta but 
a Gaussian with standard deviation 0.01. 
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Figure 5.16: Representation of the real and imaginary parts of the energy- 
distributed coherence TZe[P-\ (e,i)] and Im[P^ (f;*)] which are central 

quantities in our new theory. The plots on the right have been obtained 
from the exact von Neumann equation and the ones on the left from the non- 
Markovian version of our new equation. The remarkable similarity between 
the right and left graphs is again a strong support for the validity of our 
new theory. The conditions are here the same as in figure E^^I 



5.3 Spin dynamics 

5.3.4 Average versus individual realizations 



123 



An interesting point is the comparison, for the subsystem dynamics, be- 
tween the averaged (random-matrix ensemble averaged or microcanonicahy 
averaged) dynamics and the dynamics of an individual reahzation within 
the statistical ensembles. This latter corresponds to a dynamics generated 
by an initial condition that corresponds to a pure state and without any 
random-matrix average (x = !)■ 



We see in figure l^^TTT cl-ff) the dynamics of a system with a small energy 
spacing A = 0.1 for different values of A^A^. The solid line represents the 
random-matrix and microcanonicahy averaged dynamics. The dashed lines 
depict some of the individual members of the random-matrix ensemble cor- 
responding to an initial pure state inside the total unperturbed energy shell. 
We see that the larger X'^N is, the closer the individual trajectories are from 
the averaged trajectories. We therefore have, when X'^N is large enough, 
that the individual realizations are self-averaging in the random-matrix and 
microcanonical ensembles. In order to quantify this behavior, we plotted 
in figure IS.lTf a) the variance between the individual trajectories and the 
averaged trajectories as a function of time for different values of X'^N. We 
observe that this variance decreases as X'^N — > oo. In figure IF.iyf b). we 
show that the asymptotic value of the variance decreases with a power-law 
dependence with respect to X^N. We again see the connection with the SOE 
regimes. In the localized regime (figure l5.17r c)). the dynamics is governed by 
very different individual trajectories oscillating with a very few frequencies 
that differ from one individual trajectory to another. This is a consequence 
of the fact that the perturbed levels are still close to the non-perturbed 
ones and are only slightly affected by neighboring non-perturbed levels. In 
the Lorentzian regime (figures l5.17r d)-(e)). each individual trajectory fol- 
lows roughly the averaged trajectory and contains a very large number of 
different frequencies. This shows that the interaction "mixes" many of the 
non-perturbed levels, deleting the discrete structure of the spectrum. There- 
fore, we can say that our new master equation derived in chapter |3] holds 
with a given accuracy for a majority of individual trajectories if A is large 
enough to satisfy A > CN~^ with v < \ and a constant C > 0, in the limit 
N ^ oo. 
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Figure 5.17: In all the figures A = 0.1, = 500, e = 0, and 6e = 0.05. (a) 
Variance between individual trajectories and the averaged one (x = 100) as 
a function of time for different values of A^A^. (b) Power-law dependence 
between the equilibrium value of this variance and A^A'^. (c)-(f) Individual 
trajectories of the ensemble (dashed lines) and the ensemble averaged tra- 
jectory (solid line). In (c) A^iV = 0.1, in (d) A^A^ = 1, in (e) A^A^ = 10, and 
in (f) A^A^ = 100. 
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5.4 Equilibrium distributions 

In this section, we want to study the properties of the distribution reached 
on very long time in the spin-GORM model. In other words, we want to 
study the spin-GORM equilibrium properties. In subsection l5.4.H we study 
the equilibrium properties of the spin. In subsection l5.4.21 we investigate the 
conditions under which the spin thermalizes. Finally, in subsection 15.4.31 
we study the thermalization of the total system. 



5.4.1 Equilibrium distribution of the subsystem 

As we have seen in the preceding sections, the spin coherences always decay 
to zero on long time scales due to the environment interaction. Therefore, 
the only element left in the spin density matrix on these long time scales 
are the populations which are described by the spin observable (t^. This 
observable evolves in time according to 

z{t) = Tt p{0)e'^'a,e-'^' 

= ^{a\p{0)\a'){a'\a,\a)e'^^'^-^'^'^K (5.62) 

a,a' 

The time-averaged value of z{t) is obtained performing the following time 
average: 

1 

z{oo) = lim — / dt z{t) 

= ^(a|/5(0)|a)(a|a^|a). (5.63) 

a 

We see that this time averaged value clearly depends on the initial condition. 
An important result of our numerical simulation, which holds everywhere 
on the parameter space and for all kinds of initial conditions, is that the 
observed equilibrium value given by the exact von Neumann equation cor- 
responds very well to the time-averaged value ()5.63() . The study of the 
equilibrium properties of the subsystem is therefore equivalent to the study 
of the time averaged quantities ()5.63|) . 

We will use an initial condition for the total system which is the product of 
a microcanonical distribution at energy e for the environment with a spin 
up state for the subsystem 

^(0) = |1)(1| '-^^ = Y. (^M) 
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We have therefore that 



{a\m\c^)=Yl 



n 



6{En - e) 
n(e) 




(5.65) 



An interesting point is to understand when the Markovian version of our 
new weak-couphng kinetic equation predicts the correct population equihb- 
rium value. We want therefore to compare the exact equilibrium distribution 
given by the time averaged populations z{oo) of ()5.63|) with z{e, oo) in 1)5. 42() 
(because for the initial condition (|5.64j) . z(0) = 1 and therefore the equi- 
librium distribution predicted by our kinetic equation in 1)5. 36() is such that 
z{oo) = z{e, oo)). 

This comparison is done in figure I^TTHl where we plotted the time averaged 
value z(oo) as function of the initial energy of the environment e. The dif- 
ferent curves in figures I5.18f a') and I5.18f b') correspond to different values 
of A. We compare these curves to the A independent curve given by our 
new weak-coupling kinetic equation and corresponding to a microcanoni- 
cal distribution inside the energy shell of the total system. As expected 
from our preceding study of the dynamics, we find again that, when A is 
too small, the "mixing" between the levels inside the total energy shell is 
not sufficient to reach equipartition in this shell and our kinetic equation, 
which assumes equipartition, overestimates the equilibrium values. If A is 
too large, our kinetic equation again gives bad results as expected from the 
restriction to weak couplings of our equation. As already observed in the 
previous section, the prediction of our new kinetic equation are correct in 
a characteristic region of A. The beginning of this region corresponds to 
A values such that A > yl/iV at which the transition from the localized 
to the Lorentzian regimes occurs in the SOE. Remember that this critical 
value goes to zero when — > oo. The end of this region corresponds to A 
values that cannot be considered anymore as perturbative (when the width 
of the Lorentzian SOE becomes larger than Se the typical energy scale of 
variation of the environment density of states). Figures l5.18r c) and l5.1^ d') 
show very explicitly that z{oo) scales like A^A^. This again confirms our 
previous analysis which pointed out the importance of the quantity A^A^. 

5.4.2 Thermalization of the subsystem 

One of the important questions is to understand the conditions under which 
the subsystem thermalizes under the effect of a weak interaction with the en- 
vironment or, in other words, under which conditions the subsystem relaxes 
to a canonical distribution corresponding to the microcanonical temperature 
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Figure 5.18: Comparison between between the equilibrium value of 2(00) 
given by the Pauh equation which corresponds to a micro canonical distri- 
bution on the total system and the exact values given by the time-averaged 
value for different values of A. These equilibrium values are depicted as a 
function of the initial microcanonical energy e of the environment. 5e = 0.05 
in all figures. The parameter values are (a) A = 0.01 and = 500; (b) 
A = 0.5 and N = 500; (c) A = 0.01 and iV = 200 - 2000; (d) A = 0.5 and 
iV = 200 - 2000. 
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of the environment. 



We begin by recalling these conditions in the general case of a small sub- 
system weakly interacting with its environment. The isolated total system 
is composed of the subsystem and the environment and has the total energy 



tot 



e + e. 



(5.66) 



e is the energy of the environment and e the energy of the subsystem. We 
suppose that the contribution of the interaction energy between the envi- 
ronment and the subsystem is negligible compared to the total energy. The 
microcanonical environment entropy is defined as 



SBie) = klnnBie), 



(5.67) 



where ri_B(e) is the number of states of the environment available at energy e. 
This number can be related to the density of states of the environment nB(e) 
using the fact that ^b{^) = n-B{e)de, where 6e is a small energy interval but 
contains many states of the environment. The microcanonical temperature 
of the environment is given by 



1 



^ dSBje) 

TB{e) de ■ 

It can be expanded in the subsystem energy as 



Tsie = E, 



tot 



Tsie = E, 



tot) 



dTeie = Ei 



tot) 



de 



+ 



(5.68) 



(5.69) 



because we suppose that the subsystem energy is much smaller than the 
environment energy. The specific heat capacity of the environment is 



1 



dTBie) 



Cvb{^) de 

If the condition 

\CvB{e)\ > 

is satisfied, the temperature expansion can be truncated as follows 
TB{e = Etot - e) = TB{e = E^ot)- 



TB{e) 



(5.70) 



(5.71) 



(5.72) 



Therefore, we understand that if equation (|5.71|) is satisfied, the environment 
plays the role of a heat bath because its temperature is almost not affected 
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by the subsystem energy. 

Now, we suppose further that the interaction between the subsystem and the 
environment, even if small, is able to make the total probability distribution 
microcanonical on the total energy shell at energy i^tot- We suppose also 
that the energy levels are discrete and, therefore, that ii^tot = Es + 
where s and b are discrete indices, respectively, for the subsystem and the 
environment. Therefore, the probability Ps{Es) for the subsystem being at 
energy Es is given by 

Ps[Es) = o ^ ' ^^-^^^ 

where rtB{Eh) is the number of states of the environment available at energy 
£";,, and Vt,tot{Etot) the number of states of the total system available at 
energy Etot- Using equations (|5.68|) and (|5.71j) . one gets that 

^tot(-E'tot) ^tot(-E'tot) 

Using the normalization Ps{Es) = 1, one finally gets the well-known canon- 
ical probability distribution for the subsystem 

e kTg(EB=Et^t) 

PsiEs) = , (5.75) 



where Z = Y^s ^ '''^B(EB=Etot) . 

We conclude that two conditions are necessary in order to thermalize the 
subsystem to a canonical probability distribution due to the interaction with 
the environment: a large heat capacity of the environment |C^,B(e)| ^ Ir^ejl 
and a microcanonical distribution on the total energy shell. 

Let us apply this result to the spin-GORM model. 
The density of states of the environment is in this case 



Therefore the microcanonical temperature of the environment is given by 

Tsie) = ^ (5.77) 
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and the heat capacity is 



ke 



e2 + i' 



(5.78) 



The first condition (|5.71j) in order to thermalize the spin to a canonical 
probabihty distribution becomes 



\C,Bi€)TBie) 



e2 + i 



> A 



(5.79) 



and is depicted in Fig l5.19l The second condition to have a micro canonical 
distribution on the total energy shell is satisfied (as we discussed in section 
I5.4.1|) when A^A'^ > 1. In this case, the populations of the spin obey the 
equilibrium distribution predicted by our new kinetic equation (see (|5.36() ) 



z[oo) = z[e, ooj h z[e 



z(0)-l 



(5.80) 



5.4 Equilibrium distributions 



131 



where 



z(e, oo) 



Z{e, oo) 



Vi-(^)^-Vi-(^+A)^ 



(5.81) 



The canonical distribution of the spin populations, at the micro canonical 
temperature of the environment, is given by 



If the two conditions for the spin thermalization are satisfied, then z(TB(e)) 
and 2;(oo) should be equal. 

The comparison between these two quantities can be seen in figure l5^2fll for 
different subsystem energies. We see, as expected, that the smaller the spin 
energy is the better the comparison is. 

We can therefore conclude that under these two conditions (| ^-ij^y^ I > A 
and A^A^ > 1), the random matrices of the spin-GORM model can model an 
environment that behaves as a heat bath and that thermalizes the spin. The 
particularity of this environment due to the semicircular density of states 
is that the microcanonical temperature can be negative (when the initial 
energy of the environment is positive). Usual environments have positive 
microcanonical temperatures due to their monotonic increasing density of 
states. 

5.4.3 Thermalization of the total system 

Until now, we have always considered initial conditions of the environ- 
ment corresponding to microcanonical distributions at a given energy like 
in 1)5. 64p . We now want to consider initial conditions of the environment 
corresponding to a canonical distribution at a given temperature. We also 
suppose that the spin is up, so that the total initial condition reads 




(5.82) 




(5.83) 



We therefore have that 




(5.84) 




Figure 5.20: Comparison between Eqs. (|5.82() and ()5.8H) for different values 
of A. The parameter values are (a) A = 0.01; (b) A = 0.05; (c) A = 0.1. 
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It is interesting to ask how the probability distribution looks hke at equihb- 
rium when the interaction is on. 

One can clarify this point by plotting {a\p{0)\a) versus energy. One uses 
the following energy representation 

P{e) = Y,6{Ea-e){a\m\a)- (5-85) 

a 

If the total system thermalizes and reaches a canonical distribution for the 
total system at an effective temperature P~g , one would have that 

P{e) = (5.86) 

^tot 

because 

(a|p(0)|a) = — . (5.87) 

^tot 

As we shall see, it is the case if, again, A is large enough to induce "mixing" 
between the states: A'^A^ > 1. 

Indeed, one sees in figure I5.21f a) that, for A = 0, the states of the total 
system corresponding to the upper level of the spin are exponentially pop- 
ulated and the ones corresponding to the lower level are not. When the 
interaction is turned on and increased, one can notice that the probability 
distribution starts to accumulate around a mean effective canonical distri- 
bution. As expected, this accumulation becomes significant when A^A^ > 1 
and, in this case, the total system can be considered as having thermalized. 
One can calculate the final effective temperature that the total system has 
reached after interaction. This effective temperature is depicted in figure 
I5.21f b) as a function of A. The correlation coefficient indicates whether the 
exponential fit of the final effective temperature is good or not. The effective 
temperature obeys the following law: P^s = t+^- We also show in figure 
I5.2ir c') the comparison between the time-averaged value of z{t) ( (|5.63|) with 
H5.84() ) and its canonical average computed with the effective temperature 
f (|5.63p with (|5.84|) V One sees that, when A^A^ > 1, both coincide. 
It is important to realize that this thermalization is not statistical in the 
sense of the equivalence between the ensembles. It is an intrinsic thermal- 
ization due to the complexity of the interaction between the states. This 
thermalization appears at a critical value of the coupling parameter when 
the interaction term becomes of the order of or larger than the mean level 
spacing of the total system. 
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Figure 5.21: In all the figures A = 0.01, N = 500, and /3 = 2. (a) Probabil- 
ity at equilibrium P(e) = Dist(e) of being at the total system energy e when 

-l3Hg 



starting from the initial condition p{0) = |1)(1| 



Zb 



with ^ 



0.5. (b) Ef- 



fective temperature of the equilibrium probability distribution obtained by 
fitting the data of (a) with a canonical distribution, (c) Comparison be- 
tween the exact equilibrium population value f H5.63() with H5.84() ') and the 
canonical one f (|5.63|) with (|5.87)) ) computed with the effective temperature 



Chapter 6 

Quantum diffusion 



In this final chapter, we will focus on a particularly important quantum ki- 
netic process: quantum diffusion. 

Quantum diffusion is a quantum transport process where the transported 
quantity is the position of a particle on a spatially extended subsystem. Un- 
derstanding quantum diffusion is of huge importance, especially in nanosys- 
tems. When the particles transported on the nanosystem have a charge, 
the transport properties of these particles determine the electric conduction 
properties of the nanosystem. And one of the most common ways to study 
nanosystems such as quantum dots, nanotubes, or single linear molecules, 
is through their electric conduction. Furthermore, the technological appli- 
cations of these electronic devices are really promising. 
If the subsystem is invariant under spatial translations, the Bloch theorem 
allows us to associate a Bloch number or wavenumber with each mode of the 
evolution superoperator. The dynamics of such a translationally invariant 
system is diffusive on long time scales if the dispersion relation, connecting 
the smaller real part of the mode eigenvalues (corresponding to the slower 
exponential relaxation modes) to the Bloch number, is parabolic around 
zero wavenumber. The proportionality coefficient between the eigenvalue 
and the squared wavenumber is then the diffusion coefficient. Such a dis- 
persion relation cannot be found in an isolated quantum subsystem where 
the real parts of the mode eigenvalues are always zero. Therefore, we want 
to see if such a diffusive dispersion relation can occur in a translationally 
invariant subsystem weakly interacting with an environment. Our study 
is based on the Redfield equation (see section 13.5(1 and the environment is 
therefore considered infinite. 

Various attempts have been made in the literature to understand diffusion in 
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quantum subsystems. Many of the early studies focused on polaron motion 
123 linj • Later studies used a kind of Quantum Brownian Motion model ex- 
tended to periodic potentials [nHHHSHliEnilHIllHllinil- However, it seems 
to us that none of them have avoided the presence of a non-translationally 
invariant term in the interaction Hamiltonian between the subsystems and 
its environment. Furthermore, they almost always use a path integral for- 
malism, which has the advantage of sometimes giving non-perturbative re- 
sults but is not easy to use, especially for numerical simulations. Our study, 
which relates to the full diagonalization of the evolution operator and there- 
fore gives access to the complete dynamical information, seems to us to be 
a particularly interesting approach to diffusion in quantum systems. 

The plan of this chapter is the following. Our translationally invari- 
ant model is defined in section 16.11 by weakly coupling a one-dimensional 
tight-binding Hamiltonian with a single energy band to a delta-correlated 
environment. The dynamics of this system is ruled by a Redfield quantum 
master equation. The long-time evolution of the model can be studied in 
terms of the eigenvalues and the associated eigenstates of the Redfield su- 
peroperator, as explained in section The eigenvalue problem is exactly 
solved for a finite chain in section 1^31 We show in section that diffusion 
dominates the long-time dynamics if the chain is large enough. The Liouvil- 
lian spectrum of the infinite chain as well as the temperature dependence of 
the diffusion coefficient are given in section [6.51 We finally summarize the 
results in section HTHl 



6.1 Defining the system 
6.1.1 Subsystem 

Let us consider a translationally invariant one-dimensional subsystem de- 
scribed by the following Hamiltonian 
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represented in the site basis |/), where / takes the values / = 0, 1, . . . , A^— 1. N 
is the length of the chain. We have here chosen periodic (Born- van Karman) 
boundary conditions. 

Such a Hamiltonian, called tight-binding or Hiickel Hamiltonian, is a sim- 
ple model of spatially periodic subsystem and has therefore a wide range of 
applications. The tight-binding Hamiltonian is for instance commonly used 
in solid state physics to describe the electronic band structure of weakly 
interacting solids in polymer physics to describe electronic conduction 
along the polymer backbone or polymer optical properties |^|^, as well as 
to understanding the conduction properties of carbon nanotubes jHO]. This 
Hamiltonian describes a process of quantum tunneling from site to site. 
The parameter A is given in terms of the overlap between the wavefunctions 
localized at the sites and is proportional to the quantum amplitude of tun- 
neling. 

The stationary Schrodinger equation of the tight-binding Hamiltonian is 
given by 

Hs\k)=ek\k), (6.2) 

where the eigenvalues are 

ek = Eo-2Acosk— (6.3) 



= ^e^'^t, (6.4) 



and the eigenvectors 



with k = 0,1, N — 1. The closure relation is given by 

N-l 

^ ^ e^K^-nt = 4,. (6.5) 

1=0 

Accordingly, the energy spectrum of the Hamiltonian (|6.1|) contains a single 
energy band of width 2A and the motion of the particle would be purely 
ballistic without coupling to a fluctuating environment. 

6.1.2 Coupling to the environment 

We now suppose that the subsystem is embedded in a large environment. 
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The Hamiltonian of the total system composed of the one-dimensional 
chain and its environment is given by 

^tot = Hs + Hb + XY, SnB^ , (6.6) 

where Hb is the environment Hamiltonian, the subsystem coupling oper- 
ators, the environment coupling operators, and A the coupling parameter 
which measures the intensity of the interaction between the subsystem and 
its environment. 



The dynamics of the total system is described by the von Neumann 
equation 

= LtPtotit) = -l[Ftot, Ptot(i)], (6.7) 

where £tot is the Liouvillian superoperator of the total system. The reduced 
dynamics for the density matrix p{t) = TrB/5tot(i) of the subsystem is known 
to obey a Redfield quantum master equation for weak coupling to the en- 
vironment. We have systematically derived this equation by second order 
perturbation theory in section 1531 from the complete von Neumann equation 
for the total system (|6.7|) . On time scales longer than the correlation time of 
the environment, the Redfield quantum master equation is Markovian and 
reads 

dp j, 

Tt - ^^^^^ 

= -'^[Hs.p\ + ^,Y.^f^PS^ 

+Sipfl - S^f^p - pf^Si) + 0(A3), (6.8) 
where C^ed is called the Redfield superoperator and where 

r« = V/ dTa,,,(^)e"^^^"S,'e^^^". (6.9) 

The correlation function of the environment which contains all the necessary 
information to describe the coupling of the subsystem to its environment is 
given by 

a,,.(r) = Trsp^'et^^^^^e-*^^^^,. (6.10) 
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where is the canonical equilibrium state of the environment. Rigorous 
results show that the long-time dynamics of the subsystem is indeed ruled 
by the Redfield superoperator according to 

^ = TTBLtPtotit) tnedPit) + 0(A3) (6.11) 
for weak coupling to its environment |441 145| . 

Let us now specify the interaction of the subsystem with its environment. 
We define the subsystem coupling operators as 

{l\Sn\l') = 5.1 5w . (6.12) 

They are diagonal operators in the site basis of the subsystem taking the 
unit value if the particle is located on the site k and zero otherwise. These 
operators have the properties: 

= 5, , (6.13) 

for n = 2, 3, ... and 

" " " . (6.14) 



We need now to specify the environment operators by the choice of the 
environment correlation functions. We make two assumptions: 
Assumption 1: The environment dynamics has a very fast evolution com- 
pared to the subsystem dynamics. Therefore, the environment correlation 
functions decay in time so fast with respect to the characteristic time scales 
of the subsystem evolution, that they can be assumed to be Dirac delta 
distributions in time. 

Assumption 2: The environment has very short range spatial correlations, 
much shorter than the distance between two adjacent sites of the subsys- 
tem. Accordingly, the environment correlation functions decay in space so 
fast that they can be assumed to be given by Kronecker delta in space. 
These two assumptions mean that each of the environmental fluctuations at 
the different subsystem sites are statistically independent. The environment 
correlation functions are thus given by 

a^^.{T) = 2Q5{T)5^^. , (6.15) 

and the operators (|6.9p in the Redfield equation 1)6. 8|) therefore become 

f^ = QS^. (6.16) 
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Because of the fast decay of the temporal and the spatial correlations 
(|H.15|) and the properties of the subsystem coupling operators (|H.12|) . the 
Redfield equation we have to solve takes the form 

^ = ^RcdP (6.17) 
at 

= -'-[Hs,p\ + ^QY.{2S^pS^ - SIp - pSl) + 0(A3) . 

K 

It can easily be verified by projecting this equation onto the site basis that 
it is translationally invariant (shifting all the site indices appearing in the 
projected equation by a constant does not modify the equation). Further- 
more, this equation preserves the complete positivity of the density matrix 
because it has the Lindblad form [SH] which is the result of a coupling with 
delta correlation functions j3Uj . 

6.2 Diagonalizing the Redfield superoperator 

The eigenvalues Sy and associated eigenstates p'^ of the Redfield superoper- 
ator are defined by 

-^Rcd P^ = Su p", (6.18) 

where v is a set of parameters labeling the eigestates. Because the Red- 
field superoperator is not anti-Hermitian, its eigenvalues can be complex 
numbers with a nonzero real part. The eigenvalue problem of the Red- 
field superoperator is important because the time evolution of the quantum 
master equation can then be decomposed onto the basis of the eigenstates 
as 

p{t) = e^^^^'piO) = Cu{0) e'^'p". (6.19) 

The dynamics is therefore given by a linear superposition of exponential or 
oscillatory exponential functions. Since the reduced density matrix of the 
subsystem has A''^ elements, there is a total of A^^ eigenvalues and associated 
eigenstates. 

6.2.1 Bloch theorem 

Since the system is invariant under spatial translations, we can apply the 
Bloch theorem to the eigenstates of the Redfied superoperator. Thanks to 
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this theorem, the state space of the superoperator can be decomposed into 
independent superoperators acting onto decoupled sectors associated with a 
Bloch number, also called wavenumber 

Let us define the superoperator Ta of the spatial translation by a sites 
along the system (a is an integer) as 

i'^a)ll',mm'Pmm' = Pl+a,l'+a J (6.20) 

m,m' 

where we use the notation 

{l\p\l') = pw . (6.21) 
This superoperator must have the group property 

= = %.+a' ■ (6.22) 

Because of the translational symmetry of the system, the translation super- 
operators commute with the Redfield superoperator 

faC-R.cdP{t) = CKcdTaPii) ■ (6-23) 

Therefore, the Redfield superoperator as well as the translation superoper- 
ators have a basis of common eigenstates. If T(a) denotes the eigenvalues of 
the translation superoperator, we have that 

t = r(a) p-" , (6.24) 

where, because of the unitarity of Tr, 

|r(a)p = 1 . (6.25) 



Equation (|6.22)) implies 

T{a + a') = T{a) T{a') . (6.26) 

Because of equations (|6.25p and (|6.26() , we find that 

r(a) = e^"?*^ , (6.27) 

where q is the Bloch number, also called the wavenumber. Finally, we find 
that 

pr+ar+a = e^'^Pw ■ (6.28) 
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A useful consequence is that 

Pu' = e'''Po,i'-i ■ (6.29) 

In order to determine the allowed values of the Bloch number, we write by 
using equation (|6.4|) that 

pu' = j^Eif^inky^''-'''^'^ (6.30) 

k,k' 

and 

= (6-31) 

k,k' 

Due to equation (|(i.28|) . we also have 

pr+i,«'+i = e'^pf,. (6.32) 

Multiplying both sides of equations (IOT|l and ((02l) by {l'\k"'){k"\l), taking 
the sum J2i i' of it, and identifying them, we get 

e'liklp-'lk') = e'^''-''">^{k\p''\k'). (6.33) 

We can now notice that if q {k — k')^, then {k\p^\k') = 0. Using finally 
the periodicity 

P/W = Pu', (6-34) 
Ph'+N = Pw > (6-35) 

we find with 1)6. 28() that the Bloch number takes the values 

27r 

q = jj^, where j = 0, 1, . . . , - 1 . (6.36) 

Accordingly, the Redfield superoperator can be block-diagonalized into 
independent blocks, which each contains N eigenvalues as we shall see in 
the following. 

6.2.2 Simplifying the problem 

We will now formulate the eigenvalue problem of the Redfield superoperator 
in each sector labeled by a given wavenumber q. 
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For this purpose, we project equation H6.18() onto the site basis using the 
exphcit expression (|6.17|) of the Redfield superoperator. We get 

A 

^vPlV — \~Pl-\,V ~ Pl+l,l' + Pl,l'-1 + Pl,l'+l) 

+ ^ {Sw - I) P^i' ■ (6.37) 
Using equation (|6.29|) and replacing — Z by we get 
2QA2\ ^ 2A f . q 



-r 

Making the change of variable 



+^Poo ^01 . (6.38) 



Poi = , (6.39) 

we obtain the simpler eigenvalue equation 

Pufi = Soifi - i(3{fi-i + fi+i) , (6.40) 



where 



and 



6.3 Finite chain 

6.3.1 The eigenvalue problem 

We can write the expression ()6.4U|1 in matrix form (we no longer write the 
index v to simplify the notation) 

/// = Wf, (6.43) 
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where n are the eigenvalues, f = {fo, ■ ■ ■ , /n-i) the eigenvectors of size A^, 
and W the fohowing N x N matrix 



/ 1 

-iP -if3 

-i(3 -i(3 



iP 

-i/3 
-iP 



We look for eigenstates of the form 

Solving equation ()6.43|) with (|6.45|) gives 

• for < Z < Af- 1: 

/i = —2i P cos I 

• for / = 0: 



-iP 




(6.44) 



(6.45) 



(6.46) 



(1 - f^){A + B)- iPiAe'" + Be- 
-i/3r^e^^f (Ae^^(^-i) +i?e-*^(^-i)) = 



for 1 = N-l: 



^e*e^ + Be 



-idN _ j^N^-iN 



■{A + B) = . 



(6.47) 



(6.48) 



Solving the homogeneous linear system of equations (|6.47l) and (|6.48|) and 
replacing n by 1)6. 46(1 . one gets the characteristic equation 



with 



2iPsme [R{q) - cos ON] = sin ON, 



Using equation 1)6. 36() . we find: 



for N odd : R{qj) = , 

for iV = 4/ : R{qj) = {-ly , 

for N = AI + 2: R{qj) = -{-ly 



(6.49) 
(6.50) 



(6.51) 
(6.52) 
(6.53) 
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with / integer. From now on, we shall speak of even (respectively odd) q, 
if q corresponds to an even (respectively odd) integer j in equation 
Therefore, the characteristic equation H6.49() becomes 

• for N odd: 

2il3sme = -tan9N ■ (6.54) 

• for N = AI and q even or for N = + 2 and q odd: 
either 

ON 

2i/3sin0 = cotan-^ , (6.55) 

or 

cos 617V = 1; (6.56) 

• for N = 4:1 and q odd or for N = 41 + 2 and q even: 
either 

6N 

2if3smB = -tan—, (6.57) 

or 

cosON = -1. (6.58) 
We solve equations (|6.54j) - (|6.58|) as follows. 
6.3.2 The diffusive eigenvalue fi^^^ 

We first look for an eigenvalue /x which is real and should correspond to a 
monotonic exponential decay. With this goal, we suppose that the angle 6 
is complex 

e = (,+irj , (6.59) 

so that the eigenvalue 1)6. 46(1 becomes 

/i = —2i(3 cos 9 = —2(3 sin sinh rj — 2il3 cos cosh ij . (6.60) 

This eigenvalue is real under the condition that cos ^ = which is satisfied 
for 
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in which case 

// = -2 /? sinhr/ . (6.62) 

We notice that the condition cos ^ = is also satisfied for ,^ = ^ but it can 
be shown that this other case leads to the same eigenvalue as (|6.61() . If we 
introduce the conditions (|6.59jl and (|(161|1 in equations (|6.54|1 . (|6.55jl . and 
(|^37|) . we get 

• for N odd: 

2/3 cosh r/ = coth A^?7 ; (6.63) 

• for N = 41 and q even or for = 4/ + 2 and q odd: 



Nr] 

2/3 cosh 77 = coth ; (6.64) 



for N = 41 and q odd or for = 41 + 2 and q even: 



2/?coshr/ = tanh-^. (6.65) 

In the limit N ^ oo, the right-hand side of these equations tends to 
unity if a non- vanishing solution rj ^ exists. In this limit, this solution is 
thus given by 

rjo = arccosh — , (6.66) 



which exists only if /? < i. Because 



1 Y 



sinh7?o = y j (6-6'') 
the corresponding eigenvalue should be 

fio = \/l - (2/3)2 . (6.68) 

However, for a finite chain with N < oo, we expect a correction 6r] to the 
solution r] = rjQ + 5rj. Replacing this correction in equations (|6.63j) . (|6.64|) . 
and H6.65() . we obtain by Taylor expansion that 
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for N odd: 



— 2Afarccosh ^ 

e 

5?7 ~ 2 — , ; (6.69) 



for N = 4:1 and q even or for = 4/ + 2 and q odd: 



— A^arccosh ^ 
g 2/3 

<577 ~ 2 , ; (6.70) 

x/1 - (2/3)2 

• for N = AI and (7 odd or for = 4/ + 2 and q even: 

— Afarccosh i 

e 2p 
(5r/ ~ -2 ^ ; (6.71) 

VI - 

up to corrections of 0{5rf'). Using the expression (|6.62j) . we finally obtain 
the eigenvalue 

^S^^ = 2/3 sinh 770 + 2/3 cosh % S-q + 0{5rf) 

= Vl - (2/3)2 + 5r/ + 0(5r/2) , (6.72) 

where drj is respectively given by equations 1)6. 69() . (|6.70|) . and 1)6. 71(1 . Ac- 
cordingly, the correction 5r] to the eigenvalue decreases exponentially fast 
with the size N of the chain. 

Using equations 1)6. 41() and (|6.42l) . we finally obtain the eigenvalue 



s 



(1) 



2 



2QA 

\2 



2QX' (2Ah . qy 2QA2 

= V^^^"2j -n^ + 0{6v). (6.73) 

The eigenvalue s^^^ = corresponding to a vanishing wavenumber q = 
is always in the spectrum of the Redfield superoperator. The associated 
eigenstate describes the stationary equilibrium state. At low wavenumbers 
(;,/?—> 0, we recover the dispersion relation of diffusion 

s(i) = -Dq^ + 0{q^) , (6.74) 

with the diffusion coefficient 

D = ^, (6.75) 
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which justifies cahing the diffusive eigenvalue. 

We notice that the diffusive eigenvalue no longer exists beyond the crit- 
ical value /?c = 5- Since the matrix ()6.44|) has a total of N eigenvalues, we 
expect further non-diffusive eigenvalues in a number of — 1 for /? < ^ and 

for /? > i, as confirmed in the following subsections. 

6.3.3 The eigenvalues ji^'^^ 

Beside the diffusive eigenvalue, we expect eigenvalues corresponding to the 
chain-like structure of the matrix 1)6. 44() . To obtain these eigenvalues, we set 
tanx = 2if3sin9 so that equations (|6.54|) . (|6.55jl . and (|6.57j) can be written 
respectively 

sin{eN + x) = , (6.76) 
cos + = , (6.77) 

sin ( — + x]=0. (6.78) 

Therefore, we obtain 

i{m:-eN) = arctanh(-2/3sin6') , (6.79) 

i (^mr + I - = arctanh(-2/?sin6') , (6.80) 

i ^nvr = arctanh(— 2/3sin0) . (6.81) 

We now expand arctanh(— 2/3sin0) around /3 = 0: 

arctanh(-2/? sin 6) -2/3 sin 9--P^ sin^ 6 + 0{p^) . (6.82) 

3 

If /? = 0, the solutions of equations 1)6. 79|) . (|6.8U|) . and 1)6.81(1 are respec- 
tively given by 

TITT 

0Q = where n = 1,2, ...,Af-l, (6.83) 

(2n + l)7r , N , , 

9o = ^ where n = 0, 1, 1, (6.84) 

2mr N 
Oo = —, where n = 1,2,..., y-1, (6.85) 



6.3 Finite chain 



149 



Notice that n = is rejected in equations ()6.83|) and ()6.85() . It is due to the 
fact that 6 = does not correspond to an eigenvector because it can be seen 
that // = ^ + -B 7^ in equation (|6.45|) cannot be an eigenvector of equation 



Using the expansion H6.82() in equations ()6.79|) . ()6.8U() and (|6.81l) with 
6 = 00 + 60, we find 

M 3M 

where M = for equation (|6.79j) and M = ^ foi^ equations (|6.80j) and 
H6.81() . The eigenvalue (|6.46l) is now given by the expansion 

H = -2iPcos{eo + 59) 

° -2i/3 cos 00 + 2i(3 sin Oq 69 

cos 9o 66^ + {66^). (6.87) 

Using equation (|6.86j) in (|6.87j) gives 

^(2) -2z/3cos0o + o(z^ 

' sin2 00 + ^ sin^ 9^ + (^] . (6.88) 



M 3M \M 

Consequently, we have 

• for N odd, using equation (|6.88l) with ()6.83|) : 

(2) /3-*0 nvr 4^2 nvr 16/?^ . 4 nvr 
' = —zipcos sm sm — , 

where n = 1,2, ... ,Af - 1 ; (6.89) 

• for N = 4:1 and q even or for iV = 4/ + 2 and q odd, using equation 
(OHl) with (IFTHlj) : 

-2z/3cos^^^^ 



8/32 2(2n + l)vr 32/5"^ 4(2n + l)7r 
■ sm — 1 — sm 



N N 3N N ' 

N 

where n = 0, 1, . . . , 1 ; (6.90) 
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• for N = AI and q odd or for N = 41 + 2 and q even, using equation 
(jHl^ with (IHl^ : 

(2) 2n7r 8/32 . 2 2n7r 32/34 . 4 2n7r 

u'- = —2in cos sm sm , 

^ ^ N N N 3N N ' 

N 

where n = 1, 2, . . . , y - 1 . (6.91) 

6.3.4 The eigenvalues /i^^^ 

The solutions of equation (|6.56|) are simply given by 

2mT N 
= —, where n = 1, 2, . . . , y - 1 . (6.92) 

We reject n = and n = ^ because the corresponding eigenvector does not 



exist in these cases. Similarly, the solutions of equation (|6.58() are given by 
^^(2n^l)7r^ where n = 0, 1, . . . , y - 1 . (6.93) 
Consequently, we have the further eigenvalues: 



• for N = 41 and q even or for N = 41 + 2 and q odd, using equation 
I^^J^ with (jn22I): 

/u(^) = -2ip cos (^^^ , where n = 1, 2, . . . , y - 1 ;(6.94) 

• for N = 41 and q odd or for = 4/ + 2 and q even, using equation 
(lOHl) with (IIM : 

(-^^ ^ /(2n + l)7r 

= -2i/3cos(^ ^ ^ ^ 

iV 

where n = 0, 1, . . . , — 1 . (6.95) 
6.3.5 The eigenvalues jj^'^^ 

An important observation is that, for /3 < |, the expansion (|6.82|) which 
implies 2(3 sin 9 < 1 is satisfied everywhere, i.e., for all the values of 9 and 
therefore for all the eigenvalues. However, for (3 > ^, the Taylor expansion 
around /3 = in equation (|6.82|) is only valid if 2(3 sin 9 < 1. Therefore, a 
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transition zone exists around sin 6 ~ 1/(2/3). According to equation (|6.46|1 . 
this transition corresponds to the critical value of the eigenvalue given by 

/ie = ±iViW - 1 • (6-96) 

Therefore, for /3 > i, the expansion (|6.82() around /3 = is only valid 
if > |/ic|- For \fj.\ < \fic\, we should instead consider the asymptotic 
expansion of arctanh(— 2/3sin0) around /? = cxd: 

arctanh(— 2/3sin6) ^ -?° ^_ _ __ — _________ 



2 2/3 sin 24/33 sin^ ^ 

1 +o^' 



160/35 sin^e 



which leads to another family of eigenvalues existing for /3 > i. 



(6.97) 



If /3 ^ oo, the solutions of equations ()6.79|) . ()6.8U|) and 1)6.81(1 are respec- 
tively given by 

(n + ^)7r 

00 = , where n = 0, 1, . . . , iV - 1 , (6.98) 

2mT N 

Oo = -j^, where n = l,2,...,y-l, (6.99) 

(2n + l)7r , N , , 

00 = , where n = 0, 1, 1. (6.100) 

Because of the condition < |/Xc| with the critical values (|6.96() . we should 
only consider the angles in the interval 0o,c < < ^ ~ ^o,c with 

Oo,c = (6.101) 



so that the integer n in equations (|6.98|) . (|6.99() . and (|6. 100(1 is restricted 
to take the inte 
extreme values. 



to take the intermediate values n^[^ < n < rimax which do not reach the 



Using the expansion ((6.97(1 in equations ((6.79(1 . ((6.80(1 . and ((6.81(1 with 
= 00 + 50, we find 

/3->oo i i i 



60 



2(5M sin 0q 24/33M sin^ 6^0 160/35M sin^ 0q 



+0(^)+o(-J-,] , (6.102) 
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where M = for equation 1)6 .791) and M = ^ for equations H6.80|) and 
Using equation (|6.1fl2j) in (|6.87j) gives 

(4) 13^00 „ . ^ „ 1 / , 1 1 



_2il3coseo + — {l + ——^ + 



+0{-sLf)+o(^] . (6.103) 



M V 12/32 sin^ ^0 sin^ 

Consequently, we have 

• for N odd, using equation (|6.1U3jl with (|6.98jl : 

^(4) /3^oo I 

^ " ' N ^1o«2...2(n+|K 



l2/32 sin^ 



AT 



80/34sin4 



where n = 0, 1, ... ,iV - 1 ; (6.104) 

for N = AI and q even or for = 4/ + 2 and q odd, using equation 
(IfTM with dlTMl) : 

(4) /3^oo 2n7r 



-2i/3cos. 



N 



2/1 1 



\^ 12/32 sin2 ^ 8O/34 sin^ ^ j ' 

where n = l,2,...,y -1 ; (6.105) 

for N = AI and g odd or for iV = 4/ + 2 and q even, using equation 
with IMI): 

/.W ^=°° -2z/3cos^^^^ 

2 / 1 1 \ 

+ TT 1 + TTT^-^ + 



N \ 12/32 sin2 8O/34 sin^ / ' 

N 

where n = 0, 1, . . . , ^ — 1 ; (6.106) 
with the aforementioned restriction on the values of the integer n. 
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6.3.6 The eigenvalues /i^^^ 

For N = AI and q even or = 4/ + 2 and q odd, two special eigenvalues 
exist in the limit [3 ^ oo around 6q = Q and 6q = vr. They can be obtained 
by taking 9 = 9q + 69 and directly solving equation ()6.55|) to get in both 

cases 

Inserting in equation (|6.87j) . we obtain 

(5) - — 1 " 1 



^2z/3 + - + 0(^-) . (6.108) 



6.3.7 Description of the spectrum 



Rewriting the eigenvalues (|6.46)) of the Redfield superoperator with their 
explicit dependence in terms of equation 1)6. 42() , we get 

/"i/ = l^qe = —2if3cos9 = — 2i [^sin^") cos 6* , (6.109) 

where u can take N"^ different values because q and 9 take values each. 
Remembering that according to equation (|6.41j) 



Sv = SqB= ^2 (/igg-l), (6.110) 

we conclude that we have found in this section all the eigenvalues of the 
Redfield superoperator. We list them in Tables IHI and EUl according to the 
parameter regime in which they hold. 



We now discuss the main features of the spectrum when the different 
physical parameters are varied. This discussion is based on our analyti- 
cal results for the eigenvalues and on the comparison between these results 
and the eigenvalues obtained by numerical diagonalization of the Redfield 
superoperator. Since the eigenvalues fiqg are related to the Redfield superop- 
erator eigenvalues Sgo by equation H6.110|) , we notice that all the eigenvalues 
of the complete spectrum always satisfy < TZe < 1 or, equivalently, 
~2Q}? jf? < TZe SqQ < 0. The imaginary part of Sqo is simply proportional 
by a factor 2Q}? /f}^ to the imaginary part of ^qg. 
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Table 6.1: For wavenumbers q corresponding to /? < 2^ List of the 
eigenvalues of the matrix (|6.44j) given by equation (|6.1fl9j) . These eigen- 
values are directly related to the Redfield superoperator eigenvalues by equa- 
tion (irrrm. 



N odd 


N = 41 and q even 


N = 41 and q odd 




or Af = 4/ + 2 and q odd 


or = 4/ + 2 and q even 



/.(I) = yi - (2/3)2 ^(1) = y/i - (2/3)2 ^(1) = ^1- (2/3)2 

^ g-2iVarccoEih l/(2g) ^ ^ - Warccoeh 1/(2,3) ^ ^-jVarccoeh l/(2g) 

Vl-(2/3)^ Vl-(2/3)2 Vl-(2/3)2 



v.(2) ''I 



— 2i/3 cos 



- sin 



16/3* 
3iV 



nTT 
N 
2 nn 
N 



/3-,0 



sm 



4 TiTT 



+ 0(4^) 

for n = 



iv 

+ o(| 

1,2,... 



r) 

N - 



-2i/JCOS ]v~ 



4_Ml 

+o(^) 

for n = 



2 (2n+l)7r 

JV 

■ 4 (2n + l)-^ 
JV 

0,1,..., f 



.(2) 



''=°-2i/3cos2f^ 

+ JV /™ JV 
I 32/3* -4 2n7r 
+ 3]V„'5^'^ iV , 

+0(^) + 0(f^) 
for n = 1,2,. . . , ^ 



^(3) = -2j/3cos(^) m'^) = -2»/3cos( '2"+'>" ) 

forn = 1, 2, . . . , ^ - 1 forn = 0, 1, . . . , .|- - 1 
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Table 6.2: For wavenumbers q corresponding to /3 > |: List of the 
eigenvalues of the matrix (|6.44|) given by equation H6.1U9|) . These eigen- 
values are directly related to the Redfield superoperator eigenvalues by equa- 
tion mum . 



N odd 


N = AI and q even 


N = AI and q odd 




or Af = 4/ + 2 and q odd 


or N = AI + 2 and q even 



^(4) -2if3cos ^=°° -2i/3cos ^ f.^ " =°° -2i/3cos 

2 ^1 I 1 I 2 /-. I 1 



+ iv(^+ o (".i+(l/2))x +n(^+ 12/32 sin2 



12/3^ 



+o(^) + o(^) 

for riniin <n < rimax 



+0(^) + 0(^) +o(^) + 0(^) +0(^) + 0(^) 

for njnin < n < "max for 

^min ^ 71 <C "^max 



If > |^,|: 



^(2) -2i/3cos 2f ^'=° -2if3cos i2!i+l)2L ^(2) _2i/3cos 2^ 

+ jv sin ^ + iV sin -h sin 

,16^" „■ 4 , 32/3'» .4 (2n+l)7r ,32/3^ .4 2n7r 

+ 3JV/'" iV + 3iV Sin jv + 3JV/™ JV 

+0(^) + 0(f^) +o(£^) + 0(f^) +o(i?^) + 0(f^) 

for n = 1, 2, . . . , nmin for n = 0, 1, . . . , nmin for n = 1, 2, . . . , nmin 

and n = rimax, - 1 and n = rimax, -y ^ 1 and n = rimax, -y ^ 



= -2if3cos{^) = -2j/3cos{i2n+l)ZL) 



for n = 1, 2, . . . , - 1 for n. = 0, 1, . . . , ^ - 1 



= ± V(2/3)2 - 1 
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We start by studying the eigenvalues obtained by fixing the 
wavenumber q (even or odd) and varying 9. For given physical parame- 
ters {A, A, Q, N), fixing q is equivalent to fixing (3. 

For /3 < |, the analytical expressions of the eigenvalues which concern 
us are summarized in Table IHTTl Two families of eigenvalues (/U^^^ and ^*-^-*) 
enter in the discussion for N odd, and three families {fJ-^^\ s^nd /U^^^), 
for N even. The numerical eigenvalues are plotted in figures IHHT a) . IH^ a) . 
and I6.3r al and are in very good agreement with the analytical results. The 
sole diffusive eigenvalue /U^^^ has a real part and no imaginary part. The 
— 1 other eigenvalues, either belongs to the /i^^^ family for N odd or to 
the fj,^"^^ and fi^^^ families for N even. The eigenvalues fi^"^^ and fi^^^ have an 
imaginary part which extends from —2/? to 2/3 and they generate oscillations 
in the dynamics as we shall see in the following section. The real part of 
the fi^"^^ eigenvalues is small and tends to zero in the large N limit. The real 
part of the fi^^^ eigenvalues is always zero. 

For /3 > ^, the diffusive eigenvalue /i^^-* has disappeared after merging 
with the other eigenvalues and the situation is slightly more complicated. 
The situation for a moderate value of /? > | is depicted in figures IHTT bl. 
I6.2r b). and I6.3f b) while the analytical expressions of the eigenvalues are 
given in Table Since ^i^^^ no longer exists, we have the two families of 
eigenvalues fi^"^^ and /u^^^ if is odd, and the three families /i*-^-*, fi^'^^ and 
^(4) ]\[ even. Two regions of the spectrum have to be distinguished. The 
eigenvalues //^^^ exist in the region where > |/ic| while the eigenvalues 
H^^^ exist in the region where < \nc\- We observe that the extra family 
of eigenvalues //^^^ has appeared because of the collision with the diffusive 
eigenvalue We can see in figures FOT b). FT^ b). and EiStb) that the 

analytical results of Table IH?^ reproduce very well the eigenvalues obtained 
by numerical diagonalization in the two regions. Here again, the number of 
eigenvalues is equal to N for a given wavenumber q, the imaginary part of 
the eigenvalues extends from —2/? to 2/?, and the real parts of all eigenvalues 
tends to zero in the large limit. 

A special situation occurs when /3 > ^ is increased to large values. This 
situation is depicted in figures IHTlT c). lOT c). and I6.3r c). The situation is 
similar to the previous one but the region > |/ic| has disappeared so 
that the family of eigenvalues ^^"^^ corresponding to the expansion /3 ^ 
no longer exists. For A^ odd and for A^ even with q odd, these eigenvalues 
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are replaced by the eigenvalues ;u(4) ei genvalues. For N even and q even, 
we find the two eigenvalues fi^^^ beside the family of eigenvalues /i^^^ . The 
agreement between the analytical and numerical results is very good here 
also. As before, the imaginary part of all these eigenvalues extends from 
—2i(3 to 2i(3 and their real parts tends to zero in the large N limit. 

A global view of the complete spectrum of the N'^ eigenvalues of the 
Redfield superoperator is depicted in figure I^TH bv varying the wavenumber 
g in a third dimension. Here, we only consider for simplicity the case where 
N is odd. The relation between the wavenumber q and the parameter /3 
is given by equation (|6.42j) . The wavenumber q varies in the first Brillouin 
zone or, equivalently, in the interval < g < 27r. We see in figure lH^ a) that 
the diffusive eigenvalues ^^^^ exists for all the values of the wavenumber in 
the case < ^ which implies (3 < ^. However, if > i, the diffusive 

eigenvalue ^^^^ disappears as expected for some values of the wavenumber 
corresponding to (3 > ^. This situation is observed in figure IHIHb). 

For very large values of > ^, the diffusive branch of the spectrum 

is reduced to the sole eigenvalue at g^^^ = 0, as seen in figure IH^ cl. In this 
case, diffusion has disappeared from the spectrum which only contains eigen- 
values associated with damped oscillatory behavior. The diffusive branch 
can be supposed to have disappeared when its last nonzero eigenvalue disap- 
pears in equation (|6.7.S|) . Therefore the diffusive branch disappears when the 
value of (3 for the first nonzero eigenvalue corresponding to g = ^ is larger 
than the critical value /3c = |. This happens when the coupling parameter 
exceeds the critical value given by 




This disappearance of the diffusion branch can be observed in figure IHTlT c') . 
We notice that the diffusive branch always exists in the infinite-system limit 
(N — > oo) in which case Ac can be arbitrarily small. 
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Figure 6.1: Eigenvalue spectrum for = 21 and given q: (a) /? = 0.4, 
(b) (3 = 0.7, and (c) (3 = 10. /Xg denotes the exact eigenvalues obtained by 
numerical diagonalisation and /x*-*^ the eigenvalues of the different families 
given in Tables 16.11 and 16.21 
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Figure 6.2: Eigenvalue spectrum for = 20 and q even: (a) (3 = 0.4, (b) 
P = 0.7, and (c) (3 = 10. /Xe denotes the exact eigenvalues obtained by 
numerical diagonalisation and /x^*) the eigenvalues of the different families 
given in Tables W/H and 
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Figure 6.3: Eigenvalue spectrum for = 20 and q odd: (a) /3 = 0.4, (b) 
f3 = 0.7, and (c) f3 = 10. fie denotes the exact eigenvalues obtained by 
numerical diagonalisation and fj,^^^ the eigenvalues of the different families 
given in Tables W/H and IH?^ 
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Figure 6.4: Complete spectrum of the eigenvalues related to the eigen- 
values of the Redfield superoperator by equation (|6.11fl|) . The system size 
is here = 21. In (a), M = 0.4. In (b), ^ = 0.7. In (c), M = 10. 
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6.4 Dynamical analysis 

We have described in the previous section the complete spectrum of the 
Redfield superoperator. We now need to discuss the dynamical implications 
of our spectral analysis. 

6.4.1 Theoretical predictions 

The spectrum of the Redfield superoperator determines the full dynamics of 
the subsystem interacting with its environment. In fact, projecting equation 
(|6.19|1 in the site basis and using the Bloch property (|6.28|) . we can write 



We see that the imaginary part of the eigenvalues Sqg will generate oscilla- 
tions in the dynamics and that the real part of these eigenvalues will generate 
a damping. The contribution to the long-time dynamics from the modes of 
the Redfield superoperator corresponding to the eigenvalues Sqg which have 
a large real part is very small. Only the modes corresponding to eigenvalues 
with a small real part will significantly contribute to the long-time dynamics 
and therefore to the transport property of the system. Here, we want mainly 
to focus on the long-time dynamics. 

As pointed out in subsection I6.,'17| the subsystem evolves in a non- 
diffusive regime if the chain is small enough A < Ac. In this regime, there is 
no diffusive eigenvalue so that the long-time dynamics is dominated by the 
eigenvalues /x^^^ describing the decay of the quantum coherences. The lowest 
relaxation rate is of the order of magnitude of the real part of the eigenvalues 
belonging to the /x*^^^ family of Table Using equation 1)6.110(1 . we find 
that this rate is given by 



We notice that these eigenvalues have a non- vanishing imaginary part so that 
the slowest modes relax exponentially but with oscillations. This means that 
the oscillating modes relax on the same time scales as the non-oscillating 
ones. These damped oscillations are reminiscent of a similar behavior in the 




(6.112) 




(6.113) 
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spin-boson model • We can thus interpret these damped oscillations as 
due to the damping of the quantum coherences. 

For a sufficiently large chain A > Ac, these damped oscillations disappear 
because the eigenvalue of the Redfield superoperator with the smaller real 
part is now the diffusive eigenvalue fi^^\ In this case, the lowest relaxation 
rate is given according to equation (|6.74|) by 




because the wavenumber of the slowest nontrivial mode takes the value 
q = 2-it/N and is inversely proportional to the length N of the chain. The 
signature of the diffusive regime is therefore that the relaxation rate scales 
as N-^. 



This behavior is consistent with the famous Einstein relation for the 
definition of the diffusion coefficient. Indeed, we can define the position 
operator of the subsystem x by 

{l\x\l') = I 5w . (6.115) 

The variance of the position is therefore given by 

N-l 

{x^){t) = Y.l^pii{t), (6.116) 

1=0 

where, using Eq (|6.112|) . 

/,„(*) = 5] a,^^(0)e^^«*+*'^ (6.117) 
q,e 

with the newly defined coefficients ag6)(0) = Cg5i(0)/3QQ. In the long-time limit 
and for a sufficiently large chain, only the diffusive branch s^^-* of eigenvalue 
spectrum will significantly contribute to the dynamics and one can write 



Pu{t) *=°° j;a,,(0) e«'^«*+^'^ 



(6.118) 
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where Sgo = -Dq^. Using qj = j^, Ag = q^+i - qj = ^ and f{qj) = 
Nag^e{0), we get 



Pu{t) ^ E /(^.) , (6.119) 

i=o 

which becomes in the hmit of an arbitrarily large chain 

Puit) ^' = °^ ^ / ^ /(^) • (6.120) 

If we choose an initial condition of the density function centered on a given 
site, f{q) = 1, and 



Finally, we recover the well-known Einstein relation for the diffusion coeffi- 
cient 

(x2)(t)^' = ~2L>t . (6.122) 

This demonstration is also a justification for calling fj,^^^ or s^^^ the diffusive 
eigenvalue. 

6.4.2 Numerical results 

The previous theoretical predictions are confirmed by the numerical inte- 
gration of the Redfield equation. We always take h= 1. 



Let us define the quantity we computed. Equation (|6.112|) can be written 



as 



P//(i) = E^'"'*^'''"'?e(0)' (6.123) 
q,e 

where aqg{0) = Cgg{0)p'QQ. If we introduce 

N-l 

p{q',t) = J2 pu{t) e-'^'^' , (6.124) 
1=0 
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using equation (|6.5() . we find that 



piq',t) = iV^a,,,(0) 



(6.125) 



e 



The evolving quantity we have computed is 



X{t)\ = \p{qut)\ 



(6.126) 



where we recall that qi = 2ti/N . We consider that the subsystem has initial 
conditions given by 



This means that X^e^g^l^) = 1/A^. The quantity X{t) is initially equal to 
unity (^(0) = 1) and tends to zero after a long time {X{oo) = 0) in an 
exponential way determined by the smallest non-vanishing real part of the 
eigenvalues Sq^e. 

Figures EUfa) and I6.5r b) show the same relaxation curves of |^(t)| but 
for two different time scalings. The different curves correspond to different 
values of the coupling parameter A, but to the same chain size N . Figure 
I6.5r al has the scaling A^t, which is characteristic of the non-diffusive regime 
and figure IH^T b) has the scaling t/{XN)'^, which is characteristic of the dif- 
fusive regime. Therefore, the curves in figure I6.5f a') with the same time 
scaling correspond to values of the coupling parameter such that A < Ac 
(in the non-diffusive regime) and they relax with oscillations at rates given 
by equation H6.113() . On the other hand, the curves of figure 1^31 a) with 
the same time scaling correspond to values of the coupling parameter such 
that A > Ac (in the diffusive regime) and they relax exponentially without 
oscillations at a rate given by equation (|6.114|) . 

Figure l6^ shows the scaling t/{XN)'^ in the diffusive regime for different 
curves at the same value of the coupling parameter but for different chain 
sizes N. This figure numerically demonstrates the property of diffusion in 
our model. 

An important remark is that the time scaling of the non-diffusive regime 
is naturally expected from perturbation theory. Indeed, according to per- 
turbation theory, the relaxation rates are proportional to the square of the 
coupling constant A. We recall that the Redfield quantum master equation 



Pw (0) = Sill Sio . 



(6.127) 
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1 1 (XNf 

Figure 6.5: Representation of the transition between the perturbative and 
diffusive regimes for a chain of size = 20 and parameters Eq = 0, A = 0.1, 
and Q = 1. The transition take place at the critical value Ac = 0.177 given 
by equation (|6.111|) . (a) The quantity versus the rescaled time A^t in 

order to identify the perturbative regime, (b) The quantity versus the 

rescaled time t/{XN)'^ in order to identify the diffusive regime. As predicted, 
we are in the perturbative regime if A < Ac and in the diffusive regime if 
A > Ac- The quantity is defined in equation (|6.126|) . 



6.4 Dynamical analysis 167 




1 1 (kNf 

Figure 6.6: The quantity defined by equation (|6.126|) versus the 

rescaled time t/{XN)'^ for various sizes = 20 — 50 of the chain. The 
parameters are Eq = 0, A = 0.1, Q = 1, and A = 0.3. Here, the system is 
always in the diffusive regime A > Ac. 

is obtained by second-order perturbation theory in the couphng parame- 
ter from the complete von Neumann equation for the subsystem interacting 
with its environment. It is quite remarkable that the time scaling of the 
relaxation in the diffusive regime is completely different since the relaxation 
rate depends on the inverse of the square of the coupling constant. 

In figure FTTI to figure IfHH we have represented the evolution of the prob- 
ability P = pii{t) to be on a given site of the subsystem for tree different 
values of the coupling parameter. In figure UTTl the coupling parameter is 
zero and the subsystem evolution is Hamiltonian. We can notice the inter- 
ferences showing the coherent character of the dynamic. In figure 16. 8| the 
coupling parameter is such that the subsystem relaxes in the non-diffusive 
regime. The dynamic still has interference patterns, but it is damped by the 
dissipation. Let's notice too that the center of the wave packet travels along 
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the subsystem. Finally, in figure 16.81 the coupling parameter is now such 
that the subsystem relaxes in the diffusive regime. The subsystem dynamic 
is no longer coherent and the interference pattern has disappeared. We also 
notice that the center of the wave packet no longer moves but that its width 
spreads diffusively. 

6.5 Infinite chain 
6.5.1 Spectrum 

The spectrum of the infinite chain coupled to its environment can be ob- 
tained from the spectrum of the finite chain in the infinite size limit N ^ oo. 
The wavenumber q becomes a continuous parameter varying in the first Bril- 
louin zone — vr < q < tt. 

For given wavenumber, the diffusive eigenvalue ^u^^^ or s^^^ given by equa- 
tion 1)6. 73|) remains isolated. Consequently, we obtain the result by which the 
dispersion relation of diffusion is exactly given by the analytical expression 















u 


\ 2 
■ 1\ 

sm - — 


= Dq' + 0{q^) 


.(6.128) 


The diffusion coefficient 




D 


~ QA2 ' 




(6.129) 



is proportional to the square of the parameter A of the tight-binding Hamil- 
tonian and inversely proportional to the parameter QX^ of the coupling to 
the environment. The transport is therefore due to the tunneling from site 
to site, which is hindered by the environmental fluctuations proportional to 
QX^. The eigenvalue (|6.128B can be considered as a Liouvillian resonance 
|44l I45j similar to the Pollicott-Ruelle resonances describing diffusion in clas- 
sical systems |29] . 
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Figure 6.7: Representation of the evolution of the probability P = pii{t) 
to be on a given site of the system. The initial condition corresponds to a 
Gaussian wave packet of standard deviation 2 centered on the site / = 25. 
We always have: = 50, Eq = 0, A = 0.1, Q = 1 and therefore Ac = 
0.1121. The system is isolated (A = 0) and the evolution is therefore a pure 
Hamiltonian multiperiodic dynamics. The center of the wave packet travels 
ballistically around the system. 
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Figure 6.8: Representation of the evolution of the probabihty P = pii{t) to 
be on a given site of the system. The situation is the same as in figure IFTfl 
except that the system is now interacting with its environment (A = 0.05). 
The system therefore relaxes in the non-diffusive regime (A < Ac) and the 
evolution is a damped multiperiodic dynamics. However, the center of the 
wave packet keeps travelling ballistically around the system. 
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Figure 6.9: Representation of the evolution of the probability P = pii{t) 
to be on a given site of the system. The situation is the same as in figure 
16.81 but the coupling between the system and its environment is stronger 
(A = 0.15). The system now relaxes in the diffusive regime (A > Ac). The 
multiperiodicity has disappeared, the center of the wave packet does not 
move and a diffusive spreading of the wave packet can be observed. 
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For 2 Ah < QX^, the diffusive eigenvalue exists for all the values of the 
wavenumber — vr < g < vr, as seen in figure UTTnT a). 

For 2 Ah > QX^, the diffusive eigenvalue only exists for all the values of 
the wavenumber in the range 

- Qc < q < +qc , with qc = 2arcsin . (6.130) 
as seen in figure E^tb). 

Beside the isolated diffusive eigenvalue, the spectrum at given wavenum- 
ber q contains a continuous part obtained by the accumulation of the eigen- 
values fJ-^^\ and fi^^^ in the limit N — > oo. Indeed, in this limit, 
all these eigenvalues accumulate into a segment of straight line given by 

AA f . q^ 2QA2 
Sqe = (.^^"^ 2/ °^ h?~ ' (6.131) 

with < 6* < vr and -vr < g < vr [see equations ()6.1U9() and H6.11U|) ]. 
This part of the spectrum is also depicted in figure 16.101 and describes the 
time evolution of the quantum coherences which are damped at the ex- 
ponential rate —TZe Sqo = 2Q}? /ti^ with possible oscillations due to their 
non- vanishing imaginary part Im Sqg. 

6.5.2 Temperature dependence of the diffusion coefficient 

In order to obtain the temperature dependence of the relaxation rates and 
the diffusion coefficient, we need to specify the coupling of the subsystem to 
its environment. In general, the coupling of a subsystem to its environment 
is described by the time correlation function. This correlation function is 
determined in terms of the spectral strength J{uj) according to 

a{t) = / duj J (lo) { h coth — cosut — i hsinujt ] . (6.132) 

Jo V 2fcBT J 

The Fourier transform of the environment correlation function, which is 
given by 

a(cj) = / — e*"*a(t) , (6.133) 
J~oo 2vr 

is therefore related to the spectral strength by 

d(a;) = (ncoth ^ + ^) , (6.134) 
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Figure 6.10: Spectrum of the infinite chain coupled to its environment: (a) 
in the regime 2 Ah < QX^ for Ah = 0.4 and QX^ = 1; (b) in the regime 
2Ah > Q}? for Ah = 0.6 and Q)? = 1 where the diffusive branch is hmited 
to the wavenumbers \q\ < Qc = 1.97022. 
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where the upper and lower signs respectively hold for u; > and u; < 0. 
Using equation (|A.7|) . and (|;-i.H2j) . one can notice that the spectral 

strength of the environment is proportional to the imaginary part of the 
environment susceptibility 

J{uj) = 2Im[x{uj)]. (6.135) 

Empirical forms for the spectral strength are often used in the literature 
inSl in21 such as 

J{uj) = Kuj'^e'^ , (6.136) 

where coc is a cutoff frequency and K a constant. The environment is called 
ohmic if 7 = 1, subohmic if 7 < 1, and superohmic if 7 > 1. Such empirical 
forms for the spectral strength are used to specify the dynamics of the 
subsystem over times longer than l/ujc which is of the order of magnitude 
of the time scale of the environment. Therefore, such spectral strengths 
are of primary importance to characterize the dependence of a{cv) at small 
frequencies. In the zero- frequency limit, the quantity 

lim aiuj) = lim Jiuj)^— = lim KksTuj"'-^, (6.137) 

(.01— +0 u)^0 UJ uj^O 

is well defined and non- vanishing only in the ohmic case 7 = 1, which we 
consider here. 



In our model [see equation ()6.15() ]. the environmental correlation function 
is given by a Dirac delta distribution a{t) = 2Q5{t). Therefore, the Fourier 
transform of the environmental correlation function is constant 

a{uj) = — . (6.138) 
vr 

This assumption is acceptable if the dynamics of the environment is much 
faster than the dynamics of the subsystem. Therefore, using equations 
(|6.137|) and (|6.138j) . the following identificat ion can be established 

Q = TrKkBT, (6.139) 

and the diffusion coefficient can now be written with an explicit temperature 
dependence 
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The conductivity a can be obtained using the Einstein relation 

where e is the electric charge of the carriers, n their density, and ^ their 
chemical potential. Assuming a low density of carriers, the chemical poten- 
tial depends on the density according to 

^ = /(T) + /cfiT Inn , (6.142) 

so that the conductivity is given by 

Therefore, the conductivity decreases as T"^ with the temperature. We 
notice that similar dependences are also obtained for the electric conductiv- 
ity in conducting polymers HO] and in Fermi liquids |59( I72j . This inverse 
power law is due to the existence of a single conduction band in our model. 
The transport is therefore confined in this band and no thermally activated 
transport process occurs in this model. Our diffusive transport phenomenon 
can be viewed as the result of the tunneling of the particle through the po- 
tential barriers of the system. This tunneling is more and more affected as 
the temperature increases. 



We expect that a crossover would occur to a regime where the diffusion 
coefficient has a temperature dependence of Arrhenius type if one consid- 
ered a system with at least two conduction bands. In such systems, we 
should find that a crossover between the quantum tunneling regime and 
a thermally activated transport regime when ksT becomes comparable to 
the energy spacing between the two bands ^1 EHl EHl IZS EHl HSj- From 
this viewpoint, the present model with a single energy band describes low- 
temperature behaviors. 



6.6 Summary 

In the present chapter, we have defined a simple translationally invariant 
subsystem which interacts with its environment by correlation functions 
which are delta-correlated in space and time. The reduced dynamics of the 
subsystem is described by a Redfield quantum master equation which takes, 
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for such environments, a Lindblad form. Thanks to the invariance under spa- 
tial translations, we could apply the Bloch theorem to the subsystem density 
matrix. In this way, we succeeded in getting analytical expressions for all 
the eigenvalues of the Redfield superoperator. These eigenvalues control the 
time evolution of the subsystem and its relaxation to the thermodynamic 
equilibrium. Two kinds of eigenvalues were obtained: the isolated eigenvalue 
()6.128|) giving the dispersion relation of diffusion along the one-dimensional 
subsystem and the other eigenvalues (|6.131|) . which describe the decay of 
the quantum coherences, i.e., the process of decoherence in this subsystem. 
The properties of the subsystem depend on the length N of the one-dimensional 
chain, on the width 2A of the energy band of the unperturbed tight-binding 
Hamiltonian and on the intensity Q of the environmental noise multiplied 
in the combination QX^ with the square of the coupling parameter A of per- 
turbation theory. 

We discovered that for a finite chain, there are two regimes depending on 
the chain length N and the physical parameters A and QX^. 
For a finite and sufficiently small chain, there is a non-diffusive regime char- 
acterized by a time evolution with oscillations damped by decay rates pro- 
portional to QX^. The oscillations are the time evolution of the quantum 
coherences. This non-diffusive regime exists if the coupling parameter is 
smaller than a critical value which is inversely proportional to the square 



For larger chains, we are in the diffusive regime with a monotonic decay on 
long times at a rate controlled by the diffusion coefficient. In this regime, 
the slower relaxation mode relaxes exponentially in time with the scaling 



In the limit of an infinite chain iV — ^ oo and for non-vanishing coupling pa- 
rameter QX^, the non-diffusive regime disappears and the subsystem always 
diffuses. 

The diffusion coefficient is proportional to the square of the width 2A of 
the energy band and inversely proportional to the intensity QX^ of the en- 
vironmental noise. Accordingly, we are in the presence of a mechanism of 
diffusion in which the quantum tunneling of the particle from site to site is 
perturbed by the environmental fluctuations. For an ohmic coupling to the 
environment, the diffusion coefficient is inversely proportional to the tem- 
perature. By using the Einstein relation between the diffusion coefficient 
and the conductivity, this latter is inversely proportional to the square of 
the temperature. 

Because the present model is very simple and shares some basic features of 
quantum diffusion and conductivity, we believe that it can easily be used to 




t/iXN)\ 
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understand how quantum transport by diffusion can emerge out of quantum 
coherent behaviors as in a molecular wire or in a carbon nanotube in the 
form of a loop which becomes larger and larger. 
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Chapter 7 

Conclusions 



The main goal of this thesis has been to understand the emergence of kinetic 
processes such as transport, decoherence and relaxation to equilibrium in 
quantum nanosystems. The understanding of these irreversible processes is 
of fundamental importance for the potential technological applications in 
the new field of nanosciences. 



7.1 Summary of the main results 



In the two first chapters of this thesis, we have briefiy reviewed the im- 
portant results already known in the understanding of irreversible quantum 
processes. In chapter |2l we have presented the important general concepts 
(quantum statistical ensemble, projected dynamics, Markovian evolution) 
necessary to understand how an irreversible evolution can arise from the 
reversible von Neumann equation. In chapter 13 we have reviewed the well- 
known and widely used theories of quantum irreversibility (linear response 
theory, Pauli equation, Redfield equation) by deriving them from the same 
perturbative expansion of the von Neumann equation in order to show their 
similarities and differences. In particular, we pointed out the link that exists 
between the Fourier transform of the correlation functions (appearing in the 
linear response theory and the Redfield theory) and the Fermi golden rule 
(that is a central quantity for the Pauli equation). 
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7. Conclusions 



7.1.1 Emergence of relaxation and decoherence in nanosys- 
tems 

In the traditional approaches hke the Redfield theory, the study of the irre- 
versible dynamics of a small quantum system interacting with an environ- 
ment relies on the assumption of an infinitely large environment. Assuming 
that the environment is infinite means that it is not affected by the in- 
teraction with the small quantum system, and, therefore, does not evolve 
and stays in its initial state. Of course this hypothesis greatly simplifies 
the calculation but is not always valid, especially in the new context of nan- 
otechnologies. Indeed, there are nanosystems in which a subsystem interacts 
with other degrees of freedom forming a finite environment. In such cases, 
the environment is not sufficiently large to remain unaffected by the inter- 
action with the small quantum subsystem. In this respect, an extension of 
the Redfield theory to finite environment, where the subsystem affects the 
environment and of course vice versa, is necessary. This first major achieve- 
ment of this thesis is presented in chapter 01 

The central quantity of our new theory is the subsystem reduced density 
matrix, but distributed in the environment energy. This quantity is ob- 
tained by applying a projection operator acting on the full density matrix of 
the total system and which depends on the energy of the environment. Let 
us notice that the projection operator used in the standard Redfield theory 
simply leads to the reduced density matrix of the subsystem without keeping 
information on the state of the environment. Thanks to our projection oper- 
ator, it is possible to monitor how the subsystem populations and coherences 
are distributed over the energy of the environment. Our equation correctly 
takes into account the exchanges of energy between the subsystem and the 
environment satisfying the total system energy conservation, which is not 
the case if the Redfield equation is applied to a finite environment. Our new 
theory is derived from a second-order perturbative expansion of the total 
von Neumann equation and is therefore valid at weak couplings between the 
subsystem and the environment. Another condition is that the initial condi- 
tion of the environment has to be an invariant state (i.e. it has to commute 
with the environment Hamiltonian). Furthermore, only our equation gives a 
consistent understanding, based on detailed balance, of the way a canonical 
subsystem equilibrium distribution (at a temperature given by the micro- 
canonical temperature of the environment) can arise when the initial state of 
the environment is a microcanonical distribution and when the environment 
has a large heat capacity. This indicates that applying the Redfield theory 
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to systems where the environment is in an initial microcanonical state is not 
vahd and that only our new theory is correct to describe kinetic processes 
within the microcanonical ensemble. When the environment becomes very 
large compared to the subsystem, we showed that, as expected, our new 
equation reduces to the Redfield equation. 

We have applied our new theory to a two-level subsystem interacting in 
a non-diagonal way with a general environment. Assuming that the en- 
vironment coupling operator in the environment eigenbasis representation 
vanishes, we have been able to obtain analytical expressions for the bi- 
exponential relaxation of the subsystem population toward their equilibrium 
value and for the oscillating relaxation of the subsystem coherences toward 
zero. This assumption is, in particular, exact in models where the environ- 
ment coupling operators are given by random matrices thanks to the average 
over the random-matrix ensemble. 

The spin-GORM model studied in chapter [21 of this thesis is a particu- 
lar case of such models where the environment operators (the environment 
Hamiltonian as well as the environment coupling operator) are given by ran- 
dom matrices taken in this case from the Gaussian orthogonal ensemble. 
In the first part of our study of this model, we focused on various spectral 
quantities (density of states, mean level spacing, SOE, ATPK). We classified 
the parameter space in three qualitatively distinct domains and we identi- 
fied a crossover at which the various spectral quantities change qualitatively. 
This crossover occurs at the value of the coupling parameter (between the 
spin and the environment) at which the interaction between the different 
coupled states becomes of the order of the mean-level spacing between these 
states. For values of the coupling parameter below this critical point, the 
spectral quantities are almost not affected by the coupling and look basi- 
cally as if there was no coupling. But for values of the coupling parameter 
beyond this crossover, qualitative changes in the various spectral quantities 
occur. It should be noticed that when the spectrum becomes continuous, 
the mean-level spacing and the crossover, both go to zero. 
The second part of the study of the spin-GORM model focuses on the sub- 
system (spin) dynamics. We identified three characteristic regimes of the 
spin relaxation corresponding to the three parameter domains identified in 
the spectral study. The first domain corresponds to weak couplings (be- 
tween the spin and the environment) and to situations where the energy 
splitting between the two states of the spin is small compared to the energy 
width of the environment spectrum. In this domain, the relaxation is expo- 
nential and scales in time as X^t. The second domain also corresponds to 
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weak couplings, but the energy splitting between the two states of the spin 
is now large compared to the energy width of the environment spectrum. In 
this domain, the relaxation is a purely non-Markovian dynamics (typically 
rapid oscillations at the spin frequency damped on a time scale of the order 
of the environment correlation function time scale) and scales in time as t. 
The last domain corresponds to very large couplings (larger than the spin 
energy splitting and than the environment spectral width). In this domain, 
the spin relaxation scales like Xt and is governed by a Bessel function of the 
first kind Ji{Xt)/{\t). 

The study of the weak-coupling regime of this model was of special interest 
for us because it allowed a detailed comparison between the spin dynamics 
predicted by our theory, the Redfield theory, and the exact numerical re- 
sults computed with the von Neumann equation for the total system. On 
one hand, this study confirmed the necessity of extending the Redfield the- 
ory to our new theory when the typical spin energy becomes non-negligible 
compared to the typical energy scale of variation of the environment den- 
sity of states or, equivalently, when the heat capacity of the environment 
multiplied by the micro canonical temperature of the environment becomes 
of the order or smaller than the energy splitting between the two levels of 
the spin. On the other hand, this study showed that both theories predict 
the same dynamics on non-Markovian short time scales. Finally, and this is 
the second central result of this thesis, this study allowed us to find a very 
important criterion for the validity of kinetic equations in finite quantum 
systems. In fact, for values of the coupling parameter A between the spin 
and the environment below the crossover that we identified in the spectral 
study, the exact dynamics does not follow the prediction of our kinetic the- 
ory. The crossover take place at the critical value Ac = 0(1/ ^/N), where N 
is the number of energy levels in the environment. It is only beyond this 
crossover that the kinetic theory applies. As expected, the critical value goes 
to zero in the thermodynamic limit and, consequently, the kinetic theory is 
always valid. Let us notice however that the coupling parameter cannot 
be too large, otherwise the weak-coupling assumption at the origin of the 
kinetic theory fails. It is conjectured that the existence of such a crossover 
to see kinetic processes is a general result and that kinetic processes only 
occur in finite quantum systems when the typical interaction terms between 
the coupled levels becomes larger than the mean level spacing between these 
coupled levels. The interaction has to be large enough to somehow "mix" 
the coupled levels, so that the unperturbed structure of the spectrum dis- 
appears. 

This conclusion was strengthened by noting the fact that a self-averaging 
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property arises for values of the coupling parameter beyond the crossover. 
In fact, if the initial condition of the total system is chosen to be a pure state 
at a given energy of the environment, when the coupling parameter is below 
the crossover, the subsystem dynamics is strongly oscillating and does not 
follow the microcanonically averaged dynamics at this environment energy. 
However, when the coupling parameter is above the crossover, the pure state 
dynamics starts to be smoother and to follow the microcanonically averaged 
one. This self-averaging behavior is reminiscent of classical ergodicity and it 
is a very important point to understand the emergence of kinetic processes 
in quantum systems. 

In the last part of the spin-GORM model study, we focused on the equilib- 
rium properties of the subsystem (spin), but also of the total system. Once 
again, we found that the equilibrium spin distributions expected from our 
kinetic theory are valid if the values of the coupling parameter are beyond 
the crossover. We also showed that when the heat capacity of the envi- 
ronment is large, the subsystem thermalizes at a temperature given by the 
microcanonical temperature of the environment. We finally showed that the 
complexity of internal couplings (given by random matrices) can have the 
effect of thermalizing the total system without the need of an external heat 
bath. 

7.1.2 Transport in nanosy stems 

Finally, in chapter El we have defined a simple model in which a transla- 
tionally invariant subsystem interacts with an environment by a coupling in 
terms of correlation functions which are delta-correlated in space and time. 
The subsystem is a one-dimensional tight-binding subsystem made of 
sites and possessing a single energy band of width 2A. The intensity of the 
delta-correlated environmental noise is given by Q and the intensity of the 
subsystem-environment interaction by A. For weak coupling, the reduced 
subsystem dynamics is described by a Redfield quantum master equation 
which takes, for such correlation functions, a Lindblad form. Thanks to the 
Bloch theorem, we succeeded in getting analytical expressions for all the 
eigenvalues of the Redfield superoperator. By studying the Redfield super- 
operator spectrum we discovered that for a given system (^4 and given) 
interacting with a given environment {Q given), there are two subsystem 
relaxation regimes, depending on the intensity of the interaction between 
the subsystem and its environment (A). These two relaxation regimes are 
separated from each other by a critical value of the coupling parameter 
Ac = ^/2AhTr/{QN). In the diffusive regime corresponding to A > Ac, a dif- 
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fusive branch exists in the Redfield superoperator spectrum which contains 
at least one nonzero eigenvalue. In this regime, the slower relaxation mode 
relaxes exponentially in time with a scaling t/{\N)'^. In the non-diffusive 
regime corresponding to A < Ac, the last nonzero eigenvalue of the diffusive 
branch has disappeared. The slower relaxation mode of this regime relaxes 
exponentially in time but with oscillations and with a scaling X^t. 
In the limit of an infinitely large number of sites on the subsystem N oo, 
and for the nonzero value of the coupling parameter, the non-diffusive regime 
disappears (Ac — > 0) and the subsystem always diffuses. We showed that 
the existence of the diffusive branch in the Redfield superoperator spectrum 
implies an Einstein relation for the diffusion coefficient and that the inverse 
temperature dependence of the diffusion coefficient can be established if 
one considers our model of environment as the infinitely fast environment 
limit (fast compared to the subsystem dynamics) of an ohmic environment. 
Demonstrating the existence of quantum diffusion on a spatially extended 
system constitutes the third central result of this thesis. 

7.2 Perspectives 

The study of quantum kinetic equations started a long time ago and pres- 
tigious names are associated with this issue. The earlier results where first 
obtained by Pauli jHS] in the late twenties with the derivation of the Pauli 
equation and later (in the fifties and sixties) by Van Hove 91 , Montroll (6f)j . 
Zwanzig |97j and others with more systematic derivations of the Pauli equa- 
tion. These early studies only focus on the subsystem population dynamics. 
Studies focusing on the full subsystem density matrix started at the end 
of the fifties with the first derivation of the Redfield equation by Redfield 
|77j and have continued until now with various systematic derivation of this 
equations by Cohen- Tannoudji ^Hl, Kubo jSJ, Gaspard IHOl and others. In 
these studies of the Redfield equation, it is an external environment which 
generates the irreversible dynamics on the subsystem and one has to assume 
that this environment is large enough not to be affected by the subsystem. 
Our kinetic equation is a new step in the theory of the quantum kinetic 
processes, because we further extend the previous theories to include in our 
description some variables of the environment (here the environment energy) 
in order to study the irreversible dynamics of a subsystem interacting with 
an environment which can be affected by the subsystem dynamics. The 
quantity on which our theory focuses is the density matrix of the subsystem 
distributed over the environment energy. 
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Our new approach to kinetic processes in quantum systems is certainly 
an important new tool to deepen our knowledge of the foundation of sta- 
tistical mechanics in quantum systems, especially because of the growing 
interest for the understanding of irreversible behavior such as decoherence 
in nanosystems. In these nanosystems, one often encounters "structured" 
environments (i.e. environments with a density of states which highly varies 
on short energy scales). As soon as one wishes to study the dynamics of a 
subsystem interacting with such environments, one has to take into account 
the effects of the changes in the energy of the environment caused by the 
subsystem if the subsystem energy is equal or larger than the typical energy 
scales of the environment density of states. This requirement is fulfilled in 
our new theory. 

Spintronics could provide situations in which models similar to the spin- 
GORM model could apply. In spintronics, one studies the spin of an electron 
moving in a nanodevice. If the energy difference between the two spin states 
of the electron (due for example to an applied magnetic field) is not neg- 
ligible compared to the typical energy scales of variation of the density of 
states of the nanodevice, again, a correct kinetic theory should take into 
account the effect of the spin on the nanodevice. The correct description of 
the electronic spin relaxation and decoherence in such devices has a huge 
importance for potential nanotechnologies such as quantum computers. 

It should be possible to extend our new theory in order to obtain a mean 
field equation for many-body systems. Such an extension could allow the 

description of the dynamics of an ensemble of interacting particles where 
energy conservation plays an important role. In fact, we have seen that 
our equation is the only kinetic equation giving an adequate description of 
the irreversible relaxation of a subsystem interacting with an environment 
described by a microcanonical state. Such kinds of extension could possibly 
lead for example to quantum kinetic equations for Bose-Einstein conden- 
sates. 

Our kinetic equation allows us to follow the dynamics of each element 
of the subsystem density matrix as a distribution in the energy of the en- 
vironment. Our theory could therefore play an important role in order to 
understand the quantum measurement within a pure quantum scheme where 
the measured system as well as the measurement apparatus are described by 
quantum mechanics. The density matrix of the measured system would be 
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distributed over the energy of the measurement apparatus allowing possible 
interesting interpretations. 

Stochastic Schrodinger equations describing the wave function dynamics 
of a subsystem interacting with its environment have been intensely studied 
during this last decade |^ ISHl • These equations will probably play an 
important role in the future to describe individual elements of a quantum 
statistical ensemble. It has been shown that if ensemble averaged, such 
equations obey the well-known master equations. Stochastic Schrodinger 
equations corresponding in this way to the Redfield equation were also ob- 
tained . An interesting issue would be to find the stochastic Schrodinger 
equations which, if ensemble averaged, gives our new kinetic equation. 

The results in the last chapter of this thesis open perspectives on trans- 
port in nanosystems. 

Investigating quantum transport in a translationally invariant model 
with at least two conduction bands is an important extension of the sin- 
gle band model studied in this thesis. In fact, in such a double band model, 
two kinds of transport are possible. At very low temperatures, only the lower 
conduction band is filled, and the transport occurs by the tunneling of the 
particle through the potential barriers. At high temperatures, the thermal 
fluctuations coming from the environment are sufficiently strong to allow 
the particle to jump from a conduction band to another and the transport 
occurs mainly by thermal hopping. This second type of transport becomes 
effective when the temperature creates energy fluctuations comparable to 
the energy gap between the conduction bands. At low temperatures in the 
tunneling regime, we have seen that the transport is slowed down when 
the temperature is increased. At high temperatures where the transport 
by hopping is dominant, one expects an Arrhenius temperature dependence 
of the transport coefficient. Therefore, at a given temperature, a crossover 
between the two transport regimes should exist. 

In our model of quantum diffusion we have used semi-qualitative argu- 
ments in order to assert that normal transport (transport with a well-defined 
transport coefficient) can only occur if the environment spectral strength is 
ohmic. Here we may wonder if superohmic or subohmic environments gener- 
ate abnormal transport on the subsystem, as suggested by possible extension 
of our ohmic model studied in chapter El and also by Refs. |58l l92]. 
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In this thesis we have encountered the criterion X^N > Cte twice for 
the emergence of irreversible processes in nanosystems. In the spin-GORM 
model, N is the number of states of the environment. In the model of 
quantum diffusion, N is the number of sites of the translationally invariant 
system. Despite the different meanings of N in both cases, it is worthwhile 
to make the analogy. One can wonder if such a criterion for having kinetic 
processes in finite quantum systems will replace in the future the much 
stronger criterion of the thermodynamic limit (which requires AT — > oo). 
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Appendix A 

Correlation functions of 
quantum systems 

The correlation function between two operators B and A of a quantum 
system with Hamiltonian H is defined as 

aAB{t) = Trp^'^A{t)B = Tr/5^«ei^*ie-i^*S, (A.l) 

where p^'^ is an invariant density matrix such that [H, p^^] = 0. The corre- 
lation function has the property a^Bit) = a*^^{—t). 
The Fourier transform of the correlation function is defined by 

/oo ^iuit 
dt—aAB{t). (A.2) 
-oo 

The Fourier transform of the correlation function has the property aAB{'^) = 

General properties 

The correlation function can be written as 

OiAB{t) = CAB{t) + iDAB{t), (A.3) 

where 

CAB{t) = \Trp^'i[A{t),BU = liaABit) + a*Mt)) (A.4) 

and 

DAB{t) = -'-Tvp-'i[A{t),B] = -\{aAB{t) - a\j,{t)). (A.5) 
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Notice that CAB{t) and DAB{t) are real and have the property 

Clsit) = CABit) , Dlsit) = DABit) (A.6) 

CABit) = CBAi-t) , DAB{t) = -DBA{-t). 

The Fourier transform of the correlation function can be written as 

otABiy^) = Cab{!^) + iDAB{(^) = a*BA{i^), (A.7) 

where 

CIbH = CABi-uj) , D*ABi^) = Dab{-uj) (A.8) 
C'asM = C'sa(-w) , -Dab(w) = -£>s^(-u;). 

If we define 

C%Bi^) = iCAB{io) + CBAiu;))/2 (A.9) 
C%s{u) = {Cab{uj)-Cba{co))/2 
D%s{^) = {bAB{oj) + bBA{oj))/2 

D\s{u) = {I)AB{^)-bBA{i^))/2, 

we can separate the real and the imaginary parts of the Fourier transform 
of the correlation function in the following way 

OiABd-o) = TZe[aABi'^)] + iIm[aABiuj)], (A. 10) 

where 

UeiaABH] = C%B{^)+iD%B{oj) (A.ll) 
Im[aAB{co)] = b'XB{'^)-iC%s{u;). 

The case of a discrete spectrum 

If the spectrum of H is discrete {H\n) = En\n)), the correlation function 
reads 

aAB{t) = p^?,et(^''-^"')U„„,B„,„, (A.12) 

nn' 

and the Fourier transform of the correlation function reads 

aAB{uj) = ^ pt^n^{uJ + Unn')Ain'Bn'n, (A.13) 
nn' 
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where hwnn' = En - En'- 



For the canonical case (f^ = e /Z where Z = Tre the correlation 
function reads 



nn 



and its Fourier transform 

aABiP,Uj) = 2^ ^ S{UJ +UJnn')Ann'Bn'n- (A.15) 
nn' 

For the microcanonical case where p^^ = S{e — H)/n{e) with n(e) = Tr5(e — 
Hq), the correlation function reads 

aAB{e,t) = ^^^^ei(^-^"')*A„„,B„,„, (A.16) 
and its Fourier transform 

aAB{e,Uj) = Y^^^—T^^6{u + Unn')Ann'Bn'n- (A.17) 

, n(e) 

nn' 

For the pure case p'^'^ = |n)(n|, the correlation function reads 

ocABiEn, t) = Y, e^^^"-''"'^'Ann'B„,n, (A.18) 

n' 

and its Fourier transform 

aAB{E„,L0) = ^ + UJnn')Ain'Bn'n- (A.19) 
n' 

The case of a quasi-continuous spectrum 

The spectrum of H is quasi-continuous if the discrete eigenvalues are so 
densely distributed in energy that one can use a continuous description for 
the energy {H\e) = e\e)) and that the operator of interest can be considered 
as a smooth function of the energy in the eigenbasis representation of the 
system. 
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In this case, taking the quasi-continuous hmit means that 

En ^ e (A.20) 
^ y (ien(e) 

Ann' ^(e,e') 
Bnn' — -S(e, e') 

where n(e) is the density of states of the system at energy e. 
The correlation function of a quasi-continuous system reads 

aAB{t) = J de j de"n{e')n{e")p^^e')A{e',e")B{e",e')eT^^''-'"^\ (A.21) 

and its Fourier transform 

OAsH = f^jde'j de"n{e')n{e")p^'i{e)S{fiaj + e' - e")A{e' ,e")B{e" ,e') 
= de'n{e')n{hio + €')p^'i{e')A{e' , ho + e')B{hu; + e' , e'){A.22) 

In the canonical case, the correlation function of a quasi-continuous system 
reads 

aAB{P,t) = de' de"n{e')n{e")^^A{e' ,e")B{e" ,e')eTi^''-'"^\ (A.23) 

and its Fourier transform 

aAB{0,co) = h / de'n{e')n{nu; + e')—^A{€',rko + e')B{nw + €',e').{A.2A) 

In the microcanonical case, the correlation function of a quasi-continuous 
system reads 

J J n{e) 
= J de'n{e')A(e,e')B{e',e)ei^'-''^\ (A.25) 

and its Fourier transform 

aAB(e,w) = h j de'n{e')A{e,€)B{e ,€)5{e-e' + nw) 

= hn{e + nw)A{e,e + hLv)B{e + lTW,e). (A.26) 
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We notice that 

-f3e 



aAB{f3,t) = / de^^ aAB{e,t), (A.27) 



and therefore that 

-/3£ 



aABW,uj) = de^^ aABie,uj). (A.28) 

KMS property 

The usual KMS (Kubo-Martin-Schwinger) property, which is easy to verify, 
is that the canonical correlation function satisfies 

aAB{l3;t) = aBA{P;-t-il3/h). (A.29) 

Taking the Fourier transform, on gets that 

dAB{f3;oj) = ef^f^d^eip'-u;) (A.30) 
= ef'^dBAW;-u;). 

This implies that 

Cab{I3; uj) = 2i^^DAB{P; co) (A.31) 
where we have defined 

£;^H = — coth^. (A.32) 

The microcanonical analogue of the canonical KMS property is also easy to 
verify and is given by the relation 

dAB{e,i^) = + ^) d\^{e + hw, -lu) (A.33) 
n(e) 

n{e + nw) ^ 

= r-\ — aBA{e + nw,-uj). 

n{e) 

In fact, by averaging ()A.33|) over energy with a canonical weight as in ()A.28|1 . 
one recovers (|A.30|) . However the KMS property is not as simple and there- 
fore as useful in the microcanonical case as in canonical case. 
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Appendix B 



Systems with a large heat 
capacity 



Here we consider an isolated system with a dense spectrum described by a 
density of states n(e). The statistical state of the system is described by a 
microcanonical distribution at energy e. For such a system, the microcanon- 
ical entropy is given by 

S{e) = hlnQ{e) =hlnn{e)6e, (B.l) 

where f^(e) is the number of states in the energy shell which is given by the 
density of states n(e) multiplied by the small width of the energy shell 5e. 
The microcanonical temperature is defined as follows 

' -"'^'^ (B^2) 



Tmic{e) de 
The microcanonical heat capacity is defined as 

1 dTmic{e) , . 

Using HB.1|) . we can connect the density of states of the system at two 
different energies as follows 

ra(e + A) _ e^i^+^Ukb^ 

n{e) ~ ~^SW^- ^ ^ ^ 

Now, we suppose that the energy A by which the energy of the system 
has been increased is small enough to perform an expansion of the entropy 
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S{e + A) in A around e. Using HB.2|) and (|B.3|) we can write 



A 1 A^ 

= S{e) + —^—--——^— — + .... (B.5) 

This means that if the heat capacity of the system is large enough to neglect 
the second order term in (IB.5I) 



a(e) » (B.6) 



(IB. 41) can be written as 



n(e) 

where (3mic{^) = ^/kbTmici^)- This shows that a Boltzmann distribution, 
characterising a canonical ensemble, can naturally emerge from a micro- 
canonical description. We also notice that the micro canonical temperature 
of a system satisfying (|B.6|) can be considered constant on energy scales of 
order A 

Trmc{e + A) = r„,,(6) + ^^^A+... (B.8) 

de 

= T_(e) + ^A + ... 



Appendix C 

Redfield equation for 
two-level subsystems 



In this appendix we apply the Redfield equation (|3.98|) to a two-level sub- 
system interacting in a non-diagonal way with its environment. The Hamil- 
tonian of the total system is given by (|4.55|) . The non-Mar kovian Redfield 
equation becomes 



{±\Mt)\±) = -^X{mT\ps{t)\±) - {Mps{t)\T)) (C.l) 

{ -{±\psit)\±)d{u;)iei^^^-'^> + e-i(±^-Mr) 



for the population and 



(±l^5(t)|T> = T^A(±|p5(t)|T)-^A(S)((T|p5(t)|T>-(±|/55(i)|±)) 

+ dT duj (C.2) 

{ -(±|p5WlT)a(u;)(et(±^-M- + ei(±^+'-)^) 
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for the coherences. Performing the Markovian approximation, we get 



{±\psit)\±) = -'-X{B){{T\ps{t)\±) - {±\psit)\T)) (C.3) 
~27r^a{±A/h){±\ps{t)\±) 

+27:^d{TA/h){T\ps{t)\T) 



for the population and 



{MMt)\T) = T^A{±\ps{t)\T) - ^X{mT\Mt)\T) - {±\Psit)\±)) 

-^{±\ps{t)\T){+7T{aiA/h) + ai-^m 
, A /■ , _ a(uj) 

+ ^(T|ps(«)l±){+f (SCA/fi) + 5(-A/S)) 



-4 / -"wfc' '^■^> 



for the coherences. Now we use the Bloch variables given by 

zit) = Trm^z = {+\Psit)\+) - {-\psit)\-) (C.5) 

= TYp{t)^, = {+\Psit)\-) + {-\ps{t)\+) 
y{t) = Tvp{t)^y = i{{+\Psit)\-)-{-\ps{t)\+))- 

Assuming as we did in section ^31 that {B) = ^, we get 

z{t) = 7(z(oo) - (C.6) 
xit) = -j-yit) 

y{t) = (^+r)x(t)-7y(t), 



^Notice that even if (B) 7^ 0, the order A term in the Redfield equation can vanish 
using the redefinition I3.10H and 13.1021 . 
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where 



7 = 27r^(a(A/n)+d(-A/;i)) (C.7) 



A2 - (hw)^ 



dj-A/h) - a{A/h) 
= a{-A/h) + a{A/hy ^^-^^ 

The solutions of this equation are given by 

z{t) = z(oo) + (z(0) - z(oo))e-^* (C.9) 
= f]^ Sin fj^(^.r)- (2)^1 e-i^ 



Vt(t + r)-(F 



+x(0)cosU/|(^ + r)-(|)2i e-i^ 



x(0)(4+r)-^(0)l . / /A A 7,2,^ -it 

y(t) = , ''^ ^ sm W-(- + r)-(-)^ne 2* 



+y(0)cos( \/|(| + r)-(|)2t e-i* 



It can finally be noticed that 

V 2 /•oo 



^ = / dr cos(At/;i) C(r) (C.IO) 

JO 

r = / dr sin(AT/;i) C(r) 

^2 /"OO 

72(00) = 4-2 / dr sin(Ar/^) D(t), 
Jo 

where we recall that a{t) = C{t) + iD{t) where C{t) and D{t) are real 
functions. 
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Appendix D 

Gaussian orthogonal random 
matrices (GORM) 



A Gaussian orthogonal random matrix (GORM) Y is characterized by M, 
the size of the matrix, and by the parameter Cy, which enters the Gaussian 

probabihty distribution P{Y) = Ce 2"^^(^ ) of the whole matrix. The 
statistical properties of a GORM are preserved under orthogonal transfor- 
mations. Because the matrix is symmetric, each non-diagonal element Yij 
is equal to its transposed Yji. The ^(^+^) independent matrix elements of 
Y are Gaussian random numbers of mean zero. The standard deviation of 
the non-diagonal matrix elements (tJjj^ and the standard deviation of the 
diagonal matrix element aj) are related to by 



^^4d = W— • (D-1) 



The density of states of the GORM Y is defined by 

M 

d{E)=Y,S{E-Ei), (D.2) 

i=l 

and the smoothed density of states by 

1 f^H 

d{E) = lim - / d{E) dE, (D.3) 

where e is a small energy interval which is large enough to contain many 
states in order for d{E) to be smooth. The averaged smoothed density of 
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D. Gaussian orthogonal random matrices 



states is an ensemble average of x realizations of the GORM. Such an en- 
semble is called Gaussian orthogonal ensemble (GOE). It is well known 
inn ESI in] that the ensemble averaged smoothed density of states {d{E))^ 
obey the Wigner semicircle law in the limit x ^ oo: 



{d{E))c 







if \E\ < ./2M 
•'II \/ a,> 



if 



\E\ > 



(D.4) 



The domain of energy where the eigenvalues are distributed (i.e., the width 
of the semi-circle) is TY = ^J'^-- Notice that when M — > oo, d{E) — > 

{d{E))oo, and therefore d{E) follows the semi-circle law. The following no- 
tation is used in the present paper: n'^{E) = {d{E))oo- 
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Perturbation theory for the 
spin-GORM model 

There is no analytical way of getting a general form of the eigenvalues 
of the total system, but the three terms in Eq. H5.6|) have different orders of 
magnitude, depending on the value of the parameters A and A. The system 
and the environment Hamiltonians are of order A and 1, respectively, while 
the coupling term is of order A. Therefore, we can examine the different 
extreme cases that can be treated perturbatively. 

A,l > A : 

When the system and the environment Hamiltonians are larger than 
the interaction term in Eq. H5.6() . we can treat the interaction term in a 
perturbative way, taking the system and the environment Hamiltonian as 
reference, 

Ho\sb) = E'^Jsb) , (E.l) 

where we replaced the index n by the two indices s, b. The perturbed energy 
is given to the second order by 

We notice that the first nonzero correction to the non-perturbed eigenstate 
is of order A^. 
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E. Perturbation theory for the spin-GORM model 



A > 1, A : 

When A is large compared to A and 1 in Eq. H5.6|) . it is possible to 
consider the interaction term as the reference Hamiltonian and to treat Hs 
and Hb as small perturbation. Transforming ()5.6() by a unitary matrix 
acting only on the system degree of freedom, we get 

H = + Hb + Xa,B. (E.3) 
The non-perturbed reference Hamiltonian is therefore, 

^0 = A<7.B. (E.4) 
Let Ef^fj and {kt]) = |K)(g)|r/) be respectively the eigenvalues and eigenvectors 
of Ho: 

HoIkv) = = XE^r]\Kr]) = XKEn\K.ri), (E.5) 

where t] = 1, ^ and k = ±1. The energy of the perturbed Hamiltonian 
is thus given to the second order perturbation in j by 

- = kE, + -{,\Hs\v) + ^ E j,o _ EO + ^ • ^^-'^ 



nr] K'rj' 



1 > A,A : 



In this case, the environment Hamiltonian is large compared to the sys- 
tem Hamiltonian and the interaction term so that they both can be consid- 
ered as perturbations. We get 

= E,, = ^s + E^ + X^Y: + OiA') + O(A^). (E.7) 

^ b'^b^b -Ey 

A » 1,A : 



We now suppose that the system Hamiltonian taken as reference is large 
compared to the environment Hamiltonian and the interaction term, so that 
these last two terms can be considered as perturbations. We then get 

E^ = E,,, = ^s + E^ + s^Y. K^'l^l^)!' + 0{l) + 0{X'). (E.8) 

Two additionals situations, 1, A ^ A and A, A ^ 1 , could be consid- 
ered but cannot be treated perturbatively because no reference basis exists 
in which Hb and B are simultaneously diagonal. 
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